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16238 SOIL FAILURE AND UNDRAINED CYCLIC LOADING 
| KEY WORDS: Anisotropy; Axial strain; Cyclic loads; Failure; Pore pressure; ue vag 
Sand; Shear stress; Soil dynamics; Soil porosity; Soil tests; Strain; Stress; dacs 

ABSTRACT: The failure modes of saturated sand in undrained cyclic loading are 
explored. A typical soil element subject to dynamic load exists initially under a 
nonhydrostatic state of stress. Anisotropically-consolidated samples in triaxial tests 
with a superimposed cyclic deviator stress represent such elements, except for the — 
rotation of principal stress directions during cycling. The response of these samples — 
depends on the degree of anisotropy and whether shear stress reversal occurs. For 
isotropically-consolidated undrained tests on saturated sand samples, the pore pressure 
increases with each cycle, thus decreasing the effective confining stress. The axial 
cyclic strain also increases as the effective stresses approach the failure state. However, 
the static axial strain should remain close to zero. For anisotropically- -consolidated 
undrained tests with no shear stress reversal, the pore pressure generally, but not _ 
always, increases. In this case, the static axial strain will increase with each _ and — 
the cyclic strain will remain smaller than the static strain. 
REFERENCE: Selig, Ernest and Chang, Ching S., “Soil Failure Modes 
Undrained Journal of the Geotechnical Division, ASCE, 


16241 THE STEADY STATE OF DEFORMATION 
KEY WORDS: Deformation; Deformation analysis; Flow characteristics: i 
Liquefaction; Particle shape; Residual shear strength; Sand; Shear stress; usage T es 
Steady state; Stress strain characteristics; Stress strain relations 
ABSTRACT: The steady state of deformation for any mass of punta is that state in 
which the mass continuously deforming at constant volume, constant normal 
effective stress, constant shear stress, and constant velocity. The steady state of 
deformation is achieved only after all particle orientation has reached a statistically a > 
steady-state condition and after all particle breakage, if any, is complete, so that the — iv 
shear stress needed to continue deformation and the velocity of deformation remain — 
constant. The similarities and differences between steady-state deformation and the — 
current use of the term critical state are described. A special undrained triaxial test on i 
a sand is presented to demonstrate clearly that a special flow structure exists during ; 
steady-state deformation, which is quite different from the initial structure, and which | : 
is credited to a nonegndom, ie., oriented, arrangement of the sand grains. 
REFERENCE: Posie, ‘Steve : J., “The Steady State of Deformation,” Journal of the a 
Geotechnical Engineering Division, ASCE, Vol. 107, No. GTS, Proc. Paper 16241, 
16234 CYCLIC STRENGTHS AND SAMPLING TECHNIQUES IBOLT tate 
KEY WORDS: Capillary pressure; Cyclic loads; Cyclic strength; Freezing; 
Liquefaction; Pitcher samplers; Sampling; Sands; Soil properties; Triaxial 
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ABSTRACT: Two distinctly different field sampling and 
procedures, used to obtain similar sand samples for cyclic testing, are documented. 
3 comparison of the cyclic strengths of the samples is ‘resented together with 
| examinations of optimum techniques in undisturbed sampling and sample preparation. 
t One procedure used included Pitcher barrel sampling, followed by laboratory trimming 
by mechanical extrusion into smaller diameter brass molds. The second procedure — 
included hand-carved block sampling followed ‘by either further hand-trimming, 
freezing and trimming, or saturation and trimming. The block- -sampling and hand- 
trimming preparation procedure resulted in significantly higher cyclic strengths 
-! corroborating the concept that disturbance effects can cause major reductions in the — 
¥ 
| REFERENCE: _ Espana, Carlos, Chaney, Ronald C. , and Duffy, Dennis, “Cyclic : 
Strengths Compared for Two Sampling Techniques,” Journal of the Geotechnical _ 
Engineering — ASCE, Vol. 107, Proc. Paper 16234, 1981, pp. 


wig 

| 
— 
q 


SS 


“ABSTRACT: The of equivalent uniform stress in soil ‘dynamics to 
| the study of soil behavior during and after an earthquake is explored. The actual i. 
amit time histories produced by an earthquake can be represented by uniform 
amplitude cyclic stresses, although there may be a considerable amount of uncertainty i 
associated with them. The stress level of 75 percent of the maximum is suggested for _ 
| such conversion, since in this case the uncertainty in the normalized soil-strength ¥ 
| curve has a minimum effect on the value of the Neg versus M relationship. A 
| relationship between Neg and the earthquake Tnagnitude is proposed here 
based on results available in the literature. The Neg could be estimated adequately by _ 
considering the component of excitation containing the peak acceleration. The Ne, 
versus M relationship proposed is somewhat different from the relationship suggested 
by Seed and Idriss. This discrepancy in Neq values may not yield "significant — 
REFERENCE: Haldar, and Tang, Wilson H., “Statistical Study of Uniform 
_ Cycles in Earthquakes,” Journal of the Geotechnical Engineering Division, ASCE, Vol. ; 
107, No. GTS, Proc. Paper 16239, May, 1981, pp. 577-589, 
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16242, ELECTRICAL DISPERSION AND SOIL PROPERTIES ate 


WORDS: Additives; Clay minerals; Clay structure; Dielectrics; 
Electrical properties; Electric networks; Expansive soils; Fine-textured soils; “Gg 

Frequencies; Heterogeneity; Hydraulic conductivity; Kaolinite; 
Montmorillonite; Permeability (soils); Soil cement; Soil engineering; 


_ ABSTRACT: In an effort t to develop p a nondestructive method for studying and 
predicting engineering properties of fine grained soils, a simplified electrical model is — 
_ employed to represent the complicated soil system. The parameters of the model are 
both compositional and heterogeneous (structural). For a particular soil system, 
 mametionl values of the model parameters are determined by computer curve fitting _ 
_ theoretical electrical dispersion equations of the model to the soil experimentally 
_ determined radio frequency dispersion curves. The fundamental nature of “radio 
_ frequency electrical dispersion in fine grained soils and the numerical parameters of the 
model offer a technique for studying the engineering behavior of soils in a fundamental 
_ quantitative manner. The following subjects are summarized: the model’s ability to 
-monitor structural changes; the relationship between one model parameter and the 
_ swell potential of soils; and the relationship between another model parameter and the | 


REFERENCE: Smith, Scott S., and Arulanandan, K., “Relationship 
Dispersion to Soil Properties,” Journal of the Geotechnical Engineering Division 
_ ASCE, Vol. 107, No. GTS, Proc. Paper 16242, ney, 1981, pp. 591-604 ah ete 


16243 LIQUEFACTION ANALYSIS AND SHAKING 

KEY WORDS: Case reports; Computer programs; Earthquakes; Layered 
soils; Liquefaction; Sand; onsition tests; Shaking; Shear stress; Soil cyataies 


ABSTRACT: method of is presented for ‘the evaluation of the seismic 


response and liquefaction of horizontally-layered level ground under multidirectional cz 
shaking. A new material model is introduced for representing the behavior of sand 
under cyclic biaxial shear stresses. The responses to multidirectional shaking are 
demonstrated through a series of case studies. These case studies were performed by 
using a new computer program LASS-III which was developed for analysis of multi- 
directional shaking of horizontally-layered ground. The results of a number of analyses — 
are presented in which the influence of the vertical component base acceleration and 
the interaction between the two horizonal components of - pasting are evaluated. 
REFERENCE: Ghaboussi, Jamshid, and Dikmen, S. Umit, “Liquefaction Analysis for 
Multidirectional Shaking,” Journal of the Geotechnical Engineering Division, ASCE, 
Vol. 107, No. GT5, Proc. Paper 16243, May, 1981, pp. 605-627 
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_ 16275 BULKHEADS: HORIZONTAL AND SLOPING ANCHORS — 
_ KEY WORDS: Anchored bulkheads; Anchors; Bulkheads; Cohesive i 9 
Heaving; tests; Sinkage; Slope stabilization; Stability; 


ABSTRACT: flexibility ‘number (DFN), pd, which can be to 
govern similarity between laboratory plates and prototype for anchored bulkheads, is _ 
_ developed. For free earth conditions, the DFN, as it is shown, must be less than five 
_ for loose soils and less than one in dense soils. For designing anchored bulkheads, the | 
DEN, pa, must be analyzed. If pa<5 or pa<1, depending on soil, the bulkhead has 
essentially a free-earth support. In general, properly designed anchored bulkheads in 
z majority of cases conform to free-earth support conditions. Bulkheads with sloping 
anchors in cohesive soils must be designed to the following three conditions involving 

a single stability number, 4c/ H: (1) to avoid overturning, 4c/ H>N)j; (2) to avoid [ 
of the bulkheads caused by sloping abchor, 4c/ H>N,; and (3) to avoid heave 
at the line or bottom of the excavation, 4c/ H>N3. 


be: Anchors,” louhdl of the Geotechnical a Division, ASCE, Vol. 107, No. 

GT5, Proc. Paper 16275, May, 1981, pp. 629-645 
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accordance with the October, 1970 action of the ASCE Board of Direction, which ‘stated 
_ all publications of the Society should list all measurements in both U.S. _ Customary and i> , 
‘SI (International System) units, the following list contains conversion factors to enable readers 

to compute the SI unit values cf measurements. A complete guide to the SI system and its x 
use has been published by the American Society for Testing and Materials. Copies of this — 
= = =68©6-—s publication (ASTM E-380) can be purchased from ASCE at a price of $3.00 each; orders must — 
All authors of Journal papers are being a asked to prepare their papers in this dual-unit format. 
To provide preliminary assistance to authors, the following list ae conversion factors and guides 
a are recommended | by the ASCE Co the ASCE Committee on Metrication. ; 


ra 


inches (i (in. 4 


millimeters 

feet (ft) meters (m) 

yards of meters (m) 
miles(miles) 

= inches (sq in.) square millimeters (mm 

square feet (sq ft) square meters (m 

square yards (sq yd) Square meters (m 0. 836 

acres (acre) hectares (ha) 0.405 

 -- inches (cu in. cubic millimeters 

cubic feet (cu ft) cubic meters (m 0.028 


cubic yards (cu yd) 7 cubic meters (m 0.765 
“pounds (Ib) kilograms (kg) J 
fons (ton) mass kilograms (ke) . 


pound force ~ newtons (N) 
kilogram force (kgf) = mewtons(N) 
pounds per square foot (psf) (Pa) 
z pounds per square inch (psi) kilopascals (kPa) 
U.S. gallons (gal) liters(L) 
acre-feet (acre-ft) cubic meters 
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FAILURE Mopes IN UNDRAINED 


‘Cycuic LoADING 


Ernest T. Selig! M. ASCE and Ching S. Chg? * M. “ASCE 
A typical element of soil in situ to dynamic loading exists under 
4 nonhydrostatic state of stress prior to the occurrence of the dynamic load. _ 
The soil element is considered to be consolidated under these static stresses. 
The principal directions of the superimposed dynamic stresses will generally 
not be the same as those of the static state. Thus, rotation of principal stress 
- directions will occur during dynamic loading. However, the cyclic triaxial test 
is frequently used to obtain information on the behavior of this soil element 
az these conditions. This test has the obvious restriction ‘that cyclic shear 
stresses cannot be applied to the principal planes. Also to simplify the test, 
the confining pressure is usually held ¢ constant and only the axial deviator eae 


_ subcategories, one in which the cyclic deviator stress is greater than the static | 
7 _ deviator stress and one in which it is less. In the latter case, shear stress reversal — 


does not occur on any plane. 
For undrained cyclic tests on oil samples, the strain and pore 
_ pressure behavior of the soil i is dependent upon which of these three categories. 


differences and show “how to the sample failure modes. 
saturated-undrained soil conditions are considered. Experimental data are includ- _ 


ed from other references to illustrate the principals examined. However, the 
ay 
interpretations and conclusions are those of the writers and no mS those of the - 


Prof. of Civ. Engrg., Univ. of Massachusetts, Amherst, Mass. 01003. 
; Asst. Prof. of Civ. Engrg., Univ. of Massachusetts, Amherst, Mass. 01003. 
ote.—Discussion open until October 1, © extend the closing date one ne mont 
~ Note.—Di il October 1, 1981. T d the closing d h, 
a written request must be filed with the Manager of Technical and Professional Publications, _ . 
ASCE. Manuscript was submitted for review for possible publication on June 3, 1980. : 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — 
of the American Society of Civil Engineers, ©ASCE, 107, GS, 
0093-6405 /81/0005-0539/$01.00. 
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i. main categories, one in which the deviator stress is symmetrically varied about _ 
a hydrostatic (isotropic) stress state and one in which it is symmetrically varied 
— 7 about a nonhydrostatic (anisotropic) stress state. The anisotropic case has | 
&§ 
4 


=< 
In this paper, cyclic triaxial tests will be termed as either isotropic or illite 
Isotropic cyclic triaxial tests (ICT) are those in which the sample is consolidated | 
under a hydrostatic stress condition. To the initial hydrostatic effective stress — ed 
state, a cyclic deviator stress is added which alternates between positive and — 
negative values of the same magnitude [Fig. 1(a)]. Thus the maximum cyclic 
shear stress alternates symmetrically between equal positive and negative values 
and has a mean value of zero during each Anisotropic cyclic triaxial 


_ Stresses. The Tatio of maximum, Gres to minimum, G5, effective : stress after 


positive ond negative v value is added to the initial or static deviator stress, producing 


a cyclic shear stress ; which is symmetrical about the initial value (Fig. 1(b)]. 
When the magnitude of the cyclic deviator stress is less than the initial deviator 
stress, no shear stress reversal occurs during the cycles [Fig. l(c)]. In Fig. 
4 1, the axial principal | stress is designated as , to distinguish it from the maximum = 
principal stress, o,, because during part of each cycle in Figs. l(a) and 1(5), 
the axial stress is the minor principal stress. on! 
Cycling stress on a soil sample usually causes a tendency for volume reduction - : 
- of the soil skeleton. When the sample is undrained and saturated with a a 
fluid which is incompressible relative to the soil skeleton, a net pore pressure ; 
increase will occur for each cycle instead of a volume change. The effective - 
stress on the will also” ‘This is by a pore 


2 540 
lion 1S Olten designated by A... AN allermaung Gevialor suress Ol equa 
| 
| 


pressure ratio which is the ratio, ~ expressed as a percent, of the change of 
_ pore pressure, Au, to the initial effective minor principal stress, G,., at _ 
‘The stress-strain behavior of sand samples in the ICT test is ; shown in . Fig. 

2. In a strai -controlled ‘test [Fig. 2(a)], the cyclic strain amplitude is held 
“constant. Because cycling is symmetrical to a hydrostatic stress state, the residual 
pore pressure will be greater (in compression) after each cycle so that the sample 
will become less stiff and the stress required to maintain the strain will decrease — 
each cycle. In a stress-controlled test [Fig. 2(b)], the cyclic deviator stress 4 
= is held constant. Thus the stiffness reduction that accompanies pore 


pressure buildup causes increasing strain n with each cycle. worl 


FAILURE» LINE 


STRAIN 
ike 


STRESS 
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FIG. of Stress-Strain 3.—Stress Paths for Tests 
CurvesfromiCT Tests j|§ (Based on Data from Ref. 


paths for the controlled-strain and controlled-stress triaxial tests 
are shown in Fig. 3 for comparison. The parameter g is equal to (o, — o)/2 
and p is either (o, + o,)/2 or (G, + &,)/2, as appropriate. The dashed ine 
represents the initially applied total stresses. These diminish in the strain-con- 
7 trolled test and are constant for each cycle in the stress-controlled tests. The 
residual pore pressure from cycling is equal to the diffe ference between the initial 


b effective hydrostatic stress, G,., and the effective stress at the end of any 


cycle when the stresses are again hydrostatic, 
¥ In both cases, the pore pressure ratio increases until, after a sufficient nl 


of cycles, ‘it becomes 100%. This condition hi has been termed initial al liquefaction 


= 


| 
| 
AG, 5 — ton 
PB 
| (a) STRAIN — 
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A por pore pressure ratio of 100% means zero s. As Fig. 3 shows, 
when the failure lines extend from the origin (gq = 0, p = 0), zero effective _ 
stress can only occur for a hydrostatic state of stress. In the stress-controlled P 


test, » when the deviator | stress is is applied after the first occurrence of a 100% 


a ‘shearing resistance equal to the applied shear stress so that the stress path 
remains on or within the failure lines. This pore pressure reduction results — 
from a tendency for the soil skeleton to expand on deformation. If the sample 

: is too loose to cause this reduction, then no shearing strength develops and 

no deviator stress can be applied. Instead, the sample will undergo unlimited 


flow deformation and the sample is in a state of liquefaction (2). vate 


a Failure for either a controlled-stress or controlled-strain ICT test can be defined i 
—_ as occurring when the pore pressure ratio first becomes 100% during a cycle. { 
For the controlled-stress test, this leads to the definition ‘‘cyclic shear resistance,’ 

which is the level of cyclic stress q.,. required to roduce liquefaction (pore , 


INITIAL EFFECTIVE CONFINING PI PRESSURE, 


* — COMPRESSION a3 


‘pressure ratio = = 100%) in a of loodieg cycles Q). ‘In an 
test, a equals one-half the cyclic deviator stress. Note that the condition 
_ of 100% pore pressure ratio only exi exists as the stresses pass through the e hydrostatic — 
Because the peak-to-peak strain increases with increasing number of loading 
cycles in a controlled-stress test, failure can also be defined as the level of 
cyclic stress required to produce a given amount of peak-to-peak strain in a 
specified number of loading cycles. This } cyclic s stress is also called cyclic shear 
_ The condition of failure represented by a pore pressure ratio of 100% is 
a unique state in which the sample has zero strength. However, it also” “a q 
no shearing stress to resist and thus strength is no problem. The important | 
consequence is the large strain that can occur upon reapplication of the shear 
stress. The magnitude o of this ; Strain is is a fun function of the dilatancy characteristics < 


— 
FIG. 4.—Isotropic Cyclic Triaxial Test on Loose Sand (3) 

i 


during each cycle between the values ash, and c. ‘These limiting. pore — 
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values correspond to the State, the peak compression state and the 
extension state as show 


Anisotropic Cycuc Tests 


Soil behavior in anisotropic cyclic traixial tests is onesie quite different — 


from that t in isotropic cyclic tests. Some interesting examples are given m 
Schroeder and Schuster (6). ‘T hey used a pulsating load triaxial "apparatus in 

3 (o, > o;) and then loaded by cycling both the axial stress, Ao,, and the a 
stress Ao,. These cyclic stresses were related by Ao, = 0.72 Ao,, based on 

_ the authors’ estimate of the cyclic stress state in the field caused by earthquake _ 
loading. The soil tested was fine Ottawa be 


An example of strain and pore pressure development in an ICT test on loose 
sand is given in ip ne nore pre ratio in ‘mie Pout 

a FIG. 5.—Stress Path in Cyclic Triaxial Test - 
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typical of axial deformation and pc pore pressure 
shown in Fig. 6. In these tests, the mean pore pressure corresponds to the 
stress state during the cycle when Ao, and Ao, are zero. Axial ienpeien 
ain continued to develop at an approximately constant rate while ‘the | mean 
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FIG. 8.—Effects of Void Ratio and ant Stress State on Pore Pressure Development in 

pore | pressure increased during the first portion of the t test, and then became . 

constant. Of particular importance is the fact that the peak-to-peak cyclic strains 

are small and the axial (permanent) strain is large. This is the opposite of the 

end observed for the ICT test. Results for ¢ a range of void ratios, atleast in 


4 
NUMBER OF CYCLES 
Etfect of Void Ratio on Accumulat 
rated Ottawa Sand (6) ‘ 


“Fig. 7 that the magnitude of pore pressure and 2 axial 

in a given number of cycles both increase as the void ratio increases. Belt de _ 

_ The maximum generated mean pore pressure is shown as a function of initial 
void ratio and initial effective stress ratio (K. = G,/G,) in Fig. 8. . With a 
- back-pressure of Uo = 5 psi and an effective confining stress of o,. = 5 psi 

in all of these tests, the pore pressure ein to zero effective me r 
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34) 9.—Cyelic Stress Paths 6 
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- or 100% pore pressure ratio is 10 psi. In no case did the pore pressure reach 


{ 10 ) psi. Perhaps more surprising is the fact ‘that pore pressure decreased in 
q 


: two cases, one a dense sample and one a loose sample, and remained unchanged * 


a - stress path plots in relation to the failure line (Fig. 9). The four different — 
total stress stress oo (TSP) used are compared to the range of n mean effective 7 


& for another dense sample. The reason for this trend can be seen by examining 


| 
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. after 100 Cycles with Stress State (Based on ; 
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MAY 981 
‘stress at ‘the end of the test for all void ratios. le general, ‘Figs. 8 and 9 show | 
that the denser the sand and the larger the initial K. ratio, », the less the generated — 


represents ‘the order of stress states s during cycling ‘that | are increasingly close 
The accumulated axial strain after 100 cycles for a void ratio of 0.68 are % 
correlated in Fig. 10, with the initial consolidation stress ratio qg/p, and cyclic 
stress ratio The strain increases with both sti stress ratios. As q/Po 
approaches zero, which represents the ICT test, the axial strain approaches 
- zero, consistent with the trends observed in the isotropic test. This is consistent _ 
with the observation by Pyke (5) that a 44 stress is needed for permanent a 
Itis clear from this example that failure for anisotropically consolidated samples 
/ 4 cannot be represented by the condition of 100% pore pressure ratio or a cope’ 4 
: _peak-to- -peak strain. The 100% pore pressure ‘ratio may never be reached and 
large cyclic strains may not occur. Instead, large permanent strains will develop. a 
OF AND Amsornome Tests 
axial strain ‘and pressure in | 
“cyclic be subdivided into (sometimes called static) and 
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OF 


FIG. 12. Stress Paths for ICT and ACT Teste nie: 


FIG. 11.—P nd Cyclic Pore Pressure 


eyelic components ‘Fig. “1D. ‘The t corm: strain or pore pressure is 
used to represent the strain or pore pressure at the end of each cycle when 
the cyclic stresses are zero. The term cyclic strain or pore pressure designates — > 

the cyclic variation about the Static value. The peak- ‘to-peak cyclic: strain, ¢.g., 
is the difference in axial strain between the maximum and minimum deviator 


stress conditions. Cyclic stresses may not be symmetrical about the consolidation © = 


_ gtress state. However, for convenience, the static stress state will be taken 
as the midpoint of the cyclic s stress ge path, although this is not necessary for P@. 


the the concepts to be valid. Gy by 
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FIG. 13.—Cyotte ~ Residual Behavior of Pore Pressure and Axial Strain for iCT 
and ACT Tests (3) AS 

sd Results for Oosterschelde sand (3) will be mete to further illustrate the 10; 
_ differences between isotropic and anisotropic cyclic test behavior. Fig. 12 shows ul 

= the total stress paths followed os tests on saturated sand samples having an 


aa for sample bi was : 0. 45 kg /cr cm? . The values of total stress ‘shown by the 
stress paths are those in excess of the backpressure at the end of consolidation. 
_ The axial strain and pore pressure response for the two tests are shown 
. Fig. 13 for 200 cycles of stress. For the ICT test, the pore pressure increased — sed 


| 
| 
— 


- more rapidly than for the ACT test . After 150 cycles, the pore pressure | ratio 
for the ICT test ‘rapidly 100% 6 (Au = 2 kg/c cm while 
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AIG. 14.— Etiective Paths for Shear Stre 
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ratio for the ACT. test approached a value. The permanent 


 » 
(b) STRESS REVERSAL 
g ss Reversal and No Shear — 

_ 
: 


4 increased rapidly near pay of the test. In contrast, permanent: axial 
strain in the ACT test increased at an increasing rate during the test, but the 
" \ Stress path a represents complete reversal of the maximum cyclic shear s stress 
- because the stress path is symmetrical about the hydrostatic axis. However, 
es path ‘‘b’’ represents no shear stress reversal because the stress path ve : 
does ‘not cross the hydrostatic axis. The maximum permanent pore pressure — 
increase in the ACT test is u, in Fig. 12. It is obvious that with no shear = 
a = reversal, as represented by sample b, zero effective stress cannot occur. | 
os stress path with an anisotropic initial state can have shear stress reversal . 
the deviator stress is greater than the static deviator stress. Such a 


FAILED AT LARGE PERMANENT STRAIN 
@=FAILED CYCLIC STRAIN 


~ REVERSAL sid 
ils botsizogal af of 1 | 
7 _ ZONE OF STRESS 
REVERSAL 


> 
‘FIG. 16. —Effect of Stress Reversal. on n Failure Mode 

«eg In the ICT test, t, the c cyclic ‘strain predominates | and in the ACT test with | 

no shear stress reversal, the permanent strain predominates. An examination 

of results for ICT tests and for ACT tests with and without shear stress reversal 
7 was made to deterimine under what stress conditions the transition in strain 

failure modes occurred. The data were obtained from Ref. 3 and represent a 


va a range of initial void ratios. The parameters used to represent the stress state = 
are defined in Fig. 15. ‘The deviator Stress at failure, Fas, ti is a function of 


| 
bom sulist 
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s The failure modes are indicated in Fig. 16 for a variety of si static and cyclic sd 
_ stress conditions. The solid circles indicate the tests in which failure occurred 
by cyclic strain, and the open circles indicate the tests in which failure occurred ‘ 
bag permanent strain. Lie transition was precisely defined by the | condition of = 
permanent axial strain ¥ was as greater than the cyclic axial strain; when shear stvens ba 
reversal occurred (04, = o4,); the cyclic strain was greater than the permanent ee 
Ina cyclic triaxial test, failure can be defined by a specified magnitude of be ; 
‘Strain developed in the sample. Two distinctly different failure modes —_ 7 
one in which the permanent axial strain reaches the specified value first and 
one in which the cyclic strain reaches the specified value first. Test — 
ee that shear stress reversal was the governing factor distinguishing these — a 
two modes. When n no shear stress reversal occurred, then the permanent strain 
mete 
_ The largest pore pressure increase involves the condition of shear ane 
“reversal. Only for the case of shear stress reversal can 100% pore pressure 
ratio and thus zero effective stress occur. However, the condition of shear 
- stress reversal or no reversal is directly associated with the condition of whether 
or not the cyclic deviator stress is greater than the static deviator stress. When : 
‘the cyclic deviator stress is greater than the static deviator stress, cyclic strain ; 


failure and shear stress reversal occurs. 


Although the same magnitude of strain may be used to define failure, whether 


the mode of failure is cyclic or permanent strain, these two modes: are ‘quite 
different and should not treated interchangeable. 


The research on which this paper is based was partially supported by funds 


_ from the National Science Foundation Earthquake Hazards Mitigation program. 
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purpose of thi this paper is | paper is to define and the of steady -state 
in particulate media, with special reference to soils. 
The steady s state of deformation (8) for any mass of particles is that state 
in which the mass is continuously deforming at constant volume, constant normal 
_ effective stress, constant shear stress, and constant velocity. The steady state 
of deformation is achieved only particle orientation wad a 


> _ The steady state of deformation occurs during undrained tests on fully saturated, J 
loose sands after liquefaction has been achieved. During ‘the steady- -state 
a - def ormation, the original structure of the specimen has been completely destroyed ; 
and reworked into a new “‘flow structure,” which was originally proposed by 
~ Casagrande (4); i.e., due to the process of deformation, the particles become 
oriented such that the shear stress needed to continue deformation eventually — 
reaches a constant value. A test to demonstrate flow structure in sands will 
Steady-state deformation also occurs in n drained tests on sands at large arge strains. 
The velocity of deformation in the usual drained test is slow, on the = 


of deformation 1 remain constant pone | 


of 0.1% /min-1%/min. Nevertheless, when large strains are reached, a “flow ; 
structure” develops. If the strains are stopped, the specimen is no longer in 

the ‘steady state of deformation. ‘However, the steady = can be achieved 

again if large enough shear strains are again applied. 

7 4 Steady-state deformation occurs in clays in drained shear v when the “‘residual’” : 
shear strength is reached at very large strains (1,6,14). At this stage the particles 
in the shear zone are oriented as well as they can be under the applied stress 


_ State, such that the specimen exhibits a minimum resistance to continuous shear 
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{The steady state of deformation is also ‘aii in undrained tests in soft 
clays at very large strains. The remolded strength of undisturbed clays as 
determined by the vane shear test is nearly i a steady-state shear strength. In 
- the shear zone a , statistically constant structure is achieved if the vane is rotated _ 
at a constant rate. The ratio of the undistrubed peak strength to the steady-state - 
_ Strength is its sensitivity (13). Some drainage takes place away from the shear 
zone during vane tests on soft clays. Therefore, the strength measured in that 
test may be somewhat higher than it would be without any drainage. Thus — 
_ the steady state of deformation is not truly reached in such a test unless the 

_ The steady state of deformation can occur in any particulate mass and for. 
any loading and drainage conditions that can break down the original structure. 6h 
and put it into a new ‘“‘flow”’ structure. The process of structural breakdown 
of a particulate mass under the influence of shear stresses should be viewed _ 

a as a continuous process that starts under the first small increments of shear — 
3 stress and is complete only when the steady state of deformation is achieved — 
and all memory of the original structure is lost. 

_ The initial structure of the specimen, and the corresponding stress-strain oENe 
are dependent upon the way ‘in which the specimen was prepared. It may be- 4 
an “‘undisturbed’’ specimen, which was prepared in nature, or it may be a 
 laboratory-compacted or sedimented specimen. In contrast, the flow structure, _ 

and the corresponding shear strength, for a given soil or other particulate mass 

j s not dependent on initial structure but are postulated to be dependent only 

7 4 on: (1) The effective normal stress during steady-state flow; and (2) the velocity — 
of deformation. Thus, for a given soil and a given velocity of deformation, — 

7 there is a unique relationship between void ratio, effective normal stress, and 

Shear stress during steady-state deformation. a 

_ It is vital to understand that the steady state of deformation is not a static 

condition. The steady state exists only during shear deformation. If shear 

“deformation is stopped, then the particles move into a different structure than 

_ that which exists during steady-state deformation. If shear is started once again 


: under the same test conditions, the specimen does not again reach the steady © 


7 ‘state until the strains are large | enough to create the flow structure again. — i» 

al The unique relationship between void ratio and effective minor principal stress " 
7 during steady-state deformation is shown for one sand by the dashed steady- iste 
dine in Fig. 1(a). Consider a specimen that is compacted to a void ratio and 
effective normal stress which plot on the steady-state line, point A in Fig. - 
l(a). That specimen is not in the steady state of deformation. The specimen, 

as prepared, has a certain structure which n must be destroyed by the shear 

‘ stresses to create the new flow structure that exists only during constant velocity 

_ deformation. Thus, the steady-state line is a locus of points at which the soil 


can deform continuously at the steady state. It is not a line | line connecting static 


; For clays, which are composed chiefly of platey-shaped grains, it is straightfor- 


ward to visualize flow structure, i.e., that grain orientation which occurs an 
large ‘Strains in drained tests when the residual is reached, which 
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- flow structure ‘in sands because the bulky shape of the grains of a typical 
sand do not seem to lend themselves to any significant degree of orientation. _ “3 
_ Although the enormous effects of structure, i.e., the arrangement of particles, 7 
on the stress-strain curves of sands is difficult to visualize, it has nevertheless 
been shown by several authors €. g., 4, 11) that such is the case. It is ce ertainly 


clear that if one is shearing a fully-saturated sand in undrained shear, then = 
rather minor changes in the arrangement of the could lead to 


changes in pore pressure. 


_ For the purpose of observing th the consequences of flow structure in sands, _ 
undrained test was performed on a fine sand called Banding 
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deformation 


4 


. 
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a very effective consolidation pressure (3,. = 
conditions were chosen so that the point representing the state of the specimen a 


0.2 kg/cm). These te 


after consolidation, point Cc in Fig. \(a), would lie below the steady-s “state 


4 


chosen t to > observe whether the specimen “would ‘dilate to the steady-s -state line, 


x or whether the steady-state line would be crossed. [Note: The selections of 

void ratio and effective stress were based on the known steady-state line as 

* measured by Castro (4). The steady-state line shown dashed in Fig. 1(a) was 
drawn parallel to Castro’s Steady-state line but passing through the point S, baa 


which represents the steady state for this specimen 


The pertinent test conditions 


— 
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g 2. Structure—Compacted at 5% water content in 10 layers using a load of 
884 gm on a 0.5-in. (12.7- mm) diam rod with 18  tamps per li layer. The procedure _ 
Initial State—@,, = 2.8 psi (19.6 kPa), y,. = 99.2 pef (1,590kg/m*) (Percent 
of maximum density = 89.1%), 
4, Test Conditions—Axial compression. Strain control at 0. 42% min. Skemp- 
ton’s pore pressure parameter, B. = 0.95. Back pressure, u. = 114 psi iad 
7 KPa). Specimen diameter = 1.41 in. (35.9 mm), height - = 3.46 in. (88.0 mm). — a 
___ The results of the triaxial compression test are shown in Fig. 1. The specimen _ 
: first slightly contractive, in that the pore pressure increased under ia q 
first small increments of stress (5). At large shear stresses the specimen became 


dilative, the induced pore pressure dropped, and the effective minor principal 
stress increased. Note i in n Fig. Me) that the early part of ‘the s stress. strain curve 


minor principal stress started i increasing. 4 
_ During the portion of the test when this specimen was dilative, 6 an > 
on the state diagram crossed the steady-state line at point S’ in Fig. l(a). At 
_ point S’ the void ratio, the shear stress, and the effective normal stress puna 
on the steady-state line in Fig. (a). Yet, steady-state deformation did not occur 


‘The explanation for this behavior is that the structure of the specimen at 


The structure had been dilative and it was still internally-braced (A. Casagrande 
has frequently referred to “‘internal-bracing’’ within a dilative sand during his 


=e S’ was not that flow structure which would permit steady-state deformation. 
lectures at Harvard University), so to > speak, at this stage. As the shear stress 


- it reached point Min 1 Fig. 1(d). Then the effective stress was 1s high enough — 4 
so that the skeleton started trying to decrease in volume. The specimen sf 
become contractive. Therefore, the pore pressure increased and the effective 

> decreased correspondingly. As the strains increased and the effective 

was dropping, the particles had enough freedom during | ‘Shear so that 


they could become re- -oriented under the influence of shear strains, ;, thus causing 
the process of pore pressure buildup and grain orientation to continue. Finally, 4 
when the strains were large enough, a steady state was reached at point S, 
after which there was no further change in shear stress, effective normal stress 
or pore pressure. The velocity, i.e., the strain rate, was essentially constant 
throughout this test. Thus the steady state of deformation had been reached © 
at a strain of about 14% and this steady state persisted between 14% and 19% _ 
strain, at which point the test was terminated. Again, at points S and S’, the a 
void ratio, effective normal stress, and shear stress were identical. Yet the 
steady-state deformation occurred at S but not at S’ Son solaediina ina ine 1 bioy 
The only variable that might be different at points S’ and S in Fig. | is 
the structure of the soil. Thus t the | difference | in structure at these two points “| 


must be the cause of the major difference in behavior that is evident at the 
corresponding points on the stress-strain curve. 


Soil Tested—Narrowl Osubangular 
4 


‘alle the state of deformation is not until ‘particle 
breakage, if any, is complete. Marsal (7) has indicated that particle breakage _ 
probably causes changes in the peak shear strength. Also, the flow structure | 
< the corresponding steady-state strength can be expected to be affected 
by particle breakage. In fact, it may be nearly impossible to reach the steady _ 
3 _ State in the laboratory if particle breakage continues, e.g., in a narrowly- — 
angular sand at high effective stresses. This ara t of steady- -state deformation 


Sreapy State ann Camcat ‘Stare—Bruer History 
The kernel of ideas is the concept that ot soll other granular 
materials, if continuously distorted until they flow as a frictional fluid, ne 
will come into a well-defined critical state... 


Critical state D 


ment:in 


: FIG. 2.—Relationship between Critical State and Residual nest of Clay (after 
See Ref. 12, p. thyast 100 Uiw | 


.* oe tear, rub, scratch, chip and even bounce against each other 


y we keep our distance and see flow. At close Tange 
be 4 would expect to find many complicated causes of power dissipation and 
some damage to particles; however, we stand back from the small details 
and loosely describe the whole process of power dissipation as ‘‘friction,”’ 
‘neglecting the possibilities of degradation or of orientation of particles. aa 
q The foregoing quotation indicates that the authors view critical state as an ideal ‘Se 
state in which the particulate mass deforms continuously but has no preferred 
structure. Also, particle breakage is neglected. 
This interpretation is supported in Fig. 2, which is a reproduction of Fig. 
- 18 from Schofield and Wroth (12). The critical state as defined by Schofield - 
6 and Wroth is shown to lie between the peak and residual, i.e., at @ point at 
Which steady-state deformation i is not occurring. wis ib 
The writer submits no mass continuous us deformation — 


a 
4 
“during the process of continuous derormation. If the mouon were Viewe 


at a constant state of effective stress, and at constant 
ee without having a preferred orientation of grains, i. e., a flow structure. 
addition, the shear resistance at large s' strains, during steady-state deformation, 
is affected by the of this steady-state 


end, or terminal, state of the shearing ‘process. fo 
The steady state of deformation as defined herein is ‘the as 
envisioned by Casagrande (2), when he proposed the existence of a critical — 
void ratio, or critical density. Quoting from Casagrande (Ref. 2, p. 262), who ‘F 
was describing drained direct shear tests ¢ on sands (material in brackets added 
During this drop {from the peak strength) i in shearing stress, the sand % 
continues to expand .. . finally reaching a critical density at which a 
continuous deformation is possible at the constant shearing stress S, . . . 
Therefore the curves representing the volume Changes during . 


density when the stationary condition is established. — 
At the time of this writing Casagrande had not considered of 
_ “flow structure,” which he later introduced. Also, he made ‘no mention of 
th the velocity of deformation and did not apply th the idea to clays. 
‘The steady state is similar to the concept “presented by Roscoe, Schofield, hg 
and Wroth (Ref. 10, p. 28) in the following description of critical state: vue 
. In a drained test the critical voids ratio state can be defined as s the 
‘= hat stress of a sample at which any arbitrary further increment of — 
. shear distortion will not result in any change of voids ratio. In any series — 
of drained tests the set of critical voids ratio points so defined can be 
Bhs <2y to lie in or near a line on the drained yield surface, 
_ In an undrained test the sample remains at a constant voids ratio, but 
— effective stress p will alter to bring the sample imo an ultimate z. 
— that the particular voids ratio, at which it is compelled to remain 
+ during shear, becomes a critical voids ratio. . In any series of undrained | ry 
tests the set of critical voids ratio points so defined can be expected — 
- to lie in or near a line (not necessarily the same as above) on the undrained = 
the results of drained and tests show that there is, 
fact, one unique line to which all loading paths in (p,e, q) space converge, 
a a then this will be called the critical voids ratio line (CVR line, for short). aa <3 
ta The drained and undrained yield surfaces will then coincide along the 
CVR line. These two surfaces may be identical and, if a, the common — 


= difference between the steady-state line at the critical voids ratio (CVR) a 


i 
: vital to recognize this behavior (8,9). In many types of laboratory equipment , 
of 2 currently, used, one cannot reach the steady state. This fact does not diminish 
fw 
=¢ 
a 
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"STEADY STA STATE 
does not exist unless | deformation occurs at constant velocity. In contrast, the 


—- of the description of critical state in the aforementioned drained 
_ case is that, if shear is stopped, after the critical state is reached, the specimen — 
remains in the critical state. These authors did extend Casagrande’s original _ 
concept to clays. However, it is the writer’s interpretation that in applications — 
_ they did not think in terms of a steady-state structure of oriented — 
_ but rather considered the structure on the CVR line to be random. 
_ Rowe (Ref. 11, p. 520) describes constant volume shear as follows: © 


= The ultimate state of continual remoulding at large strains for an actual Biv 


- while the other half have failed and are making contact. An additional : 
7. energy loss results from friction on the inclined planes as compared with > 
the flat plane, even though the mass as a whole moves as if it were a 
sliding on such a plane at n No volume change. rin 
This description of constant volume shear is essentially the same as the description 
_ given herein for steady-state shear. However, in this description the velocity ; 
of deformation is not considered. Nor is there any explicit statement that a 
_ statistically oriented structure, as » as opposed to a random structure, might a 
_ Im general, the concept of steady-state deformation has been embodied in 
the previously quoted writings. However, it has not generally been recognized | 


¥ 


j that continuous deformation is a necessary condition of the steady state. Nor 


has it been recognized clearly that a flow structure develops at the steady 
- state in all soils—clays or sands—and for both drained and undrained conditions. ; 
* These generalizations of the original concept are ‘important to the proper qj 


understanding of the ‘stress-deformation properties of particulate media. 

Pe ‘The term steady state of deformation was selected because: (1) The term Sow 
_ “steady state’’ is used i in analogy to steady-state (nonaccelerative) flow in liquids; 
and (2) the term “deformation” is included to emphasize that the steady state 
does not exist unless deformation is ongoing. If the velocity of deformation 

7 zero or changing, the specimen is not in the steady s state _ deformation. bic’ 


term critical state is not used because: om, 

1. The ‘steady ‘on is reached only after complete orientation and breakage Py y 

: of grains, or ‘‘flow structure,’’ has been reached, which is not the case in his 
the current (12) use ofthe termcritical state, 
_ 2. The term ‘“‘critical’’ usually is used to refer to incipient breakdown of : = 
‘some sort. In contrast, at ‘the steady state the original st structure has been completely — 


deformation is really a noncritical state. hh bra Bow 


of “state | deformation is postulated to be to 


wen it has reached the steady ‘state of | deformation. The steady state ‘of | 


j 
1 

| 


all, rock, corn in bins, pulp in | 


chemical wastes, coal piles, etc. jo 
_ Several practical applications of the dae of steady-state deformation are 
described by the writer (8). One of the principal advantages of the —— 
is that it permits rational organization of stress-strain data. _For example, 
determination of the effects of various parameters, such as soil composition, 
_ void ratio, sample preparation technique, and effective stress, on stress-strain | 
curves can be more easily understood when viewed in terms of the concept _ 
- of steady state. Also, if a particular test result does not fit the pattern shown» 
in other tests, one can question its validity and start investigating the cause a 
_A few additional practical applications are 
1. The meaning of pore pressures measured in-situ before, during and well 
after ; a failure, and the use of those measurements to analyze the failure can 
2. The selection or shear strength when analyzing the stability of and 


“3. No one apparatus | can be used to” measure the stress- curves for 
all soils and test conditions. The selection of apparatus for particular applications _ 


deformation helps one to understand the shape of stress-strain curves for various — 
test conditions, this concept is quite helpful in the selection process. 
4. The relationship between the stress-strain curves of sands and those f 


- clays can be clarified. The terms ‘ ‘normally- consolidated”’ and “‘loose”’ as ; opposed 
“‘over-consolidated”’ “‘dense”’ can be in terms of soil behavior. 
‘The steady state of deformation has been described and explained. The concept _ 
of steady-state deformation originated with Casagrande’s (2) idea of “critical 
void ratio.’’ That idea was extended later by Roscoe, Schofield, and Wroth 
(10) to include drained and undrained shear of clays as well as sands. In this 
aper the following items areemphasized: 
aig ‘Steady-state deformation requires a constant velocity of deformation and % 
the steady state exists only so long as deformation continues. | sa: 
* 2. The structure of a specimen during steady-state deformation is very special 
and quite different from both its initial structure and the structure of the specimen 
before steady-state deformation is reached. 4 
3. The steady state can be achieved only after grain orientation (in sands — 
as well as clays) and grain breakage, if any, are both complete. These conditions _ . 
‘normally can be attained at large — those that 


= 
is related to the shane of the stress-strain curve Since the concent of steady-state 
; 


“aie residual strength is reached in drained tests” 
on clays or sands are examples of steady-state deformation. Also, continuous 
_undrained shear of saturated soft clays, as in a vane shear test in-situ, and : 
the continous flow that occurs after liquefaction of a loose saturated sand . p 
undrained shear are examples of -state 
wari 
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end of consolidation of test; 


specific gravity of soil solids; 
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- maximum dry unit weight as measured in dry state on vibrato 
= axial strain (major principal strain); ye 
_&, = Major principal effective stress; and 
ic = minor principal effective stress. 
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The evaluation of the liq liquefaction at nuclear power station 

_ incorporated two distinct procedures of field soil sampling and laboratory 
preparation prior to cyclic strength testing. The initial procedure consisted of 4 
Pitcher barrel sampling and laboratory trimming of samples by extrusion into” a 

x a brass tube mold. Low cyclic strengths were obtained for those samples obtained | 
by such a procedure in comparison to strength estimates based on Standard 

. Penetration Test N values obtained for the site ‘soils. A second program was 7 

: undertaken to ascertain if more representative cyclic strengths could be obtained 
__ by a different procedure of field sampling and laboratory preparation. In this _ 

* program, block samples were hand-carved from two 6-ft (2-m) diam, dewatered, — 
auger borings. Samples were carved in into 9- 1/ 2-1 -in. - (24-cm) diam by 12- in. (30-cm) 


_ The subsequent laboratory preparation techniques included either further a 
_ carving, freezing and carving, or saturation and carving into test-size specimen. ie 
_ The following paragraphs document the two sampling and laboratory prepara- e, 
tion schemes, and compare the results of the cyclic strength tests for ‘specimens om 7 


= upper 150 ft ft (15 m) of the site soils consists 5 primarily ¢ of po to very 
dense sands, silty sands, and silty clays. Scattered gravels were also aah 
The materials are alluvial and occur in discontinuous layers, lenses, and pockets. 
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soils and their general properties is shown in Fig. Bs, 
‘Three areas at the site, all having basically the same soil conditions, were 
_ investigated. These were designated as areas Ul, U2, and U3. The central — 
portion of the site had been cultivated for many years. The associated irrigation 
created a perched ground-water condition within the upper sands and the 
underlying clay formed the aquitard supporting the perched water. Concentration 
_ of irrigation in the central portion of the site formed a perched water mound — 
- with shallower depths of water near ar the center, and greater depths beyond. 

As a result, water levels at the three | areas differed. Water levels for _ 


a4 


AVERAGE 
UNIT WEIGHT AVERAGE VOID RATIO | AVERAGE SPT N-VALUE 
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1.—Soil Profile at Site 


_ areas U2 and U3 are shown in are. | Area Ul encountered water below the 


Ss The field sampling program first consisted of Pitcher barrel sampling in 

- _ rotary-wash borings followed at a later date by block sampling i in large-diameter 
auger borings. The Pitcher barrel sampler utilized consisted of a thin-walled, 

seamless, steel tube within a ‘thick- walled cutter barrel. In operation: (1) The 

_ samples enters the hole [see Fig. 2(a)] while the sliding valve directs drilling 

_ fluid through the thin wall shelby tube to the hole bottom for thorough elses 

4 action; (2) the instant the shelby tube touches bottom [see Fig. 2(b)] it telescopes 

_ into a cutter barrel to close a sliding valve which automatically Givers fluid 


&g 
pili? 
a sample-loaded shelby tube is removed from the hole [see Fig. %e)]. The : 


sampler i is used i in dense to very dense materials. The; of the s sampler 
is that as the inner sampling tube enters the soil, the materials outside of the 
_ tube are removed by the cutter barrel thereby reducing confinement and outside - 
friction, and permitting penetration of tl the i he inner tbe without premature blocking | : 
Inner tubes were 3 in. (7.6 cm) in outside diameter, 2.875 in. (7. 30cm) in 
7 inside diameter, ¢ and 36 in. (91 cm) long. Tube tips w were designed to . provide 
“3 an inside clearance ratio of 1%-2%, and an area ratio of 10%-15%. Tube somplee 
were obtained in the sands at each of the three site areas. st” 
& _ A second field sampling program was initiated at site areas U2 and U3 to = 
_ obtain undisturbed block samples in large-diameter borings. The boring at area _ 
-: U2 was designated U2-LB-1, and the two borings at area U3 were — 


U3-LB-1 and U3-LB-2. The three were located within the respective 
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SAMPLE PROTECTED FROM 


FIG. 2—Field Operation of Pitcher Barrel Sampler: (a) Inserting Sampler in "Drill 


In general, block samples were drilling [7- ft diam], 
with a caisson drill rig using a short-flight auger with extender bars. Fig. 3 
shows the sequence of sampling and drilling by which the boring was advanced. 4 
Borings | were drilled to within approx 5 ft (1.5 m) of the water table. To protect 
personnel working in the excavation, 6.5-ft (2.0-m) diam corrogated metal pipe 7 

7 casing was lewered to the bottom of the excavation [see Fig. 3(b)]. After : 
_ being anchored at the surface, a 6.0-ft (1.8-m) boring was advanced in stages _ : 
as the sampling operations progressed [see Fig. 3(c)]. Sampling was initiated 
by drilling a 24-in. (61-cm) diam “‘pilot hole”’ in the center of the excavation. 
The pilot hole provided a place for the sampling personnel to stand while samples. - 


were carved from the floor and wall of the excavation. Well points placed 


= VALVE SEAL 4 4 aN 
3 Side the auger boring provided dewatering. 


Each block by hand into a a 9- -1/ 2-in. (24cm) di diam by 
7 a 12-in. - (30-cm) high lucite mold. The use of a cylindrical instead of cubical 2 
- mold was intended to prevent the development of stress concentrations at the ea 
- comers a and thus, minimize sample distress. By using a clear sampling container, © 
: _ the degree of disturbance around the perimeter would be visible, enabling the 
field engineer to assess the quality of the sample obtained. 
_ Each lucite cylinder was 0.5 in. (1.3 cm) thick and machined to produce 7: 
a sharp cutting edge and a beveled top. After machining, the diameter and 
- out-of-round dimensions were measured, and the volume of the cylinder deter- 
mined. Each cylinder had accompanying base plates, Tubber bolts, and gaskets _ 
7~ to contain the sample for transporting. These were weighed together as a unit 
: - establish the sampler container weight. The gaskets were intended to seal 
| 
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BORING CASED AND ADVANCED s 


3.—Typical of Block-Sampli 9 Operati 


the: soil within — at a constant water content to ‘insure ‘that soll 


. & 


Figs. 4 and 5 present details of the block-sampling operation. By careful — 
- trimming the soil away from the leading edge of the samples while exerting & 
a slight downward thrust on the cylinder, a soil sample was eventually enclosed 
in the lucite cylinder. During sampling, disturbance could be qualitatively 
evaluated, and noted by examining the stratification within the clear lucite 


After enclosing a sample, the base and top were trimmed flush with - 
lucite mold and sealed with the base plates which were then bolted tr ol 
The bolted sample container was also sealed with wax around both cylinder — 


to base plate contacts. Fig. 6 shows the the completed —* ready for transport. 
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CYCLICSTRENGTHS 
\BORATORY SAMPLEPREPARATION, 
a _ Both the Pitcher tube and block samples were carefully packed in foam-lined oe 
boxes and transported to the laboratory by truck. Laboratory processing of out 
samples from the tube or lucite cylinder to the cyclic triaxial equipment involved _ 
separate procedures. Pitcher tubes were cut into 8-in. (20-cm) segments, using 


—tLucite Chamber Prior to Sampling Operation 


Samples selected for testing were then trimmed from the original 2.875-in. ; 
 (7.30-cm) diam to 2.50 in. (6.35 cm) by mechanically extruding the samples ~ 
from the tube through a 2.50-in. (7.30-cm) cutting head, into a 6.0-in. (15-cm) 


~ long brass mold. The cutting head and mold were prelubricated to minimize a 


sample disturbance. The cutting head had an inner clearance ratio of zero (flush — 


walls), and an Tatio of ss 


a tube cutters for the steel and a knife for the soils. Metal burrs at the tube _ 
» ends obtained from the cutting process were reamed off. Exposed soils within 
consisted Of initially Nooding the block mangle with 

&g / 
; in 5%. Both ends of the sample were trimmed P- 


flu ush with the ends of the brass mo oid, after which the s 
weighed, and set up in a cyclic triaxial pressure cell. 
_ Undisturbed samples [2.5 in. (7.3 cm) in diameter x 6.0 ) in. (15 em) i in length] 
were trimmed from undisturbed block samples [9-1/2 in. (24 cm) in diam 7 
12 in. (30 cm) in | height) using two trimming methods: (1) Cut and shape; ; and = 
(2) freezing. The cut and shape trimming method consisted of trimming the _ 
full-size block sample, according to the method described by Lambe (2). — 


separate techniques were utilized for this method to maintain the sample integrity 


FIG. 5.— —Undisturbed Block Sample a after Sampling 
during the _ The first technique involved g the existing 
4 sample without any additional support. This was done if the sample exhibited 
enough cohesion or cementation to enable it to stand by itself. In contrast, © 
. b the sample was unable to stand by itself a second technique utilizing capillary 
tension was employed. This technique consisted of flooding the block sample — 
3 water, allowing the sample to drain, and then trimming. The ‘Temeining 
“water present in the soil matrix tended to hold the material together by capillary a 
of the and drained block sample was 


> 

| —— 


6.0 in. (17.5 cm) in length). The lucite cylinder was directly on 


undisturbed block material. Material was then removed from around the base % 


of the cylinder as ‘the cylinder v was lowered to add support to the soil column — 4 

i - being produced. The smaller test-size sample was subsequently removed from vt 

: i the tube upon completion of the trimming operation. The cut and shape method | “a 
‘ of sample trimming was used on fine soil material. 

_ The second trimming method employed on block samples was the ili 7 

= method. This method consisted of initially flooding the t block sample v with = 


and then | allowing the sample drain. After | drainage had taken the 


Block in Chamber Ready for 


liquid nitrogen as opposed to allele refrigeration \ was ; the elimination 
of the formation of ice lenses in the sample. The freezing process did cause : 
a small amount of observable volumetric expansion in the sample which probably - 
would have gone unnoticed i in a smaller-sized sample. | 
Inspection of the block sample ; after freezing indicated that the expansion — 

- appeared to be uniform. The uniformity of the expansion in the sample did 
‘hot appear to have changed the individual particle orientation or soil fabric, = 
and therefore, it minimized disturbance. After freezing, the sample was trimmed — 

ws the method described by Lambe (2), and placed in a pressure ¢ coll to thaw 

Cycuc Triaxia. Test Resucts 
_ Both tube samples and block samples were tested in accordance with maton 
cyclic triaxial as examined recent literature (4,6) -As summa 


until frozen. e advantage of the fast-freezing process by the use of 
 .... 


- ; of the static properties and cyclic strengths from both the tube and block samples _ 

obtained from depths between 30 ft (9 m) and 50 ft (15 m) for the three site — 

TABLE 1.—Summary of Cyclic Triaxial Test Results 


40-42.5 


3638.5 
—30.7-31.8 
-43.4-35.0 
34.9-35.8 
-35.9-36.7 
4343.8 
43-43.8 
 35.9-36.7 
 43.7-44.5 
44.7-45.6 


_ “P = Pitcher tube sample; B, = block samples hand-trimmed from the boring (9-1/2 ft 

and hand-trimmed into test size samples (2.5 in. x 6 in.); and B, = block samples hand-trimmed 

from the bo and frozen before hand-trimming into test size samples. _ a ak 
ring g into test size samp 


_ above and below the perched water table were tested. As will be seen, strengths — 


area } ounds| De, | asa pounds} Cycles | to 
ing | ple Depth ae uare| €,,= | 03, = ple 
Test| ing bic | milli- | cent-| o,,/ |sq 
num: ters| age | | inch | 5.0% | 1.0 = 
infeet | fest“ | | | (9) | (10) | | 
21 0.25 | 28.0 | 
3-B3 | 9 | 325-35 | 101.1 | 0.17 | 41.0 | 380 | 350) Po 
111.1 | 0.60 | 7 0.20 . 
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-B9 | 18 | 37.39-5 103.3 | 0. 
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071 | 032 | 12 | | 260 
26 |U3-LB-1| 10B 8.0} 0.19 | 17 | 0.40 | 28.0 
27 |U 0.25 | 30 
97.91 0.29} 8 | 0.25 bog, 
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Pa. were generally the same for these samples. B Block sa sins mile ne obtained 
_ Summary plots of laboratory cyclic strengths in terms of the cyclic stress i 
a ratio, 04, /203., versus the logarithm of the number of cycles of load application 
_ for two failure criteria, are presented in Figs. 7 and 8. Figs. 7(a) and 1b) 
Present data for the tube samples on failu re criteria of 5.0% double amplitude 

NUMBER OF CYCLES 


a 


A 


| Samples, 30° -60', Triaxial Strengths: (a) Sa = 5.0%; 
a3 comparison, Figs. 8(a) and 8(b) present similar information for block samples. 
| _ The position of the average curves visually drawn through the data was based | + 
os on typical shapes of strength curves developed from shaking table tests by — 
a te De Alba (1). A review of Figs. 7 and 8 indicates that similer cyclic strength | 


are produced, whether based o1 on the = = 5. 0%, or initial 


ii 
in 
| 
WOTE: WUMBER ADJACENT TO DATA POINT CORRESPONDS TO TEST NUMBER IN TABLE 1; | FT=0.305 


criteria for either the tu whee samples or or samples. 
_ The average soil characteristics (dry density, D,,, and percent of fines) for 4 ; 
both the tube samples and block samples are summarized in Table 2. Table 
_ 2 shows that samples obtained by the two sampling techniques compare very 
4 well in terms of Ds and the percent passing the No. 200 Standard ee 
Society for Testing and Materials sieve. No — in material type are — 
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‘small layers and lenses of sand and silty sand and clay within individual test 
specimen. Although the average properties of the entire specimen generally Ji 
were similar, behavior during ¢ cyclic have | been affected by peculiari- 


er 


0.8 

| 
= 
al Strengths: (a) = 50%: (6) AU/oy > 
ted. However. phvsicallv, e 


STRENGHTS 


ties within individual specimen. For t this | reason, and because ov the waa 
of more than 5% fines, relative density determinations were not made, but J 
rather dry densities were used in subsequent evaluations of cyclic ingereend 7 


‘The comparison of cyclic strengths for similar samples obtained and = 


A 


by two mapenee procedures permitted observation of Several wegte and indica- 


_ Pitcher Tube Samples.—The cyclic test results shown on Figs. 7(a) and 10) 
are consistent, and the scatter is low. The overall level of cyclic strength, a 
_ indicated by the average line visually plotted through the data, is low when ~ 
- considering the soils generally had N values ra ranging from 30-65 in the — 
50 ft (15 m). Disturbance, although suspect, was still difficult to prove or quantify. 
In addition to the low strength, disturbance could be inferred by the minimal 
che edges. THE TABLE 2.— 2.—Average Soil ‘Characteristics 
preparation per cubic asa percentage 


Sample type oom... Average Average 
98.1-1 ‘| | O41 | 738 | 0.41 | 7- 38 


have controlled response and produced the consistent, low scatter 
4 Block Samples. —By testing the block samples, it was possible to observe © 7 


on Figs. 8(a) and 8(b) show that, overall, the strengths were high, and the 
variation in peracsind and possibly percent fines (particularly the plastic fines) 


65 gm /cm contained somewhat high fine- grained 
‘material. For the lower strength curve, all but one test specimen (Test 34) — 


F 
“4 
— 
— 
7 
| 
| 
Average 
° 
— d exist among 
- the samples tested, but were not indicated by the cyclic test results. Confusing — 
. this observation, however, was the fact that individual test specimen contained __ 
a small layers, lenses, and pockets of materials that differed from the average _ 
and quantify the amount Of disturbance that resuited [rom the Pitcher sz J 
a followino la ratory nrenarati ocedu e test results s 
| 


Test 3 34 was attributed to the low fines content (6%). 
hig: success of the laboratory preparation procedure was somewhat verified 


= from block samples prepared by both the trimming method and the 
Pitcher Tube and Block Samples.—Cyclic strengths based on = 5.0% from 
“Figs. 7(6) and for Pitcher tube samples and block are compared 
_ in Fig. 9. To quantify the minimum increase in strength | obtained using the 
_ = sample procedures, a lower-bound cyclic ‘strength curve is plotted in Fig. — 
8. The corresponding average dry density of the block samples lower-bound 
_ . cyclic strength curve is about 99 Ib/cu ft (1.59 gm/cm’), as compared to 103.6 
| J —Ib/cu ft (1.66 gm/cm’* ) for the tube samples. _ Fig. 8(b) shows that the average | 
sample cyclic strength 5%-24% lower cycles-30 cycles than th the 


4 
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BLOCK SAMPLES 


AVERAGE OF 


of Thin Wall and Block “cyclic ‘Taal 
: corresponding block sample ‘cyclic. Seni A direct comparison of the strength 
_ curves shown on Fig. 8(b) with the curve for tube samples (Fig. 9) indicates 
- block sample strengths 24%-110% higher (at 30 cycles) than for tube ‘samples. 7 
Cyclic triaxial liquefaction test results have been compared for soil samples 
_ obtained by two different sampling and laboratory preparation techniques. Results 7 
indicate a reduced cyclic strength for samples that were obtained by xtison ‘i 
followed by tube cutting and trimming byt mechanical extrusion 


= 


= 
J 
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| 
hand-carved, block samples which were further hand-trimmed into test-siz 


be caused by the destruction of adhesion bonds that existed between soil gieitites “ie 
in situ. Soil sampling by utilizing block-sampling techniques tends to minimize 
_ the destruction of these bonds. The formation of the adhesion bonds in granular ~ 
material can be caused either by the e consolidation of of a sand material under 
_ high effective stresses for a short period of time as shown by Lee (3), or. 
under low effective stresses for a long period of time, as shown by Malis, 
__ The process of Pitcher barrel sampling and the subsequent laboratory prepara- 
tion combined to cause disturbance. The amount of disturbance caused indivi- 


_ Block sampling into clear large-diameter plastic cylinders allowed retention - 
of a reasonably undisturbed sample. Further hand-carving allowed selection 


_ of the undisturbed portions of the large sample, usually observed away | from a 


the edges. The laboratory preparation technique where block ‘samples were 
: saturated, drained, and then further hand-carved, or frozen with liquid nitrogen 
and hand-carved, was effective in minimizing disturbance. 
Corrections to cyclic strengths obtained by tube someting, and laboratory 
_ preparation would be useful and significant. More data such as the preceding _ 

- need to be compiled to identify ranges of corrections (increase factors) that 

could be applied to © obtain n realistic strengths for natural soil deposits. Based 
on the limited data “presented here, correction factors for the dense to very ; 

er sands tested could conceivably range from 1.05-2.0, depending on the = . 


mumbers of cyclesof loading, 
- Py og the preceding ¢ examinations should not tbe. construed as a condemnation 


s have found in further 


vi that proper m maintenance and adjustment of the samplers along with careful a 
= pas of the clearance ratio of the inner tube tips provides a very good sample 
; in dense to very dense soil deposits. The samples are quite adequate f¢ for static a 


testing 2 and comparable to most tube samplers fo for dynamic testing. ber 
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By Achintya Haldar’ and Wilson H. Ts Tang, 
4 The devastating damage caused by ear earthquakes in n Anchorage, Alaska ‘a and 3 
Niigata , Japan in 1964 created an enormous amount of interest in the research 
communities to understand the behavior of soil during and immediately following 
_ the earthquakes. Consequently, voluminous research efforts were exerted. The — 
results of of these s studies started appearing i in the | literature early as 1966 (28). 
With fewe exceptions (15,16), the dynamic behavior of soil samples in the laboratory = 
were studied under uniform cyclic loading conditions (3,5,7,10,17,19,21,22,24,27— 
32). The results thus obtained were — to study the behavior of soil 


An earthquake loading pattern is extremely iregular. Ideally, to investigate 
"the behavior of soil during earthquake shaking i in the laboratory, the specimens - 
; should be tested under irregular-patterned loading like that actually generated 
during earthquakes. The use of the actual earthquake loading ir the laboratory, 
however, may not be practical in many cases for the following reasons: (1) 
_ The exact shear stress-time history of a future earthquake at a particular site 
is unknown; (2) assuming a given time history « of earthquake acceleration, the 7 
induced shear stress-time history may vary from site to site; and (3) even if. 
_ a shear stress-time history is assumed, it is very difficult and expensive to 
_ duplicate this loading on a soil specimen in the laboratory. _ ae val 
_ Moreover, in the past, a considerable number of laboratory investigations | 
is carried out under uniform cyclic loading conditions to predict the in situ 
q 
ri The behavior under any dynamic loading, including earthquake loading conditions. 
The past research results can be utilized properly only when a successful 
correlation between the two loading conditions can be found. The number of 
equivalent uniform cycles corresponding to an earthquake motion is one of 
the important parameters input to a liquefaction potential calculation. ‘Thus, 
a reliable correlation between the ae conditions is necessary in a 
_ The purpose of | this paper is to stud / and su gest a statistical correlation 
purp pap y 
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between the irregularly- patterned earthquake loading and the equivalent uniform : 
stress cycles, and compare the statistical relationship with the relationship * 
commonly used in practice. The factors involved in the process of converting _ 
= earthquake loading to equivalent uniform s stress cycles are identified. 7 
Considering the uncertainty associated with these factors (which | may not be 
negligible in many cases), the problem is studied statistically, based on the 
results published in literature and other methods. 


a Prosiem 


The general consensus is that the on: history ot an 
~ may be represented in terms of an equivalent ‘number of uniform — a 7 


cycles, N to Lee and Chan (18), betes mer 


, Tefers to that number of uniform cycles of stress intensity t,,, which a 


m ‘applied to an element of soil in the field or a sample of the same <a 
soil in the laboratory, would have the same effect in terms of the soil \ Ud 
«a Strength or deformation as if the ; actual train of irregular cyclic shear a 


‘This concept of N,, is essentially based on Miner’s (23) damage rule. The 
method proposed by. Miner is logical, simple, convenient, and has been found 
to be reasonably - satisfactory for many conditions (1). The N, , concept in soil — 
dynamics has been described in “detail in the literature (1,11, ‘18 ,25), and used 
extensively to solve practical problems (13,20,26,27,29). For details of the 
conversion procedures from irregular to equivalent uniform tel the readers — 


record and the shear stress-time- e-history record obtained from a 
in the deposit are normalized properly, i.e., when the corresponding ordinates 
of the records are divided by the maximum acceleration and the maximum 
_ shear stress intensity, respectively, they are expected to produce similar norma- — 
— lized time history records. The values of N., from both records are expected 
to be the same. Thus, in the subsequent sections, the analysis | will be based a va 
on any normalized time- history record without specifically whether 
it is from an acceleration or shear stress-time history record. ts 
The details of the conversion procedures from irregular to equivalent uniform a. 
cycles suggest that the value of N.. _ depends on the selected intensity of the 
_ uniform cyclic stresses and on the. specific soil strength curve for the site. 
_ The soil strength | curve can be described as a failure curve representing the — 
relationship between the intensity of the applied uniform loading and the number 
of cycles required to cause the specimen to fail. The intensity of equivalent — 
_ cyclic applied shear stress, S .» is usually referred to as a percentage of the = 
maximum shear stress, which in turn is defined as the shear stress” mecded 7 


When the “uniform cycles of stress intensity “concept was in solving 
problems related to earthquake excitation, some implicit assumptions were made. = 
These include the following: (1) The ground motion is uniform at all sites, 
Le €., the oo motions used for a particular site can also be used 1 for other - 
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similar sites in the same general area; (2) the stress- time —— at the depth Re 
of interest is directly proportional to the acceleration recorded at or near the ~ 
4 ground surface; and (3) for all soils, the laboratory liquefaction test data results iy 
_ can be represented by a single _ normalized curve relating stress ratio or stress 
_ level, S,, to the number of cycles causing liquefaction, 
a These assumptions have been studied extensively (1,11 18). Due to the tectonic a 
a nature of earthquake loads, the first assumption may be reasonable. The validity 
, of the second assumption depends on the soil profile through which the . 
frequency motion coming up from the bed rock. Shallow or stiff soil deposits _ 
4 will tend to maintain the same motion shapes. The time history may change 
along the depth of the deposit, depending on the soil properties involved. To 


= some of these aspects related to N., evaluation, Lee and Chan (18) a 


Seen has" been propagated. Deep soft deposits may ‘smooth out high 
fi 


considered six different soil deposits. They concluded that “the value of N,, 
computed from the s surface or near surface time history is ‘appropriate for all 


gied edz ‘hice 


1 A 50 
ods tual 3] ext t bee ,<baol 
s, 1 —Normalized | Soil- “Strength Curves 
a other depths within the soil profile.’’ They added that for routine work, it 
should not be necessary to determine ¢ Tigorously | for more than one location 
in ‘a soil profile being analyzed. Regarding the third “assumption, , although a a bed 
‘af . general trend can be established, a wide spread in the soil strength curve has ( 
been observed (see Fig. 1). The reasons for this spread have been _— 


With these assumptions, the general procedure for converting the irregularly- j 
patterned shears stress-time history to unif orm cyclic shear appears sin simple. 


However several questions could be raised. 
1. Would the value of N,, be sensitive to the shape of the soil strength — 
What stress level, S,,8 should be use used? 


me How would Neg vary with the magnitode of the earthquake? 
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4 ; and the earthquake magnitude, which one or both of the two horizontal seahhs ra 


These questions need to be to understand the significance of N., 
in studying the earthquake motion. To address all these questions adequately 


associated with the estimation of N., 

of the risk of liquefaction. Several “probabilistic models (4, 6,8,9, 13,33) lave 

been proposed to evaluate the liquefaction potential of a site, considering the — 
- 7 uncertainty associated with several parameters. The risk of liquefaction could 

be significantly different depending on the uncertainty in the parameter involved 

_ (12,14). The incorporation of the uncertainty associated with N., in some of . 

the aforementioned probabilistic models may be “desirable. This : area will be 

explored further in this paper. The estimation of the value of N,, for any a 

future earthquake, along with the uncertainty associated with that estimation, 


Sransticat 


parameter for the estimation of an equivalent number of adit cycles, N,.. a 
i | f rom the irregular time history of an earthquake. The probiem could be simplified © 
; ~ considerably if a single normalized soil strength curve could be identified for 


all soils. The writers (13) explored this area extensively. They collected from 


: _ test and simple shear test conditions relating the intensity of uniform shear . 
stresses applied to the specimen, number of cycles required to liquefy the a 
specimen, the stress conditions in the specimen during the application of the 
uniform loads, and the soil Properties. It was observed that the laboratory test 4 


‘results could be normalized in a systematic way by introducing a nondimensional 


parameter, R, which can be defined as: 4 


‘in which + = 
= the initial ambient pressure under which the sample is consolidated; ‘and 
D,= relative density of the sample. Basically, R is a measure of the resistance 
of the soil to liquefaction in terms of the shear stress intensity, tr. Some correction | 
’ factors need to be applied to R to obtain the in situ ‘soil resistance to liquefaction. . 
They are examined elsewhere (11,13). The writers observed that this soil strength — 
q parameter, R, depends on the mean grain size of the soil specimen D,,, -— 


_N,,. The relationship among R, D,,, and N,, , was proposed in the form of 

Mag De = A + Bin(n,,) + (dso) 


@ Stalisiical approach may appropriate, since some Of the parameters involve 

in the process, namely, the soil strength curve, have a considerable amount 

| 
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Parameters A, B, and Cm the coefficients of the regression equation. The 
left side of Eq. 2 is the expected or mean value of the soil strength parameter, — 

a _R, for given values of N,, and D,,.. The scatter of the experimental results 4 

about the mean regression equation as measured by Var (R |N Dso) is assumed 

to be constant regardless of the values of N,, and mi Regression equations 

- for different testing conditions (i.e., triaxial or simple shear) and for different 
failure criteria (i.e., initial 5% or 20% double amplitude strain) have been proposed 

_ (11,13). As expected, the variance of R is found not to be negligible (11), suggesting 

- considerable scatter of the test data around the mean regression line. Furthermore, - 

_ the regression equation suggests the nature of the soil strength curve, at least _ 

7 within the limits of the data considered. A typical soil strength curve thus —n 
; , - obtained can be normalized by dividing it by the value of the R parameter dae 
a corresponding to N. eq = 1 cycle (the ordinate of the ct curve expressed in = ; 
of the percentage of cyclic stress causing failure in 1 cycle). The soil strength 

. curve (similar to Eq. in can similarly be normalized. Lee and Chan (18) gee 


1 Equations between N, and Me at M > 5.0 5.0 


Soil- “strength 
curve 
(3) 
| +75% of data . 6. 162, 
Mean —59. Ss 78.6 
+75% of data 31.92 
Mean 
—75% of data 
of data 48.82 


Mean | 66.39 
of data | 


—75% of data | 169. —58.52 5. 10° 


similar types of normalized soil strength curves; they are shown in Fig. 

3 Observing the scatter of the data, Lee and Chan (18) proposed upper and lower 

75% ranges of the available data along with the mean curve as shown in Fig. — 
1. Within the limits of the available test results for the regression analysis, : 


' the normalized mean regression | curve obtained from Eq. : and that Proposed — 


“curve can be identified. The effect of the 1 mean ‘soil strength curve and the 

5 associated uncertainties in the evaluation of N., _ will be examined in the following _ 
sections, considering the normalized soil strength curves proposed by Lee and 


varies with the duration of earthquake shaking, it is expected that some kind — 
of correlation would exist between N and the earthquake magnitude. Earthquake 


_ magnitudes expressed i in ‘Richter’ s scale are considered fee this analysis. —“—s 
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Equation Number of 
— 
= 
. has been shown that larger magnitude earthquakes are associated with a longer ae . 


| and ‘Chan (18) reported values of 37 recorded 
near the ground surface. They also included 12 artijicially generated earthquakes 
ia their studies. Only earthquakes of magnitudes greater than 5.0 were considered. 
These data on N ,, and earthquake magnitude are used in the subsequent statistical 
Observing the trend of the aforementioned data, the relationship between 
__N,, (for a given stress level, Bast and soil Strength curve), and magnitude, M, 
“may be represented by the following regression equation: | 
in which A, B, and C = regression coefficients. Eq. 3 would give an expected © ¥ 
or mean an of | Ne for a given earthquake magnitude, m. The scatter of 
the data about the n mean curve is observed to be approximately constant. Thus, _ 
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2. ‘Number “Equivalent St Stress Cycles an and 


citi > hoe sit 16 oi 
4 the variance of N 9 OF Var(N,, |M = m) is assumed to be constant. The —_ 
coefficients and variance are estimated later in this ods 
Selection of Stress Level, S, .—The value of N for a given acceleration-time 
history record depends on the stress level, S, , and on the soil strength curve. Bd 
The soil-strength curve has considerable spread, as shown in Fig. 1. The choice 
of S, is primarily subjective, . based on the literature review. Intuitively, the 2 . 
particular choice of S, could be based on the degree of sensitivity of the N.. 
J versus magnitude relationship to different soil strength curves. To find a suitable 7 
value for S,, the data reported by Lee and Chan (in a personal communication 
and Ref. 18), is considered here. For S, = 65%, 75%, and 85%, and considering a 
the average, and +75% of data soil strength curves (Fig. 1), nine sets of data 
on pairs of N., and M values can be generated for each of the 69 earthquake 
Cha histories Teported by Lee (in a personal communication), and by Lee and 
Chan (18). Regression analysis is then performed on each set of om » thus = a 
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tabulated in Table 1. The regults are plotted i in Figs. 2, 3, and 4. _ is 
In the simplified procedure, Seed and Idriss (27) assumed 
They also suggested that N.. , be 10 cycles, 20 cycles, and 30 cycles corresponding 
— to the earthquake magnitudes of 7.0, 7.5, and 8.0, Tespectively. This suggested  __ 
a, relationship between No and M is plotted in Fig. 2 for comparison. Similar _ 
_ curves corresponding to S, = 75% and 85%, based on the same suggestion 
by Seed and Idriss, may be ‘derived following the general procedure (18). These _ 
_ curves are plotted in Figs. 3 and 4. The following observations can be made 
e from Figs. 2-4: (1) The effect of variation of the soil strength curves on the f 
4 N,, versus M relationship i is minimum when S, = 75%; and (2) Seed and Idriss’s Ss 
suggested relationship lies about one and a half s standard deviations below the 


a 
= 
3 
2 
a 
a 
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3.—Relationship between Number of Equivalent Stress ss Cycles Magnitude 
; _ The first observation is very interesting. The closeness of the three curves 
in Fig. 3 indicates that if the stress level chosen is 75%, then N,, calculated — 
on the basis of the mean soil strength curve will not be significantly different 
from those based on other soil strength curves. ‘Ths effect of uncertainty in- 
; the soil strength curve on the estimation of Nw can be explored further by © 
again considering Figs. 2, 3, and 4. these “three figures, N,, versus M 
-— gelationships are plotted in solid lines for the mean soil strength curve, in 
 chain-dotted lines for +75%, and in dotted lines for —75% of the soil strength — 
curve. The relative location of the N, versus M relationship curves corresponding — = 
to the mean and +75% of the soil strength curves appears to be random although 
general nature of the curves are identical. Thus, the relative location 
_ be very critical as far as the relationship between the earthquake — a 
and N, oq 18 concerned. Thus for further Eq. 5 of Table 1 may b 


| 

obtaining nine separate r 
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as an acceptable. relationship between N,. and M. ‘The. >. 
E(N,.|M = m) = 106.08 — 36.42 m + 3.312m?; m= 5. 


and the corresponding Var (N.,|M) = 29.05. When the earthquake magnitude 


is 7.0, the corresponding expected or mean value of Neg will be about 13. * 
cycles with S, = 0.75; whereas it is about 70% probability that N,, will be 
within 8 cycles and 19 cycles. This scatter is considerable. 108 
_ The second observation indicates that Seed and Idriss’s suggested relationship — 
underestimates the value of N ., for agiven magnitude earthquake. The implication 

ae this observation may not ‘be as S significant a as it ee to be. nee will 

a ship. —The data s sets used to obtain the regression equations given in Table 


FIG. 4.—Relationship between of Stress and 


_ Linclude the earthquakes of wails between 5.0 and 8. 25. In many engineering 
problems, where strong motion earthquakes are of primary concern, earthquakes 
_ with magnitudes less than 6 may not be relevant. To study the significance 
of the lower bound earthquake magnitude on the N, , versus M relationship, 
another set of data in which M is greater than 6.0 is “considered. Again, ne ; 
sets of data on pairs of N. and M values can be "generated for this reduced — 
sample of 49 earthquakes, considering S, = 0.65, 0.75, and 0.85, and the average 
P and +75% of data soil strength curves. The regression analysis given by Eq. 
2 is performed on each set of data. The results are presented in Table 2. bye 
The regression curves thus obtained can be plotted similar to Figs. 2, 3, 
: and 4. When plotted, they exhibit the same | characteristics as mentioned in 
the previous two sections and are not presented here. In Fig. 5, the N., versus — 
M relationships based on the mean soil strength curve and S, = 75% are ’ plotted i 
for the cases M = 5. 0 and M = 6.0. The proximity of the two curves taal 
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the N., versus M relationship for the mean soil 
= 15%, shown i in Fig. 3, could be used for all practical purposes. abenange Th 
(N .,) max Versus M Relationship.—In an earthquake, three acceleration-time 
history records are available at a particular site: two in the horizontal directions e 
and one in the vertical direction. In many problems, the effect of the compression 
waves produced by an earthquake traveling upward through a saturated soil - 


deposit is is probably negligible (11, 19). Thus, for a liquefaction study, the vertical 
accelerograms need not be considered (11). The value of N ,, for the horizontal 
accelerograms can be calculated in two ways: (1) By considering the accelerogram 
— contains the maximum acceleration a,,,,; or (2) by considering the 
accelerogram which gives the maximum value of 
All previous analysis made in this section are based on the first alternative. 
4 To | ) study the second : alternative, which considers the a accelerogram giving the 4s 
maximum value of N,,, the values of N., corresponding to the two horizontal 
time histories are sadibesed for each earthquake. The larger of the two values al 
identified as (N,, Then a regression analysis (Eq. 2) is performed to 
TABLE 2 —Regression Equations between and M at M = 6.0 


 E(N,,|M = m) = A + Bm + Cm? 


Soil-strength 
curve 


data points 
+75% of data | 64 91. «187. 
Mean 16 | —87. 379 88.74 
of data 
+75% of data 
of data 
+75% of data 
Mean 
—75% of data 


find the relationship between (N .,)max and the magnitude, M. Results similar 
to Table 1 considering all 69 earthquakes can be obtained. The relationship 
between (N, eq) max and M, based on the mean soil strength curve and S, = 
75%, is plotted in Fig. 5, along with two other previously mentioned curves 
representing the N.. versus M relationship when M = 5.0 and M = 6.0 for 
comparison. As expected, the (N 2g) max Versus M curve lies above the other 


two curves. However, the difference is not significant, especially in view of 
rve 
& large values of ‘Var (N., |M = m) about the regression curve indicated in 
Table 1 for all cases. Moreover, due * the logarithmic relationship of Eq. 
2, it will take a large variation in N,, to produce a significant variation in 
 : cyclic shear stress. This was also observed by Lee and Chan (18). From 
aforementioned observations, ‘it seems reasonable that since is the 


- main design input parameter in the evaluation of liquefaction potential or similar | 
a problems, it is to N ,, considering the 


— 
number 
, 5 0.75 
6 | 0.75 
q 
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and "Liquefaction Study. —Figs. 3, and 4 indicate a considerable 
_ ‘discrepancy i in the N,, versus M relationship proposed by Seed and Idriss (27) I 
and that obtained in this study. It may create > great concern among the practicing 
between N,, and M proposed here has a considerable amount of uncertainty. 
The problem may not be as serious as it appears at this stage. When the N,, 
concept is used to study the liquefaction potential, it can be | shown (13) that 
for a saturated sand deposit of, ¢.g., mean grain size D,. = 0.2 mm (0.008 _ 
in.), a 50% variation of N,, values, say from 20 cycles-30 pore would cause 
the soil-strength parameter, R, to change about 7%, according to Eq. — 
"small change in the R value for a large variation of N,, values is due to the : 
basic nature of the logarithmic selationship in Eq. 2. “A similar observation 
was made by Lee and Chan (18). I Tn a probabilistic | approach, this variability 


7 to evaluate liquefi action potential. It must be emphasized that the mean relationship — tin 
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5 —Comparison of Number of of Equivalent Stress | Cycles Versus Ma: Versus 


_ will be systematically incorporated. The large uncertainty in the estimation of 
_ N,, will undoubtedly increase the uncertainty in the R parameter, but this increase 
may not be as significant as the uncertainty in some other parameters, ¢.g., . . 
relative density or maximum 1 acceleration estimation at the s site in the liquefaction _ 
_ study (11,13). For a liquefaction study, a large discrepancy in the estimation - 7 
of N,, may not be significant; for other problems in soil dynamics a coal 
deal of attention should be given in selecting a single value relationship between 


ub 


The applicability of equivalent uniform stress cycles in soil dynamics to the 

study of soil behavior during and after an earthquake has been studied. It is 

observed that the } ectual al irregular | time histories of shear - Stresses s produced by 
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an earthquake ca can be represented by uniform amplitude cyclic shear stresses, _ = 
although there may be a pe ET of uncertainty associated with — 
them. When the intensity of uniform stresses is selected to be 75% of the a 
Maximum stress, the variation or uncertainty in the normalized soil-strength 
curve has a minimum effect on the value of the N,, versus M relationship. 7 
Thus, S, _ = 0.75 is suggested t to | be used for such conversion. . A Statistical oS 
eee between . N,, and the earthquake ‘magnitude i is proposed here, based ie 
on results obtained from “earthquake time histories recorded at or near the ground 
_ surface. Some artificially generated earthquake time histories are also included 
in this evaluation. It is also shown that N, 5 could be estimated adequately — 
by considering the component of excitation containing the peak acceleration. 
te versus M relationship proposed here is ‘somewhat different from the 
relationship s suggested by Seed and Idriss in the course of their liquefaction 
- study. It has been observed here that this discrepancy in N., values may not ; am 
a yield significant differences in estimating the soil strength in a liquefaction study. _ 
The influence of N,, values estimated by the two methods is not known for ‘ 
other problems in earthquake engineering, but this be cevaluated 
in a manner consistent with the analyses performed in those problems. wo 
This material is based upon. work partly” supported by the National Science ia 
_ Foundation under Grant Nos. PFR-8006348 and PFR 78-19680. Any opinions, — 
‘fading, and conclusions or recommendations expressed in this publication are - 
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RELATIONSHIP OF ELECTRICAL DisPERSION- 


By Scott S. Smith,’ A. M.. ASCE and K. K. Arulanandan,? M. 


A new and ‘unique has recently been for characterizing 
ss mineralogy, state of packing, and area of particle contact in fine grained pa 


soils (2,3,5,14). The technique is based on the variation with alternating current — 


frequency of the apparent dielectric constant, e’, and conductivity, o, of saturated 
clays i in the radio frequency range (10° Hz-10° Hz). This phenomena is referred — 
0 as dielectric dispersion. A dispersion curve for a saturated illite clay is shown 
in Fig. 1. With increasing frequency, e«’ is seen to decrease and o to increase. 7 
It has been shown that the principal factors influencing the dielectric dispersion | ; a 
of fine grained soils in the radio frequency are the compositional properties — ; 
of the different phases and the heterogeneous nature of the system (2,14). The = 
_ engineering properties of fine grained soils are also controlled by these factors. — 
In an effort to exploit this mutual dependency and develop a nondestructive 
- method for predicting engineering properties of soils, the writers have employed 
a simplified three element heterogeneous electrical model (Fig. 2) to represent © 
_ the complicated soil system (13). In Fig. 2, , each zone is represented by a parallel 
circuit of resistor and capacitor. Impedance of | each zone is determined by 
its dimensions (a,b,c,d) and the specific conductivity (k,,k,) and dielectric 
== (€,,€,) of the material forming the zone. The parameters of the model 
are both compositional and heterogeneous, and the model exhibits a theoretical 
- apparent dielectric and conductivity dispersion in the radio frequency range. __ 
The electrical model shown in Fig. 2(a) considers that the flow of electrical 
current through the soil can have three paths: (1) Through solution and clusters 
g particles in series; (2) through clusters in contact with each other; (3) through 
_ solution only. Fig. 2(b) is a schematic representation of the model. The geometrical 
4 parameters a, b, and c represent the fractional cross section of the three paths - 
and d represents the length of electrical current flow through clusters in path 
(1). The impedance of this model is uniquely determined by the geometrical _ 
| ‘parameters and the dielectric constants (€,,¢,) and conductivities (k,,k,) of 
the clusters and solution, respectively. In the case of soils saturated with water, — 
= WwW 
of Civ. Engrg., Univ. of California, Davis, Calif 95616. 
- Note.—Discussion open until October 1, 1981. To extend the closing date one month, 4 
_ a written request must be filed with the Manager of Technical and Professional Publications, om 
ASCE. Manuscript was submitted for review for possible publication on May 27, 1980. 


_ _ This paper is part of the Journal of the Geotechnical Engineering Division, = 7 
is the American Society of Civil Engineers, — Vol. 107, No. GTS, May, 1981 
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< FIG. 2.—Heterogeneous Electrical Model: (a) Representation of Current Path through = 


“the pore fluid. The apparent dielectric constant, €/,,, and the apparent conductivity, ; 
o,,, Of the model have been evaluated as a function of frequency by elementary 


+ 
d(l-d)Ss 


l-d 


= the capacitance of the unit capacitor in vacuum, i.e., 0.0885 x 10-7 
farads; k, and k, are given in 1 mho/cm; w = the angular frequency i insec™’; 
_ These theoretical equations can be fitted to experimental frequency dispersions aan 
of the apparent dielectric constant and conductivity of saturated fine grained _ 
soils by computer optimization of the geometrical (a,b,c,d) and compositional — 
parameters. By systematically varying the structure determining 
factors within the fine grained soils, the ability of the parameters of the model © 
© characterize the ¢ composition and heterogeneous nature of the soil system — og 
has previously been evaluated (2,14). The c parameter was found to be a measure 
of the size and distribution (tortuosity) of the intercluster pores and the intercluster 
_ void ratio. The parameter is very sensitive to intercluster void ratio and intercluster — @ 


Ss and independent of compositional variations - Model parameter b is a 
measure of the contact area between the clusters. In the case of untreated 
- fine grained soils, this parameter is always “small. The term ¢, is mainly an 

_ average measure of the type of clay mineral present in ‘the soil and the intracluster 
r void ratio. The distinctive €, values obtained for different wet minerals is a 
result of each mineral’s characteristic water retention. The erey is thus 


the assumption that the solid phase i is composed of ag aggregates of primary particles P- 
as opposed to the assumption that e€, is the dielectric constant of the particle. a ; 
_ In such a case, €, should be independent of water content and fabric. The 
heterogeneous model is, therefore, in agreement with recent developments and = 
results reported in soil technology. These studies have shown that various particle _ 
and particle group arrangements exist in soil. Yong and Sheeran (15) recently 
» presented a scheme of fabric classification based upon the interaction of groupings 
§ particles. These fabric units are termed domains, clusters, and peds. The 
soil response behavior involves the interaction of fabric units of various sizes. 
Barden (5) covered the current evidence supporting a similarity in the bulk 


; behavior of clay and sand. This concept is predicated on the belief that groups a 


4 


4 


of particles rather than individual particles interact in clays. Ga ow 
The fundamental nature of the electrical dispersion in fine gre grained soils and 
the numerical parameters b, c, and ¢€, of the heterogeneous model present an F 
effective tool for studying the engineering behavior of soils in a fundamental a 
-quantitive manner. This paper presents a the results and 


ELECTRICAL DISPERSION = 
q | 
4 


conclusions of three such studies. The studies are considered sigutyand 
included: (1) The evaluation of the heterogeneous model’s ability to monitor _ 
structural changes during the hydration of soil cement; (2) relationship between _ 

model parameter and swell ability. between model 


ExPerRIMENTAL ProcepuRrE 


basic equipment and procedures followed in these studies have been 
reviewed in detail in previous publications (2,5,14). The electrical dispersion — 
measurements of the soils were performed with a type 250 RX meter (Booton © 
Radio Corporation, | Division of Hewlett- Packard, Rockaway, N.J.). The instru- — 

“ment is essentially a Schering bridge, with oscillator, amplifier detector, and 
null-indicator designed to measure equivalent parallel conductance and capaci- _ 
tance in the range 0.0 mhos-0.067 mhos at frequencies of 0.5 megacycles-2. 50 wy 
megacycles. The soils to be measured were loaded into an electrical cell and — 

' connected to the electrical bridge. A novel cell design and method of connecting — 


effects | that arise at radio frequencies (13). 4a, 
_ A digital computer optimization program (1,10) was employed to analyze the 


7 emperical dispersion curves. The optimization program determines the values 


of the geometrical and compositional parameters of the model such that the 
_ frequency dispersion curves of the apparent dielectric constant and conductivity 


calculated from the theoretical Eqs. 1 and 2 fit the empirical results within 

- acceptable error limits. The value of ¢, was kept constant at 79, the dielectric b 

- constant of water, during all computer runs. The initial estimates of the remaining 

_ parameter values were determined following a rational and standardized proce- 
dure. One might assume that a set of seemingly random model parameters 


could: fit almost any set of experimental dispersion curves. However, in a 


_ many y applicati ons of this model, ‘the optimized model parameter values have 
given reasonable results from our knowledge of soil technology and from the 


_ There are, in fact, six random parameters ¢, a k,, € s» and k, involved 


coal to zero and a is ‘dependent | on c, from the result a + b +e = 1. At 
first sight, it might | seem that the choice | of six seemingly random ‘parameters 
- would be enough to fit almost any set of curves and that the physical meaning 
of the result is doubtful. It is important to realize, however, that only the 
_ three independent geometrical parameters are genuine unknowns. The respective 
conductivisies k, and keys of the solid and interstitial solution and t the dielectric 


; fr rom om the point of view w of physical s science. One would expect that the conductivity 
_ of the interstitial solution to be of the same order as the conductivity of the ; 
solid pore fluid extract; the conductivity of the solid should be of the same 
_ order as the isoconductivity value (see Fig. 18 of Ref. 2); the dielectric constant 
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- computer program, the results presented in the following | show that the require- 


Sructurat Cxances or Son Cement Durinc 


| The ability of a ates to stabilize soils i is often a function of its capacity 
to affect alterations in elena the soil. Mitchell and El Jack (9) have — 
utilized the electron microscope to investigate changes in fabric and composition 
of soil cement during changes were reflected as an increase 
in strength. Accordingly, the heterogeneous model was applied to one of the © 
soil cements studied by Mitchell and E El Jack to evaluate the model’s a 


4 
to monitor these Structural changes. 


Pr 
CURING TIMES 
032 Days 
16 Days 
Days 
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UF Soil Cement at Varying Times after Mixing and for Pure Kaolinite Hydrite UF — 7 


=~ 


The soil cement was a commercially available kaolinite clay ‘(Hydrite UF) 


with 24% Type II portland cement by weight. The dry mixture of kaolinite 
and cement was mixed with distilled water to a saturated water content of 
140%. A high water content was necessary to assure that no air voids would 
result when the soil cement was loaded into the electrical measuring cells. _ , 
_ Electrical measurements were performed at different time in’ervals as a 


: proceeded: 1 day, 2 days, 8 days, 16 days and 32 days. The electrical cells 
were stored under 100% humidity conditions between readings. The measured 
dielectric and conductivity dispersions for these curing periods are given 
ie an indication of initial soil structure at time zero, electrical dispersion we 
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measurements were ads oceans oni on pure kaolinite (Hydrite UF) ‘mixed with 
distilled water to a water content of 140%. It was necessary to use pure Kaolinite 
_ instead of the kaolinite- cement mixture because ‘the electrical properties of ye 
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—Experimental Conductivity Dispersion Curves for Kaolinite UF Soil 


Cement at Varying = 
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—— Optinized Model Fit 
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FIG. 5.—Curve Fitting of Theoretical Electrical Model Dispersion “ Experimental al 
we Dispersion of Kaolinite Hydrite UF Soil Cement for 32-Day Curing Period = a 
| soil cement during the first day of curing were found to change too rapidly 
for the 1 measurement of an electrical dispersion. Although the kaolinite 
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isn not an exact measure of the soil cements’ initial structure, it does give an 
picture of the particulate structure of the: soil cement at initial 


ite hic 
a It is seen from Fig. 3 that as curing time increases, the magnitude of dielectric ~ 
dispersion increases. The conductivity dispersion curves plotted in Fig. 4 decline a 
with reference to the conductivity scale with increasing curing time. _ aio 
The curve fitting computer program was applied to the experimental curves =~ 
as previously outlined. A typical curve fitting result is presented in Fig. 5 for 
i the case of | the 32-day curing , period. The associated ‘optimized values of the 
parameters of the 1 model are in Table. for all periods. 


TABLE ABLE 1.—Optimized Mo Model Parameter Values during Hydration ot Coment 
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l-day 
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-day 
hydration 
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hydration . 55 | 0. |1.5 x 107% 


¥ _ “Model parameter €, was kept constant at 79.0 during optimization. 


The vi value of ‘the c parameter is seen to decrease with increasing hydration 


through the soil cement due to the interlocking of particles. This is in agreement 
_ with the observed decrease in the coefficient of permeability of clay soils upon 
the addition and curing of cement (11). 
_ The optimized b parameter of the soil-cement samples i is seen to be ‘significantly 
i greater than that of the pure kaolinite; indicating that the cement and clay 
- particles i in the soil cement mixture are no longer discrete. The area of contact 


between solids is quite large indicating the fabric is more of a solid network. iq 
. The relationship between 1 the b b pensapater and the curing time is plotted ed in 
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FIG. 6.— -Relationship between 5 Parameter of Electrical Model and Curing <a 7 


ok 6. The major increase in ‘- takes place during the first day of curing. 7 
After the third day, the rate of increase of b is small but measureable, and 
is seen to be almost constant. The large | increase in 5 during the first day 
can be ascribed to the hydration of the cement. Gel resulting from the hydration 
= the discrete soil particles into a more solid network. The small but 
constant rate of increase of b after the third day may be due to secondary 7 
" - cementing compounds resulting from the reaction of | the clay particles with 
the lime liberated during hydration. 
These results suggest the following sequence of compositional and fabric 
changes during the curing of the kaolinite Hydrite UF soil cement studied. 
_ The major portion of the hydration of the cement particles occurs during the 
first day of curing. The resulting gel interlocks the discrete soil particles into 
amore solid network. A ear cementing compound results from the reaction 
of the clay 
contribution of these ps oa ws bonds is much smaller ‘than that of the het + 
and proceeds at an almost constant rate. It would appear that the initial strength _ 
of the gel is quite low, but with time, the gel recrystalizes with a resulting 


increase in strength. These changes in composition and - are reflected 
terms of increases in strength of the soilcement. = 

_ The preceding description agrees qualitatively with the com sompustionsbe and fabric — 

7 changes suggested by Mitchell and El Jack. The rate of gel formation is, however, 
not in agreement. These macroscopic model results indicate that the hydration - 

of the cement particles is completed in the initial stage of curing, whereas 
_ Mitchell and El Jack conclude that the hydration ‘continues for weeks. It may 

: be that these differences i in interpretation are due to the fact thet the kaolinite 


_ kaolinite soil cement stu studied by Mitchell and El Jack was compacted at a water 4 


content of of 34%. * 4 
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the | type ype and amount sale ade nalentetinn present. Each m mineral ty type: has i its own o 
inherent capacity for water absorption. In the case of fine grained soils, the - 
€, parameter of the heterogeneous model is mainly an average measure of 
the type of clay mineral present. The term, €,, represents the dielectric a 
of a wet soil cluster. Its value should be between 4 (the dielectric | constant 
:., of adry silicate mineral) and 79 (the dielectric constant of pore water) depending - 

on the proportions of the two phases in the cluster and, thus, the water retention — - 

7 characteristics of the clay minerals present. This observation has guoviouliy) 
led to the investigation of the relationship between the dielectric constant of 
the wet soil cluster, €,, and the swell potential of of a soil « 


solidation, 


after con- 
as a per- 7 


(dispersed) 4, 65. 1 0.90 


“Model Parameter €, was kept constant at 79.0 during optimization. — een tintin, 
Natural soils from throughout the state of California; supplied by the state of California - 


al 
Division of Highways, 


= “Natural soil from Las Vegas; supplied by Federal Housing Administration. - 


_ Fourteen soils were studied in all in the most recent investigation (7). Ten 


The remainder of the soils are commercially available. anitenp: 


_ Each soil tested was mixed with distilled water consolidated from an initial 
Be consistency under a | kg/ cm’ load. A cylindrical specimen 1 in. (25.4. 


= ae and 1.4 in. (35.6 am) in ¢  Gomater was allowed to dry in a 50% 7 


GIs. ELECTRICAL DISPERSION. 
' alues for Swell Potential Soil Samples (7) | 
Optimized Model Parameters” 
39. 0.76 | 0.02 | 0.22 | 0.30 | 10. | 0.0006 | 0.0008 
106. 0.61 | 0.01 | 0.38 | 0.99 | 36. | 0.0053 | 0.0007 
ae 0.68 | 0.01 | 0.31 | 0.72 | 18. | 0.0018 0.0013 
0.87 | 0.01 | 0.12 | 0.69 | 41. | 0.0028 0.0008 
0.72 | 0.01 | 0.27 | 0.88 | 18. | 0.0036 | 0.0015 
32. 0.66 | 0.01 | 0.33 | 0.75 | 9. | 0.0007 | 0.0005 
© 0.55 | 0.03 | 0.42 | 0.94 | 8. | 0.0015 | 0.0005 
0.50 | 0.03 | 0.47 | 0.98 | 23. 0.0060 0.0007 
| 0.49 | 0.01 | 0.50 | 0.98 | 26. | 0.0046 0.006 
0.46 | 0.03 | 0.51 | 0.98 | 13. | 0.0041 0.0006 
114. 0.67 | 0.02 | 0.31 | 0.88 | 25. | 0.0007 | — 
(flocculated) 0.59 | 0.02 | 0.39 | 0.90 | 26. | 0.0036 
001 1009! 099 | 48 | 0.0076 0.0005 
; i were sampled from throughout the state. One of the soils was from =. a 


moisture room. The dried soil st specimen was then nesiatiy sontiat and allowed = 
. swell without surcharge. Water was introduced by placing the confined sample 
ona continually saturated porous stone. Thus, the definition of swelling potential 
in this investigation is ‘‘the percentage swell of a laterally confined sample 
oe soaking without surcharge after being consolidated from a slurry a 
under a load of 1 kg /cm* and dried in a 50% moisture room. ” Ra 

dispersion measurements were performed on yn each soil in the consolidated state 
before drying. ‘The curve fitting computer program was then applied to the 
- resulting empirical dispersion curves and the heterogeneous model parameters 
were calculated. Table 2 presents the optimized model parameter values and 
the final water content after consolidation for each of the soils. 
optimized values of €, are are plotted ve versus the percentage swell for each 

soil in Fig. 7. A linear curve has been drawn through the data points. A good 

correlation between increasing dielectric constant of the wet soil cluster and _ 
_ increasing percentage swell was obtained for these various soil | types. | It ae aa 
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Mite 


Flocculated 
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DIELECTRIC CONSTANT OF WET SOIL 


be concluded from the results obtained from these soils that €,isa measure 


of a soil’s water absorption characteristics and may prove rove useful for the evaluation — 
‘ The hydraulic permeability of cohesionless soils is reasonably well represented — 


- the Kozeny-Carmen equation (8). This equation predicts that for a given 
soil, the hydraulic permeability i isa a function of the void ratio 
equation fails to predict measured flow rates in clays. Olsen (12) nis shown 
_ that the concept of ‘unequal pore sizes’’ can explain these discrepancies. He _ 
suggested that the primary particles in a clay mass are arranged in clusters 
and the the total pany of the soil system is distributed among intercluster and ol 
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4 _ pores) will be considerably larger than those within the clusters between the _ =a 
- primary particles. By deriving the relationship for flow rates through a cluster _ 
model similar to that of Fig. 2 and assuming that the fluid only flows through © = 
& large intercluster pores around and between the clusters, he was able to 

explain the unique hydraulic flow characteristics of clays. 
The c parameter of the heterogeneous model represents the fraction of the 
electrical current passes through the intercluster solution. It has been found 


is established, [Olsen (12)], a limited experimental program was conducted to 

_ investigate the relationship between the c parameter and the hydraulic permeability _ 

of saturated fine grained soils. $= 

Ss Seven clay soils prepared from commercially available clay minerals were 
five illite (Grundite); one kaolinite (Hydrite-R); one montmorillonite 


in centi- | 
meters 


| 
@) 


tonite | 457 | 0.0010 


_ “Model Parameter €, was kept constan! at 79.0 during hydration. 


(Woke Clay). The illite and kaolinite soils were first made homoionic to a particular 

- cation by mixing a very” dilute suspension at 0.1 N and allowing the particle - 

to settle out. This was repeated several times to assure thorough cation a 

- replacement by the law of mass action. The final desired pore fluid concentration — i, 

i obtained by decanting the required portion of the solution left after the 
4 


particles had settled, adding distilled water, and thoroughly mixing again. The © 
montmorillonite was treated in a similar fashion. However, due to its strong ' 
colloidal characteristics, a centrifuge was employed to accelerate the particle 
 sutiling process. The chloride form of each cation was used. The kaolinite and 
montmorillonite clay were made homoionic to Na*. Three of the illite clay 
soils were made homoionic to Na*, and the remaining two to K* and Li* _— 
bi respectively. The illite clay soils were all prepared at an approximate concentration « ‘ 
After the desired | ‘pore » fluid electrolyte concentration had been obtained, ul 


4 Material jas a per-| per cen- Cation 
description | centage | timeter | type per sec €, 
Tite | 61 | 0.0014 | Na* | 7.7x1 9. | 0.0017 | 0.0083 
Iilite 0.0009 | Na* | 7.5 x |0.65 | 0.01 | 0.34 | 0.58 | 19. | 0.0011 | 0.0006 
62 | 0.0014 Na* | 1.1 x 10°” | 0.64 | 0.01 | 0.36 | 0.65 | 21. | 0.0015 | 0.0008 
Be 0.0011 K* | 1.1 x 1077 | 0.62 | 0.02 | 0.36 | 0.52 | 13. | 0.0012 | 0.0008) 
Illite | Lit | 13x 10-7 0.01 |035| 0.49 | 14. | 0.0011 | 0.0008 
sKaolinite | 60 | 0.0002 | Na* | 3.0x 10~’ |0.56| 0.01 | 0.43 | 0.43 | 10. | 0.0004 0.0002 “i 
Na* | 2.0x | 0.82 | 0.02 | 0.16 | 0.93 | 63. | 0.0043 | 0.0005 


soils: were consolidated from a slurry consistency under a load of 1 1 kg / cm? <2, “ 

fter the soils were through consolidating, they were leached by applying an 

vi 8-psi vacuum to the bottom of the sample and introducing an excess of solution 
: 3 to the top. The added solution contained the same electrolyte type and wes 
| ie Fae the same concentration as the pore fluid. By monitoring the volume 
of water leached through over a period of time, the hydraulic permeabilities 
”,- the samples were determined. Each soil was then loaded into the electrical — 
oat and the electrical measurements were performed and the model parameters 

_ The computed model parameters are presented in Table 3 along with other | 

pertinent soil data. A plot of the c parameter versus the hydraulic permeability ¢ 
_ is presented in Fig. 8. A best fit line has been drawn through the data points. _ 


; - from ‘these results that the < c parameter may offer a numerical evaluation of 


Y | 
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FIG. 8.— 


the size and distribution of intercluster pores and an | effective tool in the — 
measurement and study of the hydraulic permeability of fine grained soils. 
_ In an effort to develop a nondestructive method for studying and predicting © 
_ engineering properties of fine grained soils, the writers have employed a simplified — & sp 
three element electrical model to represent the complicated soil-water system. wr 
by model exhibits a theoretical apparent dielectric and conductivity dispersion | 
in the radio frequency range. The parameters of the model are both compositional Ca 
heterogeneous (structural). For a particular soil system, numerical values 
“1 of the model parameters can be determined by computer curve fitting the = < 


“equations of the model to the experimentally determined radio 


- frequency apparent dielectric and conductivity dispersion curves of the soil. ee 
Three of the (b, c, and €,) in have 


the number oF data points is limited, there appears to be good agreement 
4 


fine 4 
a soils. The c parameter was found to be a measure of the size ae 
_ distribution (tortuosity) of the intercluster pores. Model parameter 5 is a measure q 
of the contact area between the clusters. The term €, was found to represent 
j the dielectric constant of the wet soil cluster and depends mainly on the water 
retention characteristics of the clay minerals present. It - therefore, mainly 
2 an average measure of the type of clay mineral present. ts ‘ 
_ The fundamental nature of electrical dispersion in fine grained soils and the 
:. numerical parameters b, c, and €, of the model present an effective nondestructive _ 
| tool for studying the engineering behavior of soils i ina a funda fundamental 


structural changes during the hydration of soil cement; (2) the relationship between 


4 The three studies included the evaluation of: (1) The models ability to — 
model parameter €, and swell potential; the between 


ie _Mitchell and El. Jack (9) have utilized the ane microscope sto to investigate 


"changes in fabric | during hydration of soil cement. The heterogeneous n model 


was applied to one of the soil cements studied by Mitchell and El Jack (Kaolinite 
_ clay-Hydrite UF). The c model parameter was found to decrease with the hydration 
of the soil cement. This reflects a decrease in the unrestricted pore paths through 
the soil cement. The 5 parameter for the soil cement was found to be significantly — 
greater than that of the pure soil. This indicates that the cement and clay particles 
cr an 
in the soil-cement mixture are no longer discrete _ and have formed more of 
a solid network. The changes in structure of the soil cement, as monitored 
4 the model parameters, were found to agree qualitatively with the changes - 
by Mitchell andElJack 
‘The relationship between model parameter €, and swell potential has been 


"previously investigated using 11 ‘natural and ‘three artificial soils. A good 


Geet between percentage swell and €, was obtained. It can be concluded ares, 
from the results obtained from these soils that €, is a measure of a soil’s water - 
absorption characteristics and may prove useful for the evaluation of Swelling: 
= \ limited experimental Program was conducted | te investigate the relationship J = 
between the c parameter of the model and the hydraulic permeability of saturated 
fine grained soils. The value of the c parameter was found to increase with a, 
the value of the hydraulic permeability. This limited data suggests that the 
: parameter may offer an effective tool in t the measurement | on hg of the 


is gratefi ully 2 acknowledged. 


K., and Mitra, Ss. K., “Soil Characterization by Use of Blectrical 


\ 
q 
4 
| | 
| a 
fFelationship between electrica and mechanical properties Of sous, supporte 
72.The support J] 


Network,” J 
4 Nov., 1970, pp. 480-485. 
Arulanandan, and Smith, S. S., 


Frequency Electrical international Symposium on Soil Structure, Aug., 
4. Arulanandan, K., Smith, S. S., and Linkhart, T. A., ‘“‘Characterization of Clays by ll 
Electrical Methods, ” Proceedings of the Sixth Annual Meeting of Clays and Clay ~ 
5. Arulanandan, K., Smith, S. S., and Spiegler, K. S., “‘Radio Frequency Properties 7 
— Polyelectroyte Systems,” Proceedings, NATO ‘Advanced Study Institute, Forges- 
Barden, L., “The Influence of Structure on Deformation and Failure in Clay Soil,” 
Geotechnique, Vol. 22, No. 1, Mar., 1972, pp. 159-163, 
7 % ‘Basu, R., ‘‘Identification and Prediction of Swell of Expansive Earth Materials,’ 
a thesis presented to the University of California, at Davis, Calif., in Jane, 1972, in i 
partial fulfillment for the degree of Master of Science. 
on 8. Carmen, P. C., Flow of Gases through Porous Media, Academic Press, New York, 
9 9. Mitchell, J. K., and El Jack, S. A., “‘The Structure of Clay Cement and Its Formation,’ 
~~ Institute: of Transportation and Traffic Engineering, University of California, Berkeley, 
10. Muat, R. W., ““Driving- Point C Characteristic Approximation with R. ‘. Network,” thesis 
a to the University of California, Davis, Calif., in 1968, in partial fulfillment — 
of the requirements for the degree of Master of Science. 
. Nussbaum, P. J., and Colley, B. E., ‘Dam Construction and Facing with Soil-Cement,”” 
and Development Bulletin, Portland Cement Association, 1971. 
12. Olsen, H. W., ‘‘Hydraulic Flow Through Saturated Clays,” Proceedings, 9th National 


13. Sachs, S. B., and Spiegler, K. S., “Radiofrequency Measurements of Porous Conductive . 
Plugs, Ion- Exchange Resin- Solution System,”’ Journal of Physical Chemistry, Vol. 


14. Smith, S. S., “Soil Characterization by Radio Frequency Electrical Dispersion,” 


dissertation presented to the University of California, Davis, Calif., in 1971, in partial 


ie _ fulfillment of the requirements for the degree of Doctor of Philosophy. 
15. Yong, R. N., and Sheeran, D. E., ‘‘Fabric Unit Interaction and Soil Behavior,” 
Proceedings, International Symposium on Soil Structure, Aug., 1973, pp. 175-183. ol 
The following symbols areusedinthispaper; 
4 b,c, d= geometrical parameters of electrical model; Ins bed 
conductivity parameters of electrical model; a 
= dielectric constant parameters of electrical model; 
= theoretical apparent dielectric constant from electrical model; a 
a = 
= conductivity (experimentally measured); 
4 theoretical conductivity from electrical model; and 


o = angular frequency. 
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LIQUEFACTION ANALYSIS. FOR MULTIDIRECTIONAL 


ese considered only one horizontal component of ground shaking (1,6,11,12). 
_ Usually, the horizontally-layered ground is modeled as a one-dimensional system 
gg propagates the horizontal base acceleration as a shear wave in the vertical __ 
directional. In reality, the level ‘ground i is subjected | to ‘three 
7 f demonstrated (13, 14), that the liquefa action potential under the simultaneous action 
of both horizontal components of earthquake is higher than the liquefaction 
r potential under one component of earthquake. Also, the character of the surface © 
response may be significantly influenced by the interaction between the two. 


_ horizontal components 2 and the vertical component of earthquake. 


_ A new method of analysis i is presented here for computing the response © and 
the liquefaction potential of horizontally-layered ground subjected to three 
components of earthquake base acceleration. This study and the proposed method 
of analysis are the extensions of the method of ‘analysis presented recently 7 
in Ref. 6, which only considered the response to one horizontal component. = 
In the following sections, the theoretical basis sof 1 the proposed method of adie * 
is described and the main aspects of the response to multidirectional shaking 
are demonstrated through a series of case studies. These case studies were 7 
performed by using a new computer program LASS-III (7) which was developed | 
‘The horizontally- layered is divided of ‘ ‘layer elements.” 
Each layer element includes the portion of the soil bounded by two adjacent 
‘horizontal ‘ ‘nodal planes.’’ The whole system is subjected to three components 
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“tally- -layered | ground a are e assumed to result from the vertical propagation of the 
base motions. Thus, the response of the system is the consequence of the , 
4 vertical propagation of two horizontal shear waves and a vertical compression - 
—vave.T Thus, at each nodal plane three components of displacements are possible; be 
- uy, and u 2 and y-axes are in the horizontal plane and the z-axis is 
the vertical axis). All displacements are assumed to be only functions of z 
_ The portion of the ground below the water table is modeled as a saturated 
porous deformable solid. The general formulation for dynamic finite element 
analysis ¢ of such a two-phase medium has been given in earlier publications _ 
(2,3,5 6). In such a a ‘formulation, the displacements of the pore fluid 
- to the porous solid a are also considered as independent variables. In the horizon- — 
_tally-layered ground, the properties of the material and the motion are assumed © 
to be constant along horizontal planes. Thus, the pore water can only migrate 
ey, aes the vertical displacement of the pore water relative to the solid, — 
denoted by w, isthe only fluid variable. j= 
‘The response of the system is described in terms of the nodal plane displacement | 7 
degrees- -of-freedom and each nodal plane has degrees-of-freedom: three 
components of the solid displacement, u,, u,, u,, and the displacement of 
pore water relative to solid, w. The nodal pi remain horizontal and el 
Be displacements. Consequently, four strains a are possible: the vertical normal | 
strain, €,, the two horizontal shear strains, y,, . and Yys and the volumetric 
strain of the pore water, The Strain- ~displacement relations are as 


— 


.~ horizontal shear stresses 7, andr, ye the pore water pressure, IT. The increments 
of stresses and strains are related ome the material model which will be 
a For a typical layer element, the components of displacement are assumed — 
to vary linearly between the | nodal planes. With this assumption on variation — 
of displacements for each | element, the stiffness matrix, mass matrix, and 
: _ “dissipation resistance’’ matrix can be computed. The matrices for the system 
result from direct assembly of the element matrices. A complete development — 


; a theory is reported elsewhere (2, 3), and will not be repeated here. — vanpat ot 


LM. M,. hy Ky Ky 


aa in which u = witli of the displacements of the solid portion; w = vector 


@ the eascenuaiee’ of the pore water with respect to solid; M, = mass matrix 
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of base acceleration: two horizontal components and 4 verucal Component. In = 
ty 
| 


af 


the bulk of soll; mass. matrix for the ‘fluid M. 


are inversely proportional to ‘0 the coefficient of permeability; K,, = tangent 
3 stiffness matrix for the bulk of saturated soil; K, = stiffness matrix for the | ye 
fluid portion; K, = stiffness matrix representing the coupling between solid 
and fluid portions; ii, = vector of three components of time history o of base 
acceleration; R,, R, forces computed from stresses; and 
on = a matrix with three columns, each containing ones at degrees- of-freedom | 
_ in one direction and zeros at all other degrees-of-freedom. nee.) 
_ Eg. 2 is a nonlinear system of equations. The source of aubtinctaisy | is through 
the submatrix, K,,, which contains the nonlinear material behavior of the ~ a 
& granular skeleton. This nonlinear matrix equation can be written symbolically — 
a(t) + D a(t) Au(t) = - ML AD) xoqe od (3) 


yg which At = the time step. This nonlinear matrix equation of motion is directly 


‘Two types of e energy absorption mechanisms or “equivalent damping effects 

are present in the analytical model used here. A hysteretic damping is introduced - 

ome the system through the material model which will be examined in the next | 

- Section. The magnitude of energy absorbed with this damping mechanism is 

- dependent on the history of shear strains within each element. Another energy 

‘ absorption mechanism is provided through the dissipation resistance matrix which 
represents dissipation of pore water pressure and flow of pore water 
This damping mechanism depends on the coefficient of permeabili- 


in time bya a time-marching scheme, as described in Ref. 3. 


’ an was ‘a that the dissipative damping is more effective at higher 
: frequencies, but the hysteretic damping is uniformly effective across the possible — 
_ frequency Tange. These two damping | mechanisms do not require any additional — 


_ In the horizontally-layered ground, the state of stress consists of two shear 
stresses, T,, 7,, the effective stress, o’, and the pore pressure, Il. The 

of the effective stress and pore pressure is appropriately taken into account 

through the coupled discrete equations of motion of Eq. 2. In addition, the — 

dissipation of the pore pressure which is governed by D’Arcy flow law is also 

_ taken into account in the equations of motion. Thus, the main purpose of a 

ae model is to establish rules governing the interrelationship between the 
shear stress, the effective stress, and the corresponding strains. It is sssumed 
‘that the in situ state of stress is isotropic (lateral stresses are equal to the 
vertical stress). Furthermore, it is assumed that the vertical effective stress 
changes occur elastically. This is a reasonable assumption since even a 
vertical acceleration, although the total vertical stress changes, the effective 
stress” varies very lit little, as be shown in a later ‘Section. Moreover, 


tv 
needs to be evaluated in a later pubiation. “wihss bens 
_ Shear Stress-Strain Relations. 
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components of shear stress on vary arbitrarily. However, the stress-strain 
relations in the two directions are not independent. The incremental stress-strain _ 
relations are determined by an elastoplastic model for cyclic behavior of sands — . 
“4 under biaxial shear, which is described in detail in Appendix I. The basic aspects 
4 of the material model are briefly examined here. The material model used in 7 
_ this paper is the extension of an earlier material model (6) to include two | 
compensate of shear stresses, and is a special version of a more general material — 
model, which | has been developed to ‘Tepresent the cyclic behavior of sands 
% A yield surface has been defined and described in detail in Appendix I. The © 
yield surface has a conical shape and initially its radius is zero and its axis _ 
7 coincides with the effective stress axis. As the shear stresses are applied, the ‘a 
" radius of the yield surface increases (isotropic hardening), and its axis rotates 
_ while the apex remains at the origin of the stress space (kinematic hardening). * 
_ The rules governing the isotropic and kinematic hardening are described in 
. Appendix I. The combination of isotropic and kinematic hardening allows the 
representation of hysteresis loops in the stress-strain relations. 
All the stress changes within the yield surface result in elastic shear srsies _ 
g no changes in the effective stress other than those resulting from dissipation 
of pore pressures. All other stress changes result in an elastic shear strain, a 


q _d-y‘,a plastic shear strain, and cause changes in the ef: fective stress according 


_ The equation for elastic shear straincan be written 


_ The elastic shear strain increments are related to the shear stress increments — = 


by the elastic shear modulus, G, which is defined as on initial Slope (oft the 
== Portion of the shear Stress- curve 


ay fdr, 


= a‘ “plastic modulus,” defined in the following equation 
in which {n}7 = {n,,n, Te ee unit outward normal to the yield surface. aa 


In a unidirectional ca case (, = Ayn n, = ( 0), ‘Eq. 6 relation 


Ite can be seen that H = the tangent modulus of shear stress-plastic shear strain — 
_ relation, which has been expressed through a hyperpolic function. Such a 
hyperbolic stress-strain relation has been proposed by Kondner and Zelasko — 
(10) and has been used earlier (4,8). The only difference with the caer 


applications of — stress-strain curves is that here it is expressed in 


n, 


J 
/ 
» 
4 
___The increments of plastic shear strain are determined from the normality _ | 
| 
| | 
— 


terms of plastic shear strain rather than the total cian strain. This | 
is given in the equation and shown in 


It can be seen that M = tan is the maximum value of. s/o, which occur 


at failure and satisfies the well known Mohr-Coulomb failure criterion ac iu 


‘The initial slope of the shear stress ‘divided by the effective stress versus 
the plastic shear strain (Eq. 8) is H,, which is called the “initial ‘Plastic modulus.” 


the expression, 

ty 


the plastic modulus decreases, zero at failure : = Mo 
- itt is important to note that with the proposed model the an ‘slope of 
a unidirectional shear stress-strain curve contains both elastic and da io 
- contributions, since the | plastic deformations occur from the start of application 
of shear stresses. The elastic shear modulus is the slope of the unloading — 7 


of the stress-strain curve. e. The initial Plastic | modulus" is determined from | the 


in which which G, = = att = 0. 


 & 
M=tand...... 
4 
IG. 1. *-Strain Curve 
(10) 
| 
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stress-strain relations under biaxial shear are given ‘the 
following equation. The details of the derivation can be found one igae I 5 


a, 


acral +4 


6 this equation, n, and n, = the components of the unit ‘outwardly normal > 7 
Z the yield surface. tt can be | seen n that | in in unidirection (a = I,n = = 0), the a 
= d 


garraiig) 
Wal 


and a modulus. The equivalent modulus is a decreasing function 


of the resultant shear stress, approaching zero at failure shear stress. _ y pl : 
 ptectv Stresses.—Two factors contribute to the changes in the s effective 
_ stress. The pore pressure dissipation and redistribution during the earthquake 
_ causes changes in the effective stress. This effect is correctly included in the 
system equation of motion (Eq. 2). Depending on the coefficient of permeability - 
of sand, this effect can significantly influence the seismic c response of the layered 
ground. The second contributing factor to reduction of effective stresses is 
; the volume change tendency of sand under shear stresses. This effect is modeled - 
by prescribing an undrained stress path in the shear stress versus effective 
stress plane. This approach is similar to the method of determining effective — Y 


_ An undrained condition i is eet at each instant when the effective s stresses _ 


| 
&g 
quantity composed of 


zi ‘time step to the next, when the coupled equations of motion are solved. The me : 
_ undrained effective stress path, similar to Ref. 6, is assumed to be a quarter — 

; of an ellipse, as shown in Fig. 2. The equation for the effective stress i 

: 


a 
. The two parameters, o/ and r, completely define the effective stress path. 


The effective isotropic pressure, Mine defines the position of the effective stress 
path. The material parameter, the ‘Tatio of vertical to horizontal axis of 
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«FIG. 3. }.—Relationship be between Material Param Parameter, and Relative ve Density (6) 


_ of experiments reported in the literature. In general, \ increases with relative 
density. 


— 


increments when yielding takes place. When the stress increments. are > purely a 
_ @lastic and fall within the yield surface, it is assumed that no changes in the ne 
- effective stresses occur other than those caused by the dissipation of pore 


Initial Liquefaction and Liquefaction Criterion.—Under cyclic shear stresses, 


as the | pore pressure increases and the effective stress decreases, s, the stress 
_ fatio t/a’ increases. If this process continues, the element of soil in an undrained © 
condition will approach the failure state of t/a’ = tan . Under monotonically-in- | 
creasing stress, such an element of soil s considered to have failed. Under : 
Stresses, the material approaches a “‘near failure’’ condition but complete 


i! 
| 
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failure does not occur. A symptom of this condition is a steady increase in 
strain amplitude. This state of near failure is referred to as initial liquefaction. 

_ This definition of initial liquefaction was introduced by the writers in Ref. 6, — 
and it must be pointed out that it differs from the conventional definition of ve 
initial liquefaction. Experimental evidence (9) indicates that under cyclic stresses 
the initial liquefaction occurs at stress ratios slightly less than failure values 

4 under static loading. Thus, the initial liquefaction criterion under ene biaxial | 


7. > which a a=a number slightly less than one (a = = 0, 95-0.98). It is important — 


to note that initial liquefaction is essentially a failure criterion and does not 
directly depend on the pore pressure. The excess pore pressure at — 
liquefaction is less than the initial effective stress, and its value depends on : 
the history of shear stress prior to initial liquefaction. At initial liquefaction 
: effective stress has not vanished and the sand cannot be considered ed liquefied, 7 
but under additional cycles it is well on its way toward liquefaction. = 
| The response of sand to cyclic stresses after initial liquefaction is ‘different 
than that response prior to initial liquefaction. This change of material behavior _ 
_ is evidenced by an increase in the strain amplitude and a change i in the effective 
Stress path in any typical constant stress amplitude test | (9). yt post-initial 
_ liquefaction model is required to represent the behavior after the occurrence _ 
_ of initial liquefaction. However, the mechanism of the behavior of sand after — 
initial liquefaction is not clearly understood at present. Of necessity, the simple — 
_ stress path model proposed in Ref. 6 has been used in this study. It must — 
i Ss be emphasized that the post-initial liquefaction model is not jon usually 
4 very important aspect of the method of analysis. Complete liquefaction usually 
- follows the initial liquefaction in most loose sands. The few seconds difference — 
between the occurrence of initial liquefaction and the onset of the complete 
" liquefaction is a minor detail. However, it is important to recognize the existence 
of the state of initial liquefaction as a critical intermediate point on the path 
_tocomplete liquefaction, 
As it is commonly used, the word liquefaction refers to a State of zero or 
residual effective stress in a region within the sand mass. Obviously, in a liquefied — 
region of sand the pore pressure is equal to or slightly less than the total mean _ 
stress. The liquefied soil has very little shear resistance, and is capable of — 
q undergoing flow and large deformations. The criterion for the onset of liquefaction 
can therefore be stated as the effective stresses reaching a specified residual 
_ value (usually a small percentage of the initial ef: “fective stress). Upon liquefaction, 
= the modulus of sand in the liquefied region is ; reduced to a fraction of its 
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In reality, the level grounds are simultaneously subjected to three | components | : 
of shaking during earthquakes. Most studies of seismic response and liquefaction — 
of scenes ground in the past have considered only « one horizontal 


| 
| 
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; idealization of the problem are investigated in the case studies described in 
4 the following sections. . The case studies are specifically directed at providing —- 
_ information on the influence of the vertical component of ground acceleration 
and on the influence of the interaction of the two horizontal components of — 
_ A hypothetical soil profile was used for all the case studies. As shown in 4 
_ Fig. 4, all the material properties are uniform over the 100-ft depth of S 
"soil profile, the shear modulus which varies linearly with depth. 


« in n Fig. 4 was subjected to the vertical component and ¢ one s horizontal component i 

_ (SOOE) of the El Centro earthquake record of 1940. The same system —_ 
also analyzed for one horizontal component of the earthquake. The pore-pressure y 
time histories for two points at the depths of 17.5 ft and 22.5 ft below the _ 
_ ground surface are given in Fig. 5. The horizontal component of the earthquake 
alone causes liquefaction extending to a depth of 25 ft. However, when the 
vertical components are included, the liquefied zone extends to only 20 ft. | 
It can also be seen from Fig. 5 that at the depth of 17.5 ft below the ground 

_ surface for the horizontal component alone, the liquefaction occurs at about — 

5 sec, but for horizontal and vertical components the liquefaction occurs at 

about 7 sec. It appears the that t the inclusion of the vertical component slightly 

reduces the zone and increases the time to liquefaction. The major effect of | 


5 the vertical aye of the base apes is that it causes oscillations 


— _ The three components of the El Centro Earthquake of 1940 were used as the : 
acceleration inthe case studies. 

oF Vertica. Component or Base AcceLeraTtion 
To study the influence of the vertical component of base acceleration on 


a 
in the pore pressure. However, the effective remain unaffected from 
4 In the analysis shown i in Fig. 5 the depth of the water table was 2 ft below 
the ground surface. The same analysis was repeated for the water table at ‘ 
5 ft below the ground surface. The pore-pressure time histories for this case ; 4 
_ are shown in Fig. 6. In this analysis, no liquefaction occurred. The same onal 
: as in the previous analysis can also be observed in the absence of liquefaction. 
It can be seen from Fig. 6 that at 10 sec the pore pressures are slightly higher 
fo or the horizontal component alone, than when the vertical component is included. 
_ Overall, the analysis results presented in this section seem to indicate that — 
the influence of the vertical component of the earthquake on the pore pressures 
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FIG. 5.— Influence of Vertical Component of Base Acceleration on Liquefaction 
is minor, aside from causing some oscillation in the pore pressures. The 
developmen of liquefaction and the potential for liquefaction are not ot significantly 
affected by the vertical component of the base acceleration. if one zo | 
_ The observed lack of significant influence of the vertical component of the 
earthquake on the liquefaction can easily be explained. The vertical component 
of ground shaking produces changes in the total vertical stress which in tum 
- Cause instantaneous changes in the pore pressure. However, due to the rapidly- a } 
oscillating nature of the vertical motion, sufficient time for the dissipation of 
the oscillating pore pressures is not available - and thus, the effective a 


and the aprnaton potential remain almost unaffected by the vertical component 
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FIG. 6 —Influence of Vertical Component of Base Motion on Pore Pressure inereaee 7 , 


of ground shaking. The time histories of « excess pore pressure ratios sé 
in Figs. 5 and 6 clearly support this reasoning. = 

"INTERACTION © OF Two ‘Horizontat Components oF Base ACCELERATION 


In this section, the results of some analyses of horizontally-layered ground 
- subjected to both components of horizontal acceleration will be presented. The | 
main purpose of the analyses presented in this section is to compare the response is 


_ of the horizontally-layered ground under the simultaneous action of the two 


_ horizontal components of base acceleration with the response of the same system r _ 


_ it is subjected only to one component of horizontal base acceleration. 
‘The hypothetical s soil profile analyzed is shown in ‘Fig. 4. The water table 
& 5 ft below the ground surface. This soil profile was analyzed with the following _ 


base accelerations: (1) SOOE component of El Centro earthquake record of 2 o 


_ 1940 with peak acceleration of 0.314 g; (2) S9¥OW component of the same 


earthquake with peak acceleration of 0.214 g; and (3) both aforementioned, 


>. components applied simultaneously in a two-directional analysis. Some selected if 

= of these analyses are shown in Figs. 7-10. Stress-strain diagrams at 

wo depths are shown in Figs. 7 and 8. me extent of interaction ya 


| t 

| 


the two horizontal components of shaking is evident by comparing the stress-strain | 
_ diagrams for the one-directional analysis of the SOOE component or S9OW — * 


component with those of the two-directional analysis of both yay — 


- analysis at several depths is given in Fig. 9. It is interesting to note that at 
greater depths a reasonably well-defined main direction is present. However, _ 
near the ground surface, no obvious main direction for the shear stress resultant — 
exits, and the shear stress _components in 1 two directions are Of the same 
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FIG. 7. —Comparison of Stress-Strain Diagrams from One-Directional and Two-Direc- 7 
tional Analyses at Depth of 32.5 ft Below Surface 82) 


2 Pore-pressure time histories for the three are Shown in a Fig. 


ot 


sn Fig. 11. Also shown in Fig. ll are the time histories of liquefaction for ff 
“the nestor with the SOOE component of the El Centro 


an 
 & 
= 
— 
| 
1.2 and 1.3 to peak accelerations of 1.2 g x 0.314 g = 0.377 g and 13g | 
os x 0.314 g = 0.408 g. It is sometimes argued tiat the main effects of two-directional _ 
a / shaking can be reproduced in a one-directional analysis if the base acceleration _ 


of the resultant of the peak accelerations from two directions. In this case, 
for the El Centro earthquake, the resultant peak acceleration is [(0. 314 g)? 


‘occurs much earlier and is : much more extensive than those of the one-directional 

. analysis with the increased peak base accelerations. This is an indication that — 

- some important features of the interaction between the two horizontal components 
shaking, which increase the liquefaction potential, 


tional Aneiyess at Depth of 47. 5 ft Below Ground Sudtece 
one-directional even with the increased base acceleration 


as analyzed earlier, may indicate no liquefaction, whereas a two- ae ae 
analysis of the same system may indicate occurrence of liquefaction. 
Additional information on the differences between the one-directional and 
two-directional analyses is evident in Fig. 12, which shows the distribution of i 
the pore pressure with depth at 4.5 sec, prior to occurrence of pre meer 
in the two-directional analysis. It can be seen in Fig. 12 that the pore pressure g 
from the water table to a depth of about 35 ft below the ground surface is 
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| 


the ground surface. “Below | 55 ft from the ground surface, the pore oe 
te almost the same for two-directional and one-directional analysis. When the 
z of the base acceleration in the one-directional SOOE component is 


increased by a factor of 1.2, it can be seen from rd = that the pore pressures ‘ 


9.—Traces of of Shear Stress Resultant from Two- ‘Directional | Analysis 

r distribution i is different than that of the two-directional analysis; the pore laine : 

te still smaller than the pore pressures in the two-directional analysis to a depth 
of about 30 ft below the ground surface. Since increasing the amplitude of : 

the base acceleration does not bring the | pore pressure ‘distribution any closer 

to that of the two-directional analysis, it can be expected that the ——e 
7 patterns of the one-directional analysis with increased base acceleration and 
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i FIG. 11 —Time Histories of | Liquefaction in One- and Two- Directional Analyses: a 
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FIG. 10.—Comparison o and Two- 
Directional Analyses _ “wee 
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described earlier, bui the peak base accelerations were reduced to prevent the _ 
; occurrence of liquefaction. The peak base accelerations used in these analyses a 
0. 150 g in the SOOE direction and 0. 102 2 g in the S9OW direction. 


‘spectra. The pow "surface response spectra were computed for the “cases 
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FIG. 12.—Excess Pore Pressure Ratio Prior to of in One- and 
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ound surface response spectra computed one-directional along 
_ with those of the two-directional analysis are shown in Fig. 13. The two-directional 

- analysis gives surface response spectra which are generally similar to those 
resulting from one-directional analysis. However, some significant differences 
do exist. The two-directional analysis yields higher peak spectral accelerations - 
in both directions. Moreover, the position of the peak spectral acceleration 

4 ie 


™ jj. The extent of the response interaction between the two components of the a . 
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FIG. 1 3.—Comparison of Surface Response Spectra from One- and Two- Directional - 


produces some changes in the surface response spectra, as as compared to a 
without affecting their overall shape 


a in the SOOE direction is shifted to a lower frequency ‘in the eaiteeaient — 
= y, the interaction between the two components of shaking _ 
wivadsd od) sew lebax wea A 
to 
sigan 


A method of analysis was presented for evaluation of seismic response and 
liquefaction of horizontally-layered ground subjected to multidirectional s shaking. 
15 _ The saturated sand below the water table was modeled as a two-phase medium, 4 
= the porous granular solid and the pore water as the constituent materials. 
q new material model was introduced for the behavior of sand under cyclic 
= shear stresses. A new computer program, LASS-III, was developed and 


on the development and the pattern of liquefaction. hs ve 


2. Interaction between the two horizontal components of shaking significantly ‘ 
influences the and the pattern of liquefaction. The potential 


horizontal components 


analysis, such that the peak base acceleration equals the ‘resultant of the peak 


~~ Increasing the amplitude of the base acceleration in a one- ‘ore 


base accelerations in two directions, does not completely teproduce the om) 
_ effects present ir in a two- directional analysis. — = 
“+ Interaction between the two horizontal components. of shaking a 
surface response are different than those of one- 
tional 


This study i part. a research project sponsored by National Science 


Foundation under Grant Number ENV- 76-00626. ip 
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_ The elastoplastic material model described in this : section is a special <4 


‘for representing the behavior of sand under the specific stress condition which 
occurs in the analysis of horizontally-layered ground under multidirectional 
shaking. The sand i is assumed to be under an initial effective ‘Stress Saal 


(16 
in which M = tan @ = the failure parameter; and 4 = the internal angle of 


- friction. This failure criterion implies the shape of the failure surface in the 
stress space is a circular cone, as shown in Fig. 14(a). Moreover, it implies — 
that the shear strength is independent of the direction of the resultant shear — 


stress. Thus, under monotonically-increasing shear stresses the material is a 


— in the _ horizontal plane, which is a ens assumption for sand 


| 
| 
if “aa 
allure Sur ace.— he Mohr-Coulomb failure criterion can be written in terms 
the resultant shear stress rt = + as follows 
| 
: 


Yield ates —The yield ‘surface is s also assumed t to have € a cor conical shape, 
but the axis of the yield surface does not coincide with the o’-axis. Moreover, 

P the axis of the yield surface and its radius change during the process of deformation 

and the yield surface undergoes a combination of isotropic and kh kinematic — - 
hardening. The cross section of the yield surface on a plane of o’ = = constant — “ 
is an ellipse, as shown in Fig. 14(5). The equation for the yield ‘surface is — 


the 


Failure Surface 


peo 25 shows i 


x \ 


the: kinematic and the isotropic A stress _vector, 
can be defined as the difference between the s stress vectors, T aS shown 


The 1 for the. elliptic is written in terms of stress 


ATA 


4 = the hali 
& minor axis of the ellipse, as shown in Fig. 15. T 
Bas. 23 and 24 Fg, 25 
| 


and a = the angle between the axes The 
superposed T denotes the transpose of a matrix. 
to the stress increments by using the associated flow rule, which specifies that 7. 


FIG. 15 Yield Surtace in Plane ae 


‘the vector of the incremental plastic shear strains is sabi to the yield surface 


| 


tt can be seen ‘from Eq. 20 that the function. H = the ratio of dt p projected 


tothe yield surface and dy” 


in which \d-y’ "d = the resultant plastic str strain increment. 
it is appropriate that H be called the “plastic modulus” 
The isotropic and the } kinematic of the is by the 


| 
| 
ay 
| 


MULTIDIRECTIONAL SHAKING 

in which g and c = the isotropic and the kinematic hardening ikatend, 
Tespectively. | Eq. 24 indicates that the change of the center | of the : 

‘surface, dr., is along the normal of the yield surface, allies sr. 


ts material behavior is defined through the three cio H, ¢ 


a _ which are related through the consistency condition, df = 0 
af af aa’ 

e+ 


=O 


substituting Eqs. 23 and 24 into Eq. 25 and using my 20, the following aaiteid 
between three material functions cont a be obtained 


Tepe 
in which B = bl(af/ar,)? + (af/ar, ) is necessary to 


two of the preceding three functions. . The third function re be evaluated from a 
‘The plastic modulus i is determined from the unidirectional stress-strain 1 relation. Se 
“Here it is assumed that the unidirectional relationship between the stress ratio, 
fo’ o’, and plastic strain is given 5: a hyperbolic function, as shown in Fig. 


We 


H ‘Ho’ 


= (72 + 12)'/?. It can be seen that the initial value of H a 
= 0)isH and the failure is reached (+ = Mo’ ), H becomes zero. 
: (Of the remaining two functions, the kinematic hardening function, Cc, is 
_ specified, and the isotropic hardening function, g, is determined from Eq. 26. - 
The equation for d+. can be written in the following form by substituting Eq. 


It can be seen from this equation that c/H = the ratio between the ength = 


vil of dt,, and the projection of d+ along the normal to the yield surface, n Td. 
It has been experimentally observed (9,15) in unidirectional tests that a 


yielding occurs in one direction, the yield stress in the opposite direction remains 
_ unaffected. As seen in Fig. 16, satisf ying such a 4 oats results in the following 


| 
a 
ke 
| 
y 


function g can now be determined from Eq. 26 


Both: are to plastic which is the 
slope of the stress ratio-plastic strain curve. Here, a hyperbolic stress- -strain 
relation has been used. However, this is not an essential part of the eae - 
and other similar ‘Stress- strain relations be used. 


OF 


wal 


of Kinematic Hardening 


es 
al The wudegie shear stress-shear strain ‘relations can be derived by using 
standard methods in plasticity (16). For the nen model th the incremental 


dt, dy, 


G= the elastic p (or unloading modulus); and n, and 


he mageneEe of the unit outwardly normal to the yield surface. _ ee 
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4 The purpose of this paper is to investigate the problem of analysis and design a 
of anchored sheet pile bulkheads anchored with horizontal and sloping tie rods, 
with particular emphasis on the influence of cohesion on sheet pile behavior. 
rs Basic princi i 


principles of soil and structural statics are used to analyze stability conditions 
for the bulkheads. In developing a method of analysis of sloping anchor bulkheads, é 
a critical review of the existing methods of analysis of bulkheads anchored > 
with horizontal tie rods i is presented. Anew dimensionless parameter is is developed 


To 
Previous Work on AncHoreD 


_ Chronologically, the first salamat after World War II to investigate experimen- 
tally the flexible wall interaction with soil i is found ir in Refs. 3, 4, and 5. Regardless" 6 
of the small scale of the ex experiments, the basic facts wd the behavior of flexible — 


the war, the interest in anchored flexible bulkheads arose in the United States 
Navy Bureau of Y ards and Docks; a research program was initiated in o 
at Princeton for the study of active earth ‘Pressure | on - flexible plates hinged bed 
on the bottom. A study of both active and passive earth pressures on a flexible — 
plate representing an anchored bulkhead began in Princeton in 1947-1948; 20. 
‘tests were performed and reported by Tschebotarioff (13). The Lego og 
Princeton tests were limited to highly flexible plates of maximum height, 
= ; = = 249 cm (98 in. ds the Rowe- St. Andrews University tests (10) covered i a wide 
om (36 i in. ); and t the Browzin- Paris tests s included also ¢ a a wide range of flexibility 
4 levels with maximum plate height, L = 36.0cm(14.2in.). me 


"Research Mgr., Structural Projects, Office of Nuclear Regulat 
States Nuclear Regulatory Commission, Washington, D.C. 20555. | 
eee. —Discussion open until October 1, 


ory Research, United 
1981. To extend the closing date one month, 4 
_ a written request must be filed with the Manager of Technicai and Professional Publications, | 
_ ASCE. Manuscript was submitted for review for possible publication on August 17, 1978. 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
of the American Society of Civil Engineers, ©ASCE, ‘Vol. 107, No. GTS, May, 1981. 
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| confirmed by larger and more extensive experimentation by Rowe (10). During 
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comparison of the tests performed with flexible plates in laboratories to 
_ the behavior of actually-built anchored bulkheads of steel sheet piles or concrete _ 
: or timber sheet piles, requires a nondimensional number which takes into account — 
both soil and sheet material characteristics. The review of recommendations _ 

- which resulted from the works of Princeton and St. Andrews University, = 
well as from the writer’s research will be presented after the development of 


IDENTIFICATION OF FreE- Eantu Support Conpitions 


' is well-known that the ene modeling of soils is generally v very difficult. ie _— 


Wi wi 


Ap 


fe a to ou sul: alevel ui 
FIG. 1.—Loads Acting on Bulkheads ond for Unit of Width 


plates. In this case, a true similitude may not be established. _> > 

If one wishes to compare the results of tests on plates to the actual ‘behavior 

of anchored bulkheads, she following distinction must be used; for sufficiently 
kg small values of 2D, the buried part of the system w will move in such | a way pig 


; : elastic properties of bulkheads dif fer er - substantially fr rom phead tested laboratory ory 


that "Ee and active soil pressure will be developed, as as shown theoretically mea 
in Fig. a condition called free-earth support, and such a system 


flexible bulkheads. =f 
a \ Tee 
— _ 


be: asa ‘simple beam resting on two supports. The upper support 
is located at anchor rod, at z = 0. For cohesionless soil, the lower support 
e is located at the centroid of the triangle representing the passive soil resistance _ 
diagram, i.e., at z = pH 7 2D/3. For cohesive soils, the lower ‘support will 
be at the centroid of the passive pressure diagram which includes the rectangular 
block representing the action of cohesion. It implies that the centroid = be — a 
located slightly higher than in the case of cohesionless soil. 
7: a As soon as the depth of penetration, D, becomes sufficiently large to create 
restraint of free deformation, the system becomes statically-indeterminate. 
_ The free-earth support conditions cease to exist. Then, in addition to equations ns 
statics, equations of deformation must be introduced. The unknown redundants 
will include: (1) The rotation of the end section, point z = pH + D; and 
® the second unknown will be the translation of the same point toward the - 
Passive sic In equations . of deformation, a horizontal subgrade ‘reaction 


 owonk 

— 


= BULKHEADS 


ist 


must be used, which is and beam elastic properties. Conse- 
quently, for statically-indeterminate systems, similarity between the tested plates © 
and the bulkheads does not exist. On the grounds of available experiments 
by the writer (3) and Rowe (10), conditions delineating the statically-determinate 
_ laboratory tested plates from the statically-indeterminate tested plates might — 
be analytically established. Then, the actual bulkhead characteristics might be ha 
mathematically investigated in order to determine whether or not they may e: 
be regarded as statically-determined systems. 
7 _ With reference to the nomenclature of Fig. 2, the. intensity y of the load, w, . q 
cK 
- intensity of the reaction, r (in per to resist 


| 

i 


+2cK 


in whieh; the vertical effective stress, &, or G,, = a function of depth iow ” 
the soil surface 442 oni 


_ Notations H and wH are shown in Fig. 1; p parameters y = weight of the soil 
(in pounds per cubic inch); B = width perpendicular to the plane xy (in inches); , 
and K, = tan’ (45 — /2) and K, = tan’ (45 + /2) = coefficients of active a 
and passive earth pressures, respectively, i in which @ = angle of internal friction — 
_ of the soil; and c = cohesion (in pounds per square inch). = © 
load, w (active soil pressure), and the reaction, r (passive soil resistance), 
are known at any point, z. However, there are two points of discontinuity — 
the load functions, w(z) and r(z), at z = (2, /yKi*) — — ») and 
_z = pH. Consequently, oe for w(z) and r(z) must be established for 
From physical conditions of t the statically- -determinate beam of Figs. | a and 
24 it is seen that the maximum deflection must occur between points ; z= 
and z = pa. The load at each part of the beam might be expressed in terms 
of the maximum load intensity, w,,,,, and the reaction might be the depth 
in terms of the maximum reaction intensity, r,,,, both acting at the depth 


| _of solids is 


W max 


in which E = modetes' of elasticity of the nit of the bulkhead: I = the 
moment of inertia of the bulkhead cross section for the width, B; and S(z) 
"= variable shape factor which | provides expression for the load, w, through 
the maximum load intensity, w,,,,, i.¢., = w,,,, S(z). The shape factor takes 
- into account the location of the supports and the load distribution. The shape 
factor might be expressed by a polynomial 


a load (see Figs. 1 pest 2). 


maximum deflection, x,,,,, Will occur at z = » from the 


condition 


_ Upon Cinieen: the function w or S(z) by z four times (Eq. 5), and 
upen ‘hana z = z,, by Eq. 7, the value of S(z) at z = z,, will be obtained, 
= Si. max» in which = the Shape factor at the of maximum 


| 
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The deflection will have the expecssion 


— x(z). Canigadood between the two might be based on the ratio of the maximum = 
‘deflection, X max» to the total length, L, i.e., based on the relative deflection. = 
a P _ Taking into account that the maximum effective stress at z = »H + D is 
C= = yL, in which L = H + D (see Figs. 1 and 2), the ‘maximum > intensity, i 


between the laboratory and bulkheads expressed t by 


= 


pol 
a for similar statically- -determinate beams of the same geometry 
in 12, will be denoted, p,, is Ya 


(12) 
load distribution is the same, , and K pes the soil property, is the mn. the key 


which is dimensionless for homogeneous units. It may be called the dimensionless 
- plate and a sr perry For cohesionless soils, i.e., when c = 0, the DFN becomes 


which depends only on the geometry of the system (laboratory plate or bulkhead), oA 
and the effective weight of the = 1, 13 ond 14, 
For r pure cohesive (> = using 13 and 14, the 
for the DFN becomes 4 


flexibility number (DFN). It determines s the statical similarity between a laboratory 
{ 


1-— 


| 
| 
| 
| 
| | 
| 
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Because ‘Eq. 15 ‘contains L in ‘the second te term, a laboratory | plate cannot 
be regarded as a small structure, as is the case for c = 0 soils. In general, 
the bulkhead must be modeled in the laboratory for c 4 0 soils with the following ¥ , 


‘in which M stands for models; and P stands for prototypes. A first step in 


_ satisfying Eq. 16 would be to select the model geometry BL*/EI in such a — 
way as to make = = this results in the following equation 


. Assuming ‘the unit it weight, to be the same for both model wien prototype, 
the law of on of elastic behavior, willbe 


For fo ora 2 modeling ra ratio 1: 5, and a stiff’ clay of the prototype; 
c = I ton/sq ft, the modeling clay must be “‘very soft,’’ with cohesion, c ot 
3 0.2 ton/sq ft. Tests with cohesive soils are not available. In order to use - 
- available tests with cohesionless soils, characterized by the DFN, Po, the influence 4 
of the term 2c/yL willbe investigated. 
Available tests which determine the border line between the statically-deter- 
minate bulkheads (with free-rotating bottom support), and statically-indeterminate _ 
bulkheads (with bottom support restrained against rotation and 
the tests by Rowe (10). Rowe, with similar reasoning to that of the writer — 
the DEN, has determined a aratio = 


ome 


__ However, Rowe’ $s ratio has dimensions and is not applicable to cohesive soils. 
‘He termed p the flexibility , number (FN), in which the length, L (his notation 
‘is L = BA), is in feet, for modulus of elasticity-pounds per square inch and 
_ moment of inertia-pounds per inch‘ for 1 ft of width of the plate. Both the 
DFN, p,, and Rowe’s FN, p, might be used to compare laboratory flexible 
plates representing scaled-down bulkheads with actual bulkheads as long as" 
the bi buried portion, L is deform. These ratios are no longer valid 
, because the lower portion of D ome 
restrained, i.e., the system Statically-indeterminate. bas 
& The experiments of Rowe (10) provided a clear understanding of the behavior 
_of laboratory plates; however the flexibility numbers, p or p,, are not applicable ‘ 
Z comparisons of f laboratory plates to actual bulkheads restrained at the bottom, >, 
~ e., , having fixed-earth support. Following Rowe, both Terzaghi (11) Tschebo- _ 

_ tarioff (14) used the flexibility number for fixed-earth support conditions, and 
& P recommended Rowe’s reduction curves for bending moments based on p values. o 
Since p and p, are not applicable for fixed-earth pret conditions « deducted 
from tested plates, reduction curves cannot be used. 2 tits — 

Rowe’s FN, p, converts to the writer’s DFN as follov 


i 
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’ Following Rowe (oe Ref. ef. 10, pp. pp. 4l and 45, and Fig. 14) values of his 
‘FN, p, which delineate free- fixed-end support, i.e., statically-determinate 
from statically- systems are for loose sand, log. p = -3. 5, and 
dense sand, log p = —4. 125 (by interpolation from Ref. 10, Pp. 45, 
14). These correspond ‘to the rounded- -up values of the DFN, p, = 5.0 for — 
loose sand, y = 90 lb/cu ft (1.442 x kg/cm’ and approximately to 
= 1.0 for a dense sand of approx y = 110 Ib/cu ft (1.762 x 10~ * kg/cm* ), 
in which y, the effective unit weight = a measure of the degree of density. y 
The corresponding DFN for cohesive soils, - following Eq. 15, will be ‘smaller 
according to the cohesive term 2e/yL. 
: Sine the DFN depends only on the geometry of the s system, and the density 
- and the cohesion of the soil, the limiting length of each type of sheet pile 
beyou id which the the system ceases to be e statically determinate : can be calculated. 
Exe 


Cohesion 
inkg/cm? 


4 


Jo 

, _ Logarithm of moment of inertia pa unit of width, log (1/B) 


1/B in cm3, 1 in? = 0.0610 47 


FIG. 3.—Length of Sheet Files Limiting adh Determinate huis in ue to 


Physically, the presence of cohesion Teduces” the actuai pressure, which 7 } 
equivalent to stiffening the sheet piles. Of 4 
The limiting length of typical sheet piles will be calculated using limiting 
- values of the DFN for loose and dense noncohesive soil, c = 0, for soft clay, — 
c= 05 ton/sq ft and for medium clay, ¢ = 1.0 ton/sq ft. Eq. 13 will be 
converted for calculating the length, L, for, cohesive soils (K, as 
L which W = EI/yB; me p, = the limiting value delineating the statically- ~deter- 
: min inate from the Statically-indeterminate system. For cohesionless soils: (c = 


L 
‘ | 
q 


Calculation of the 


L, of Eq. 21, of the 


presented in Fig. 3. In Fig. 3, the limiting length, L, is shown as a function | re 


of the logarithm of pile stiffness 


per unit of width, for loose and dense soil me 


combined with cohesion characteristics, c = 0 ton /ft, 0.5 ton/ft, and 1.0 ton/ft. 
In denser soils, the limiting length is smaller, i. e., the system becomes statically- 
indeterminate with a shorter length than in looser soil. Cohesion, as expected, 
_ provides the possibility of using longer sheet piles prior to the system becoming 


‘Wes 


xed TABLE 1 —Limiting Sheet Pile Length Calculated by Limiting DFN posit... 


modulus| ment 
of wall, 
in cubic | inertia 
inches | 1/B, 
per foot | of wall, 
Manufac-| (cubic | in -ubic 
turer- | centi- inches 
meters | (cubic 


Manutac- defini- | per | 


turer | tion | meter) | meters) 


States 
‘Steel 
United 

States 


gitter 
a 


LIMITING LENGTH FOR STATICALLY- 


DETERMINATE SYSTEM 
(FREE-EARTH SUPPORT DESIGN) | 
FIG. 3, in feet (meters) 


of ratio Soil Characteristics, Cohesion, 
1/B, in in tons per square foot (kilo- 
Subic grams per square centimeter) 


meters 


“indeterminate. For intermediate values of density and cohesion, a ie 


extrapolation between the (CUrves 


density, 100 Ib/cu ft AU 602 x 10 


an be used; otherwise Eas. 21 and 22 might 


be used for the DFN between values 1 and 5 of Rowe’s tests (10). For a 


~ kg/cm’), DNF p = 1, and 90 Ib/cu ft. 


(1.442 x 10~° kg/cm’), DNF p = 5. In this analysis, DFN p = | was used 
with soil density 110 Ib/cu ft (1.762 x kg/cm®). 


__ Table 1 provides characteristics from the 1979 catalogs of the two major — 
manufacturers of sheet piling: United States Steel, and Stahlwerke Peine- -Salzgitter — 
AG, in limiting length is listed for particular sheet 


: 
| le | (17,350) | (7.8 | 
| 
j 


a wide range nge of the ratio i.e., e., from 15.4 cu in.-476 cu in. (for the 


metric conversion, see the table). For a small J /B ratio, i.e., 15.4 cu in., the i 
in Table is 26 ft for dense cohesionless soil. For a large 
: ‘tia, i e., 476 cu in., for dense cohesionless soil, the limiting length is 56 
. ft. In cohesive soils and in soils of lesser density, the limiting length is even 
ao eo In most practical cases, bulkhead lengths are less than limiting lengths, — 
a , the bulkheads are usually statically-determinate structures, satisfying free- _ 
ie the stability of the majority of bulkheads with | horizontal anchors 
‘g be correctly analyzed by simple statical principles, it is assumed that the ” 
same principles apply to the analysis of bulkheads with sloping anchors for : 
at 
The: in the of wich Ko and anchors 
-~ onsists in the fact that the sloping anchor tie rod causes a vertical sos pe eel 
acting — LASS + DIE — 
a  Bulkheads with sloping anchors in granular soils that offer as substantial friction 
on the sheet piles do not represent difficulties for construction. The ool 
component caused by | the sloping tie rod acting downward on the sheet piles _ 
can easily be transferred to the soil by friction, and the sheet piles — 
by the general principles of soil statics. A numerical example of this analysis = ; 
has been provided by the writer (7). In cohesive soils without internal friction, 
= 0, or with ‘small friction angles, analysis of bulkheads with sloping tie 
_ Three conditions of stability must be in onder to provide a sufficient 
depth of penetration for adequate stability of bulkheads. The first is the 
overturning of the structure around point 0 (Fig. 1). The second condition : 
the stability against the sinking of the sheet piles” downward by the action of — Bo 
the vertical component of the reaction caused by the sloping anchor tie; a 
the third condition i is the soil stability against heave at the bottom of the excavation. ‘ 
_ These three conditions are analyzed in this report based on the principles’ 4 
of fundamental soil properties, and accepted methods in theoretical soil mechan- 
ics. Friction angles of ¢ = 0°-10° were considered in order to demonstrate a 
the major effect of the friction angle on the shear strength of the essentially- -cohe- 
sive soil. The shear strength composed of cohesive properties of the soil at 
+ od = 0, and the internal friction  # 0 are responsible for the stability of the 
bulkheads. The division of the overall strength into the components’ cohesion 
and friction is necessary in the principles of is 
Tt was found in this rep report that all three conditions of stability might be | { 
expressed in terms of the dimensionless stability number, 4c/-yH, as that used / 
by Terzaghi for evaluating the stability of the soil in braced cuts against heaving 
of the Ref. 12, P. 193, formula 57, 


q 


will read as follows: _ 


d 


_ The bulkhead must be 2% against overturning around point 0 (see 
oe 1). The load on the bulkhead is caused by active earth pressure, pase, * 
on Fig. 1, by a surcharge and by hydrostatic differential pressure in harbors. The 


earth pressure is reduced by the action 1 of cohesion. The stabilizing reactive | 


forces are provided by passive pressure and increased by the action of cohesion. 
_ A safety factor, G,, might be introduced assuming the passive pressure ordinate 
_ ‘The cohesion might also be reduced to the value c/G, on the passive side. 
‘The equation « of equilibrium will be elaborated in order to demonstrate the 
influence of contributing factors on the stability, 
_ The equilibrium against overturning requires that the depth of penetration, — 
—D, be sufficient in order that the sum of moments about co 0 be = ey 
=a The equation of equilibrium, the aie of moments ‘men point 01 in 


K,L E L- HO + 


+——6,K,D\—D+nH 


in effective p = = (H+ D)y; & = Dy, where ¥ = effective 
weight of soil; and G, = safety factor, as defined by Terzaghi (12); the ae b 
factor, G, should have the level from 1-1.5, as analyzed by the writer (5). 
‘The last term of Eq. 24 with the plus sign, compensates for the ‘ ‘tension”” 
_ caused by cohesion (which actually does not produce | tension c on the bulkhead), \! 
shown by a triangle in the upper part of "Fig. The moment ‘caused by this 
triangle load is small and acting in the clockwise point 0, ie., 
Neglecting this term, Eq. 24, in compact mm -* resolved for the stability 


nM which the geometric characteristics of the bulkheads are: B = > L/D; en. 
+B — 1) + 1, where L, D, H, and pH are shown i in Fig. 1; andA = Kp/G, i 
Active lateral earth pressure coefficients, ¢ ,and passive lateral earth pressure 
coefficients, K,, may be selected following the Coulomb theory or or e.g., principles 
involved in the treatise of Caquot and Kerisel (8), or others. The influence 
of friction between the wall and the om may also be taken into account oe 
using the proper of K and as in any 


| 
4 
4 
fi 
7 
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/ 


/ 


anchor pull | 


overturning} 


FIG. 4.—Stability Number 4c/yH Versus Parameter 8 = 1 + H/D for Three Condi- 
tions of Equilibrium: Sinking, and Soil Heaving 
4 The 1¢ depth of + penetration, D, must be | selected to satisfy ‘the equation of — 
4 equilibrium (Eq. 24 or 25) for given H. At c = 0, Eq. 25 reduces ‘to the ‘equation 
>+3p7 a — AGa 
A for which a eeitie solution i is published by the writer (5). A method of application — 
_ to stratified soils is also available pone in which coefficients K, and K, might 


ol of the 
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for ‘ pressure diagrams, and Eqs. 24-26 used. 
The graphical solution shows a high sensitivity of | the required depth, D, to 
Ls the assumed angle of internal friction, o. The same fact is also pointed out 
by Kovacs, et al. (9). For illustration of the variation of the depth, D, the 
x stability | factor, 4c/yH, was calculated from Eq. 25 with parameters p = 4 
G, = 1, and angles of internal friction, 6 = 0°-10°. It is seen on the graphs 
in Fig. 4 that at changes i in the angle, , influence considerably the the stability f factor ; 
Acainst Sinxine of BULKHEAD By ACTION oF VERTICAL Component Causeo 


* however certain ‘Teasonings differ from theirs in analyzing the —* agains 
against sinking is adequate if the forces of ‘might 
‘sufficiently developed on faces of the sheet piles to overcome vertical load, 
_A,, caused by the vertical component of the reaction, A,, in the tie rod. The — 
_ Broms and Stille references use the term ‘‘adhesion,’’ and assume that a fraction — 
of the ‘“‘average undrained shear strength’’ may be as adhesion. Since 
the soil statics analysis requires separation of the shear strength | into components — 
7 determined by ¢ and ¢ according to Mohr’s strength theory, a fraction of the 
cohesion is regarded as adhesion in this analysis. The legitimacy of accepting 
_ adhesion in a theoretically pure cohesive soil (¢ = 0) is not examined in this — 
Equilibrium: of forces in the z direction might be written as follows: © : 


in which a = the reaction in the x ¢ direction at point 0; and =Q represents n 


the sum of friction and adhesive forces acting on the faces of the bulkhead — 
_ Sheet piles. Eq. 27 for vertical equilibrium might be written in the developed 
n fo the calculation of the avoid as 


— }(tan a — tan 15) 
K,D+ 
in which G, = same factor of as that used in the equation of stability 
against overturning; see Eq. 24; 5 = the angle of friction of soil on the material 
of the sheet piles, which may be used as 6 S o; and m and n= fractions 
expressing the adhesion in terms of cohesion. . The angle, a, is shown in Fig. 
1. ‘‘Relative adhesion,’’ as used by Broms and Stille, involved fractions B, _ 
and 8, of the unconfined compression strength without separating cohesion — 7 


. 
| 
q 
} 
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from internal friction. Since it is separated following Mohr’s strength theory 
in this his report, . friction on (terms in including factor tan 5) and adhesion — 

factors nc and mc) appear separately in Eq. 28. ;§ 

ety. Eq. 28, in compact form, resolved for the dimensionless | stability 1 factor, reads 


ari 
+ nK}?; B= BK,” (tan 
— tan 8) + Ki? (tana + tan 1/2 (m + nB); and C = B?K, 
(ten a — tan 8) - "Kp (tan a + tan 8)/G 
The ‘‘adhesion’’ and friction on the a active side are relatively small and can bs 
2 neglected for cotiay, | i.e., , to use n =0 and 8 =0. The depth, , D, for given 
_ H, must satisfy Eq. 29. For. ‘delineating ‘conditions which cause the reaction. 2 
in the anchor A, = 0, the sum of horizontal forces in Eq. 28 is equated to 
0, i.e., the five ‘terms in parentheses at 
Parameters of | of Eq. 29 thenreducedto 
Using Eq. 29, the: stability factor, 4c/yH, required for stability | against ‘sinking 
of bulkheads with sloping anchors is represented as a function of the parameter a 


tad (29. For values of 4c/yH seeing: than the the numbers indicated by lin lines foi for 
A, = 0, the anchor i is not required. Yo aur clans od rad? rad? od 
‘Examining the graphs of Fig. 4, it is seen that the stability against overturning - 
_ Tequires larger 4c/~H numbers than the stability against sinking except ect 
y = 0°, where the lines representing overturning and sinking conditions intersect 
; at 8 = 3.7. This leads to the conclusion that the stability “agama overturning 
is the governing condition for bulkheads with sloping anchors. | 
= AGAINST Heavine at Bottom oF Excavation oF OR AT Drevce Line in Harsors sult 
~ Stability against heaving is achieved by the weight of aban overburden soil 
surface below “this level. The method of analysis is based on Prandl’s theory 
for weightless piastic medium resisting penetration of a rectangular infinitely 
7 long rigid body, corrected by Terzaghi (12) for weighted soils. The resulting — 
bearing capacity equations are generally accepted (12). Consequently, this theory 
will be applied here for analyzing the stability against heaving following —— 7 
(12). The load of the soil mass behind the bulkhead, y(H + D) B a 
resistance provided by the soil overburden above the horizontal line at z= 
»H + D at the front of the bulkheads, yDB, must be in the plastic state 
4 of equilibrium. T he resistance is determined by making use of Terzaghi’s bearing _ 
capacity equation. The equilibrium of vertical forces will be written for a unit, " c0 7 


B=1, of width normal to the plane xz of of Fig. 1 as 1 as follows: 


above the bottom of the sheet piles and the friction resistance on the sliding | al 


7 
| 
| 
4 
1 


in which N. N, represents resistance 
due. to the overburden; and N, represents resistance due to the internal soil — 
friction. For small friction angles (0-10), N, = 0. 31 for the stability factor 


For up to to = 10, approximate values N. _and N, were used (12). 
The lines were traced for calculated values of the stability factor, 4c/yH, ion 
the parameter B = 1 + H/D (Fig. 4). It is seen that the curves are a 
_ substantially lower ‘hen the curves representing the stability against ov erturning — 
and sinking. It is well-known that in the absence of seepage, pressures heaving» 
in the bottom of excavations with : sheet pile walls do not occur if the other — 
stability conditions, against overturning and sinking, are assured. Heave ~o 
occur in braced walls within small depths of penetration, but not in the sheet 
pile bulkheads. Even though the investigation was made for essentially clay 
soils (> = 10°), heave of the bottom will not occur i if oe — of er 

_ Broms and Stille (1) report several failures of flexible bulkheads with sloping — 
anchors in Sweden. In an earlier publication, Broms and Stille (2) presented - 
a method for determining the stability number. It appears to the writer that _ 
this stability number is insufficient for evaluating the stability of anchored 
bulkheads with sloping or even with horizontal anchors. In this connection, _ 
it must be stated that the ; analysis of stability against overturning of the bulkheads © 
about point 0 in Fig. 1, is of major importance. Stability against overturning — 7 
takes care of the ‘ “*toe failure’’ of the Broms and Stille terminology. Next in ‘ 


importance is the verification of the stability against sinking downward, i.e., 
against ‘‘penetration failure,”’ following again the Broms and Stille terminology. 
The ‘‘slope failure’’ analysis, i.e., sliding on the ‘‘Swedish circle’? does not 
require any | special consideration in respect to bulkheads with sloping anchors 
if the anchor pull length is sufficient. 
_ The chart in Fig. 4, may show the stability conditions for a bulkhead in 
a clay of such a low strength asc = 0.1 ton/sq ft, = 5° and y = 0.055 
ton/sq ft (1.76 ton/m’). The necessary depth of penetration, D, fora bulkhead 
with H = 12. 0 ft (3.66 m), must be analyzed. The required stability number 
4c/yH = = (4 x 0.1)/@. 055 x 12.0) = 
= 5° for overturning, indicates B = 70 ead, consequently, = HIB 
. 1) will be D = 17.1 ft (5.21 m). For. an anchor of 45°, the line b= 
for sinking, indicates B = 2.04 and, consequently, D = 11. 5 ft (3.52 m), oo 


A dimensionless flexibility number, DFN, p,, which can be to 
similarity between laboratory plates and a prototype for anchored bulkheads, 
i developed. = the free- earth support contiens, the DFN, as it is shown, _ 


— 
| 


"general, properly designed anchored bulkheads in the majority of cases sodden an 
using limiting values for DFN, Py, i.e., p, = 5 for loose soil, and 


sheet pile length i is greater the limiting value, the bottom end of ‘al 
- Bulkheads with sloping anchors in cohesive soils must be designed to the 


ee three conditions Te a single stability factor, 4c/yH: (1) To : 


ig length for the sheet piles is analyzed. Whee 


of the excavation /4H = = Values of N,, N,, and are functions 
of geometric and soil parameters poe the factor of safety. Values of N,, N,, 
and N, are calculated by using the right-hand side of Eqs. 25, 29, and 32, 
_Tespectively. Charts or calculation by trial and error is necessary. = 
: Stability against sinking can be analyzed using reasonable assumptions of © 
adhesion and friction. In granular soils, the sinking conditions is satisfied because 
_ of sufficient available friction. Sinking conditions in cohesive soils are not critical _ 
_ because of low tension in anchor pulls as a consequence of a great depth of nf 
: penetration needed to satisfy the stability condition against overturning. Graphs 
for stability conditions based on Eqs. 25, 29, and 32 for essentially-cohesive 
soils are presented as an example in Fig. 4 for ‘slope angle a = 45°, safety — 
= G,=1, location of the anchor rod at » = 1, and adhesion coefficients — 


A Reema, B., and Stille, H., ‘Failure of Anchored Sheet Pile Walls,” rr of th the 
M4 Geotechnical Engineering Division, ASCE, Vol. 102, No. GT2, Proc. hope 12000 a 
2 Broms, B., and Stille, H., “Stability of Anchored Sheet Pile Walls in | Clay,” Swedish 
oe Geotechnical Institute and Royal Institute of Technology, Stockholm, Sweden. _ . 
a 3. Browzin, B. S., “‘Upon the Deflection and Strength of Anchored Bulkheads,”’ 
Proceedings of the Second International Conference of Soil Mechanics and Foundation 
Engineering, Vol. III, 1948, pp. 302-308. 
Browzin, B. S., “Approximate Computation of Anchored Bulkheads,”’ Proceedings 
7 of the Second International Conference on Soil Mechanics and Foundation Engineering, 
_ 5. Browzin, B. S., “‘Upon the Coefficient of Stability Reserve of Anchored Bulkheads,”’ 
_ Proceedings of the Second International Conference on Soil Mechanics and Foundation 
_Browzin, B. S., “Vertical Flexible | Plate Interaction with Granular | Soils: : Experiments,” : 
7. Browzin, B. S., discussion of “Failure of Anchored Sheet Pile Walls,’ > by Bengt — 
_ Broms and Haken Stille, Journal of the Geotechnical Engineering Division, ASCE, _ 
Vol. 103, No. GT4, Proc. Paper 12831, Apr., 1977, pp. 350-353. 4 altel 
& —Caquot, A., and | Kerisel, J., Traité de Mécanique des Sols, Gauthiere- -Villars, Paris, 
9. Kovacs, W. D., Martin, J. H., ., and Gerig, F. A., ‘“‘Sheet Pile Design | with Variable 
Soil Parameters,” Journal of the Geotechnical Engineering ASCE, Vol. 100, 
No. GT2, Proc. Paper 10310, Feb., 1974, pp. 200-205. 
oo Rowe, P. W., ‘“‘Anchored Sheet- Pile Walls,” Proceeding 


Engineers, Vol. 1, No. 1, Jam., 1952, 


| 
| 
4 


ll. Terzaghi, K., Anchored Bulkheads, Proceedings, ASCE, Vol. 79, Separate ‘Mev, 
ak K., Theoretical Soil Mechanics, John Wiley an , Inc. York 
G. P., “‘Large Scale Earth ‘Pressure Tests with Model | Flexible 
Bulkheads,”’ Princeton University, Princeton, N.J.,1949, 


14. Tschebotarioff, G. P., Foundations, Retaining and Earth Structures, McGraw-Hill 


A,B,C symbols in Eq. 29; 
_= force i in anchor tie rc 
A vertia component of anchor tie rod; 
_ width, transversal to xz plane of Fig. 2; 
factor ad the passive sesiotence or 


height from soil pat ae to dredge line « or depth of excavation; ir 
= moment of inertia; 
L = lengthofsheet piles; 
_M, M bending moments in beams and sheet Piles; 


c= fraction of cohesion on passive side; 
N functions of geometric and soil parameters factor of safety; 
N = bearing capacity coefficients (by Terzaghi); 
fraction of cohesion on active side; 


= intensity of reaction; 

= shape factor for simple beam analysis; 

= shearing strength of of cohesive soils; ee iis oe. 
W = section modules; WY 
an w(x) = varying load intensity on beam; + ‘ 
“coordinates (see Fig. 1); 

p(H/D) +1, = ‘Tatio; 

= effective unit weight (buoyant); 
= angle of friction, soil sheet pile 


s flexibility number with | dimensions in English units 
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ANCHORED BULKHEADS 
= dimensionless flexibility 1 number (DFN); 
limiting value of DFN delineating fr rom 
systems; 
DFN for cohesionless soils, c = 0; 
of vertical friction 
- = effective stress, function of depth z; 
= 
effective stress, function of z on passive side; and 
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_A recurring problem in geotechnical engineering is to estimate the properties a 
f 
a of assumedly homogeneous zones or strata of a soil profile for use in modeling. : 
Aside from the difficulties of measurement bias, the properties of soil themselves 
vary from point to point, and the question remains of how to combine the | i. 
7 mee measurements into point estimates or “‘equivalent’’ parameters for | 
deterministic models. This note proposes a simple procedure for reducing the ai 
5 variance e of estimators of equivalent Peneniers by using a weighted combination 
redundancy of information in closely spaced observations, the differing 


influence of regions of the soil profile due to nonuniformly — loads. —_ 


Gregory B. Baecher,' A. M. ASCE 


Letz observations of the sa same property 
homogeneous zone of the subsurface, and let 2 be the point estimate of wo4 


€ g., finite element analysis). Here z is taken to be a vector of scalar measure- 


7 _ ments, but extension to vector measurements (i. e., more than one property _ 


ig 
at each location) is straightforward. Then, the common estimator used in 
geotechnical practice is either the arithmetic mean _ 


or the mean less some multiple of the sample standard deviation, "a : 


stical point of view, these might be called the ‘ estimators. 
If the expected values of the z, equal the spatial average, —E(2), then Z,, i. 
is an unbiased estimate of E(z), in the sense that its expected value over all 
Ways the n observations might have ‘occurred equals the true spatial mean. = 


TAssoc. Prof. of Civ. Engrg., Massachusetts Inst. of Tech., Cambridge, Mass. 02139. _ ae 
7” Note.—Discussion open until October 1, 1981. To extend the closing date one month, i 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on July 8, 1980. 
‘This paper is part of the Journal of the Geotechnical Engineering Division, esaeia 
_ of the American Society of Civil Engineers, ©ASCE, Vol. hy ah No. GT5, May, 1981. 
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n the case of observations, is also ‘the minimum variance, 
in the sense that its sampling variance across all ways the sample 1 a“ have — 


in which V(z) = the variance of the individual (i.e., point) properties (2). 
_ When the components of z are not independent, but have covariances among __ 
re themselves reflected in a covariance matrix é.; then the sampling variance — 


— 


AV 2; 


in which (1/n) = vector of n element equal 1/n, and is 
the minimum of all estimators. Considerable empirical work has vetted = 


unless the spacings among observations are re large, | tie 
_ of z are indeed correlated . In this context, and depending on the particular 
deposit, “large” means 50 m-100m(1,3).} | 
Considering the class of estimators comprising sums s of observations 


=Swz,.. } 


nom ‘optimal’’ estimator of the spatial mean can be found by adjusting = 
_ weights to minimize the sampling variance of 2 » about the spatial mea 
for the purpose of estimating the spatial mean, ‘but not necessarily 
purposes, 27, can be said to be an optimal estimator. 
; = The variance of Z,, among all ways the observations could have been cen realized, 
‘ “given the covariance matrix C,,is 


“which, an algebraic manipulation, can be shown to equal 


= 


. 


in which w,< 10. 


| 
= 
that in the independent case, the sampling variance of Z, is §|§ 
— 
: eee 
in which fens, and C, is obtained from 
the autocovarie sun property. Minimizing with respect to 
w yields the optimal weights 
(wt) fe, 1]-'fo 
RO? 


mores 
in which = (1, 13 = @, 0); and A = a Lagrange multiplier to 
ensure Eq. 7. Application to the boring grid of Fig. 1 yields the ene shown aa 
and reduces the estimator variance about with 


ee For simplicity this example assumes isotropic two-dimensional autocovariance. 


Extension to three-dimensions requires only that the elements of C, be evaluated 


from vector separations, 
of ¢ course, the error in of the spatial mean are re of ess 


e prediction of engineering performance, i. e. 


= 


volmihety 


WEIGHTS 


‘FIG. 1- —Boring Array for Brent B B Platform, Ret, 1 


deformation, or flow. Therefore, one might expect observations taken at locations — 
especially important to performance to be assigned higher weights than those 
i other locations. For example, to predict settlement under a foundation, an 
observation made beneath the centerline should be given more weight than 
one outside the imposed stress bulb, almost regardless of its covariance with _ 
other observations. For weights calculated from Eq. 101 this would not necessarily 
be the case, as shown in Fig. 2. 
«Let the soil be modelled as a stationary field 


| 
| 
| — 
| 
5 


| 
= = t(x) + +e) 


-zero- mean random tert mn with autocovariance function Cc in in which r 


‘separation. Further, let the performance measure to be mie i depend 


te which g(-) will be said to be the geomechanical model. The question tion becomes, 


how to choose the weights of Eq. 5 ‘such that the | error in estimates of | y is 
That is, in using deterministic models, like most finite element 
analyses, the aoneeny spatially varying properties, fae are replaced by a single 
a 


DISTANCE ,r 
2.—Optimal Weights for Estimate of Spatial Mean : 
the expected squared error between using 2, 


the ‘spatially varying vee 


‘Then, the variance of F to be minimized is the squared error wig 
ny Expanding Eq. 14, and making the assumption that g [2@)] | isa linear function a 
of z(x), then setting the differentials with respect to the w, equal to zero yields - 


| 


NOTES 
G,,) = a vector of sediment 


The effect of the influence factors is to distribute additional to: ¥ . 
8 _ Those observations that are in locations to which the predicted performance, —s 


7 ys is sensitive, as reflected in the derivative of g(-); and (2) those observations eS 


are correlated with the soil properties of nearby locations to which is 
also sensitive, as reflected in the product of the covariance and derivative. oS 


common in site of estimating equivalent 
unif orm parameters for deterministic modeling, simple statistical methods provide | 
a vehicle for improving the way geotechnical analyses are now performed. They 
‘ do so at little additional cost, and without introducing elaborate reliability models. 
a Fay methods are no panacea; they offer a marginal improvement. The major 
er limitations to the accuracy ¢ of geotechnical predictions, of course, ‘still remai hoe 
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Som STABILIZATION WITH FINES 


Sudipta Ss. Bandyopadhyay,’ M. A 


Preheater fines are defined as the material ; produced i in a cement kiln preheater — 

- incident to the manufacture of Portland cement. Because of its ability to control 

- swelling of potentially expansive soils and increase soil strength, preheater fines 

* can be used as a substitute for lime for soil stabilization. As of ecember, 

- 1980, the average bid price in Kansas | per ton of lime and preheater fines was | 

$80 and $40, respectively. Field application of preheater fines can eonvenanty 
be accomplished in the form of a slurry. Preheater fines have recently been 

used for subgrade stabilization i ia at least one highway project in Kansas. — <a 


a XPERIMENTAL INVESTIGATIONS AND EVALUATION ered 


4 


wu 
A series of laboratory te tests was carried out on selected ‘subgrade soils to j 
examine the feasibility of using preheater fines as an alternative admixture for — 
lime. Parallel testings with 3% hydrated lime (subgrade stabilization is a 
based on 3% lime in Kansas) were also conducted to gain data on the relative | 
effect of the two additive materials and a comparative evaluation was = 

_ mainly on the basis of two criteria, namely: (1) Reduction of swelling potential; * 

} and (2) strength gain. The basic characteristics of the soils used in the — 
The variation in the plasticity index (PI) of soils 3, 7, and 8 with percent 
DF cerenoes fines is shown in Fig. 1. With the addition of a ‘maximum 7 

; preheater fines (per dry weight of soil), the reduction in PI of the three soils _ 

d ranges from 52%-65%. Within the range of the percentage of preheater fines 

used for the tests, 6% preheater fines seems to be the optimum percentage. 
| ‘The reduction in PI ranges from 53%-59% when the three soils are treated 

- with 3% hydrated lime. The effect of the addition of preheater fines (up to 
8%) on the standard American Association of State Highway ( Officials (AASHO) 

density and optimum moisture content was found to be minimal and within 

The results of the swelling tests  enlenes on untreated soils as well as : 

soils treated with preheater fines and hydrated lime are shown in Fig. 2. For | 

‘Sr. Geotechnical Engr., National Soil Services, Inc., 5814 Heffernan St., Houston, : 

Note.—Discussion open until October |, 1981. To extend the closing date one month, - 

i a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on September 11, 

1980. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — 

_ of the American Society of Civil Engineers, OASCE, Vol. 107, No. On, May, (1981. 


TABLE 1.—Soil Characteristics 
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FIG. 1 —Effect of Preheater Fines a and Lime o on Index 


| 
| pounds | 
— per cubic | moisture miasa 


a the laterally confined acai. ‘4 in. diam x 2 in. diam (10.2 cm X 5.1 cm) 
at 96 h soaking under 1- -psi (6.9-KN/m7) surcharge after being compacted to 
92% of (AASHO) standard density at optimum moisture content (of treated © 
or untreated sample, as the case may be). It can be seen from Fig. 2 a 
the addition of 6% preheater fines resulted in reduced swelling of three soils a 
to an acceptable level of 1% or less. Comparable results were obtained when 

_ the three soils were treated with 3% hydrated lime. The percentage reduction — 
in n swelling for the three soils 1 ranges from 79%- -86% with 6% preheater fines, = 
 91%- 98% with 8% preheater fines and 89%-91% with 3% hydrated lime. a 
a _ Triaxial compression tests were performed to determine the modulus « of 
deformation of untreated and treated soils. The curing time used was 48 h> 
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PERCENT PREHEATER FINES - LIME 


‘FIG. 2. —Reduetin in | Swell with and Lime 
before. and 7 days” after molding into a 2.8 in. x diam 
(7.1 cm x 15.2 cm) specimen. A traffic coefficient of 10/6, which corresponds 
to a range of 301-450 equivalent 18,000-lb (80-KN) axle load per day, and = 
_ tainfall coefficient of 0.9, which corresponds to an average annual rainfall of 
35.0 in-39.0 in. (88.9 cm-101.3 cm), were used in the determination of the : 
modulus of deformation. A complete account of the testing procedure and the ; 
evaluation of the modulus of deformation is available (1,2). 
_ The increase in the modulus of deformation of soils 9, 10, and 11, when | 
treated with preheater fines and hydrated lime is shown in Fig. 3. The percentage 
_ increase in the modulus of deformation of three soils ranges from 66-131 with — 
6% preheater fines, 138-307 with 8% preheater fines, and 73-91 with 3% ol 
lime. Again, the seems obtained with 6% preheater fines are comparable or 
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NOTES 


" 


—PREHEATER FINES 


PERCENT PREHEATER FINES OR LIME 
‘ FIG. 3.—Effect of Preheater Fines and Hydrated Lime on ‘Modulus of Deformation — 
better than the results obtained with 3% hydrated lime. 
Limited tests on soil 12 were conducted to determine the improvement in | 
_ the unconfined compressive strength when the soil was treated with preheater _ 
: fines and lime. The percentage increase in the unconfined compressive strength - 
of soil 12 was 50% with 6% preheater fines and 25% with 3% hydrated lime. q 


a ‘From the results of a series of laboratory tests carried out on selected et ll 


om, it was found that the addition of preheater fines to a soil sme —— 


g 6% preheater fines are comparable or better than that of 3% hydrated lime. 
The low material cost of same fines makes it economically attractive in 


‘The work reported in this note was done while the writer was with ‘the 
bs Geotechnical Section of the Kansas Department of Transportation. Any opinions, _ 
findings and conclusions or recommendations expressed in this note are those 


| 
have, 


of the writer who is solely responsible for the accuracy of the data | presented 
_ and do not necessarily reflect the official views or policies of the Kansas 
Department of Transportation or the Federal Highway . Administration. | 


. “Design of Flexible Pavements Using the Triaxial Compression Test,’ 

a Highway Research Board, State Highway Commission, Kansas, 1947. _ 

; (2. Worley, H. E., “Methods Used in Kansas to Determine Thickness of Concrete and 5 

Pavement, Bases and Subbases,’ report to to Assistant 
Commission, Kansas, Nov., , 1959. 


By s. Bandyopadhyay,’ M ASCE 


of for natural soils and the identification of 
q potentially expansive natural soils at the outset of an investigation is of essential 
_ necessity in the design of highways, building foundations, hydraulic structures, 


underground utilities, dams, and earth | retaining structures. In the United States. N 


alone, the cost of the damages caused by expansive soils amounts to some 

- $2.3 billion annually—more than twice the damage from floods, hurricanes, — 

tornadoes, andearthquakes(1), 
_ Though there is a considerable amount of knowledge currently available dealing — 

: # with the nature and behavior of expansive soils as a result of elaborate laboratory 

- tests, builders of small structures such as single-family residential buildings 4 
are, not surprisingly, often reluctant to enter into a full-scale investigation 
regarding the evaluation of the swelling potential of foundation soils. In cases 
like this and others as mentioned in the preceding, prediction of swelling potential 
for natural soils based on simple and inexpansive index tests is very — 

"Sr. Geotechnical Engr., National Soil Services, Inc., 5814 Heffernan S ., Houston, 

z Note.—Discussion open until October 1, 1981. To extend the closing date one month, 

a written request must be filed with the Manager of Technical and Professional Publications, 7 
ASCE. Manuscript was submitted for review for possible publication on August 12, 1980. 7 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings _ 


of the American Society of Civil Engineers, OASCE, Vol. 107, No. GTS, May, 1981. a 


The avetege erter m the the is 6. » 


@ A number « of empirical approaches, based on the index tests such as Atterberg 
- limits, activity, percentage of clay, linear shrinkage, and so forth, have been 
_ developed as an aid in this regard. However, as Seed, et al. (2) pointed out, 
y these approaches for the most part have been based on limited experimental 
data and have provided only approximate correlations between index properties 
and swelling characteristics. The results reported herein were arrived at 
be test data of 259 subgrade soil samples collected from 22 counties in Kansas. on 
_ Each of the 259 soils was tested in the laboratory to determine its liquid 
— its plastic limit, its grain s size distribution, and its expansion under a surcharge 
approx 1 psi (6.9 kN/m7’). Placement densities for compacted soils range 
90% 100% of American Association of ‘State Officials (AASHO) 


a . 92% based on on a a standard AASHO compaction test was adopted. The sample _ 
size used in the swelling test was 4 in. diam x 2 in. diam (10.2 cm x 5. * 
cm) height. This leads, for the purpose of the present study, to the following — 
definition of swelling potential. Swelling potential is the percentage swell of 
_ laterally confined sample at the end of 96 h on soaking under 1-psi (6.9- -kN/m* 3 a 
7 surcharge after being compacted to 92% of standard AASHO density at optimum > 
potential obtained for each was then expressed in the fi 
‘inwhich S =the swelling p potential as as defined denotes the percentage 
“of clay sizes finer than 0.002 mm; and K = 


_ type of soil. The activity, A, of each soil defined 
- 72 computed. (Plis the plasticity index. ) Since it has ld been established 
> 


by Seed, et al. that a linear relationship exists between log K and log A,a 


regression analysis to fit a curve of the form 


t the fi ll 
=3.28x 107; b= 2.259; 
# 20) gn > 
2.259 3.44 
‘ Population = 133; a=2.40x 107; = 2.573; 


659 
F 4 
@ 
1 
= It was found that a better correlation can be obtained if the soils were grouped . 
into certain PI ranges. The results of the regression analysis along with the — &g 
| 


<PI< 60; Population = 14; = 0.72 x 107°; 2.669, rw095; 


TABLE 1 —Compured ‘and Experimental Swelling Potential 


| Swell, as 
-:percentage 


& 


t 


35 heat 0.867 


Dickinson 


Bs 


47 


on 


together, then the values of a, bya ,and rare 2.2 x 10~°, 2.437, and 0. 78, respectively. | 


The accuracy y of predicting the swelling potential by the developed equations — 
can amply be demonstrated in Table 1 in which the predicted swelling potentials 
4 


_ Eqs. 4, 5, 6, and 7) for 26 soils obtained from eight counties in Kansas 
"were compared with the swellings determined in the laboratory. These soils — 
were not included int the le regression analysis. be Seen from Table 


| 
| 
ne hig value OF correlation icient in each case Suggests a hig egree 
| Liquid | per- | 
&g County number limit PL centage Ae PI/C ested 
= | @ | |] © | 
Anderson | 1C-2 | | 21 | 38 
ais | om | 4 
| 25 | 33 
9 | 2 | 100 | 18 | 18 
Sh 29 | 39 | 
Lyon | a1 | 46 
! 


there i is an excellent between the swelling and and the 
swelling. The The average err error in the prediction is 6.9 9%. 
‘Summary ano CONCLUSIONS 
_- Using test data « of 259 naturally occurring ‘alta obtained from 22 counties _ 
4 in Kansas, a well-defined relationship was established between the percentage < 4 
clay sizes present in a soil, the activity of the soil and the percentage of swell 
5 under a 1-psi surcharge of | a sample compacted to 92% of standard AASHO 
density at optimum moisture content. It has been shown that there is an n excellent . 
- "agreement between the predicted swelling and the actual swelling. The average 
error in predicting the swelling potential for the 26 soils was found to be 6.9%. 
_ prediction equations will serve as a valuable tool in the design of any — 
structure on c acted soils. i od! 
a The work reported in this note was done while the writer was with the 
Geotechnical Section of the Kansas Department of Though 


Jones, D. E., and Holtz, w. ‘Soils—The Disaster,” 


“a Engineering, ASCE, Vol. 43, No. CE8, Aug., 1973, pp. 49-51. 
‘ ~ Seed, H. B., Woodward, R. J., Jr., and Lundgren, R., **Prediction of Sweling F Potentiol 


v4 ba 


a - _ solely responsible for the contents and conclusions of this note and the accuracy bil _ 
7 Z _of the data presented. The contents do not necessarily reflect the official views — 
— 
| 
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oF SEEPAGE IN — Dams 
— | low ¢ oi 
Yang H. Huang,' F. ASCE 


_ mathematica! solution for determining the line of seepage through an earth — 
: Er an impervious base was proposed independently in 1916 by Scheffernek = 


in Austria and Iterson in Holland and was described by Casagrande (1) as applicable 
an  outslope not over 30°. However, the method is valid only when the 
impervious base is horizontal and does not apply to earth dams on an inclined 


ledg<. The latter case is particularly prevalent in eastern Kentucky where the 

overburden is thin and earth dams are mostly constructed on an inclined ledge. 
The purpose of this paper is to present simple charts for determining the line __ 

of seepage in earth dams on inclined ledge. After the line of seepage is determined, _ 
4 


Fig. l an on an inclined ledge. The downstream face of 
the | dam has a slope of S:1 (horizontal:vertical) ; and the Soec 

; base, has an angle of inclination, a. The height of the dam measured from 
the pool level to the toe is h. Following Casagrande’s procedure, it is assumed 

_ that the theoretical line of seepage starts from the pool level at a distance = 
_ of 0.3 4 from the dam, in which A = the horizontal distance shown in Fig. | 
od . Therefore, when the toe of the downstream slope is used | as the origin of 

7 coordinates, one point on the line of se seepage > with x - dandy = his eee. a 
Assuming that the coordinates of the exit point and the mid point are (aS, a) 

; and [(aS + d)/2, b], respectively, it is now necessary to determine a and 


—é, so that the coordinates of these two points can be located. Through the 


'Prof., , Dept. of Civ. Engrg... Univ. of Kentucky, of Engrg., Lexington, Ky. 


Note. —Discussion open until POEs 1, 1981. To extend the closing date one month, > ; 
a written request must be filed with the Manager of Technical and Professional Publications, ~ 
- ASCE. Manuscript was submitted for review for possible publication on November 7 
1980. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, ‘No. GTS, May, 1981. — 


ne 
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and tangent to the downstream silane. the theoretical line can be slightly | a 
to fulfill these boundary requirements. 
7 _ Based on Dupuit’s assumption that in every point ona vertical line the Getic 
ad gradient is constant and equal to the slope, dy/dx, of the ope} of seepage, 
the the can be expressed t by s law as a 
=k(y xtana)—— . 


in which g q= discharge; k = permeability; « a= of of base; 


ons x, = = coordinates. For integration purpose, Eq. 1 can be written 


#3 


which x is considered a a variable; and yis an 5 


‘ad variable. E 2 can be readily integrated a: 


The constant of in substitutin 


y = h into Eq. 3. At the pc 
ka(li— Stana) 


=- 


mas 3c can be written as 


«FIG. 1.—Earth Dam on Inclined Ledge 
ec. 
x = d and 
(4) 
~ ona) 
x= 


SS 
02 
= 


a 2.—Charts for Determining Point of Exit ‘ 


ral 


06 
| 
_ 
COON 
| 


7 can be solved numerically by: means of Mueller’s iteration scheme 
of successive bisection and inverse parabolic interpretation (2). The results are 


wide’. The lack of Pit the sans is puch dams 


MIDPOIN 


05 


4 


resented in Fig. 2. After a/h is Quen b/h can be determined from Eq. — 


| 
- 
| 
_ | }s-25) HH a 


tan f tana 


\h 


ott can be seen from Figs. 2 and 3 that the line of seepage for earth oa 
= an inclined ledge is than that on a horizontal ledge, 


a application of these two charts can be illustrated by the following ppunengee: oa 
Fig. 4 4 shows the cross section of a dam. As ‘sown in the d= 


d= 390 ft (8.9m) 
FIG. 4. Example 
ft (118.9 m) and h = 150 ft (46.7 m), thus ae 2.6. From Fig. 2, when 
‘= 15, d/h = 2.6 anda = = 6.7 7°, so a/h =0. 29; when S = 2, ash = 0. 
‘Pee the ‘outslope « of dam is 1.75, a/h = (0.29 “i 46)/2 = = 0.38, or 
= 57 ft (17.4 m). From Fig. 3, when S = 1.5, b/h = 0.73; when S = 
q b/h = 0.77; so b/h = (0.73 + 0.77)/2 = 0.75; or b = 113 ft os 4 of dtd 
coordinates of the three points, as indicated by the cross, the line of 
of only approximate results are required, the line of seepage can be determined — 
by connecting a straight line through the two end points, while neglecting the en, 
midpoint. This will result in a line of seepage which is too low and is on the 


| 
1 imple Charts are presented for earth dams | 
4 


« an impervious: base. The charts can be applied to a ew base as well 
_ as a sloping base. The laiter case is of particular si significance because in a rugged 
4 terrain many earth are constructed on inclined ledge. Since the 
_ promulgation of the Surface Mining Control and Reclamation Act of 1977, it © 
,% ‘ is necessary to determine the factor of safety of every coal refuse dam. These > 
_ dams are mostly constructed on an inclined ledge with no internal drainage 
+ provided. The lack of a simple method to determine the seepage in such dams 
ot has posed a great difficulty j to engineers in making stability analysis. It is hoped 
. that the simple method presented herein will be found eats for the stability 
of these dams. 3 
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By Massimo o Magni’ and Ale ex P. Mi 


INTRODUCTION 
bulb i 
The discovery of small | reserves of g gas does not, in general, ‘economically 
: justify the investment required for construction of conventional offshore plat- — 


"Mgr. of Servizio Offshore, Agip S.p.A., Milan, Italy, 
Vice Pres. and General Mgr., D’Appolonia Societe Anonyme, Boulevard Du Souverain 
7 Note.—Discussion open until October 1, 1981. To extend the date one 
a Written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on June 10, 1980. 
_ This paper is part of the Journal of the Geotechnical Engineering Division , Proceedings 
of the American Society of Civil Engineers, O©ASCE, Vol. 107, No. GTS, May, 1981. 
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platforms for drilling and producing gas wells, in water depths of the order 
of 100 ft (30 m), have been described by Magnanini (3). A collar caisson was: 7 
considered in the design of such monopodes, in order to decrease excessive 
a deflections at the mudline level when subjected to lateral loads (Fig. 1). The 


a 


= 


= herein was out in order to assess: ‘the influence of the 


collar caisson on the lateral deflections at mudline. ory es? 

™! case investigated is a typical monopode platform | for clay soil conditions | 


which are characteristic of the Adriatic Sea bottom, offshore Italy. The cntons “a 
is driven in the soil to a depth of 148 ft (45 m) and acts as a foundation 


q 
q 
View of Monopode with Collar Caisson | 
— i 


ys 772 


m | 


‘ (2 


14.0x0. 


bi 

@ 13x0.062 ft 
(4.0x0.025m) 

13x0.096ft 
(4.020.03m) 


91320. ustt 
(4.0x0.035m) 


(30m) 
75.5 ft 


a 


an 
16.4ft 


13x0. ust 
(4.0 x0 .035m) 
fol 


§13x0.082 
4. 020.025m) 


= 


(4.020.02m) 


= 


= 


soil strengt 
(8 m) depth of the caisson, , and then to increase 


a TECHNICAL NOTES 


2 2. .—Monopode and Soil Properties: 


ig fry ed) od] 
- S. was assumed to be constant for the first 26 ft 


of 


A 
fos 


STRENGTH 


PROFILE | 


bai, (KN Dodi 101 


linearly with depth, reaching 


N's pile. The diameter of the monopode is 13 ft ( ; -ft (8-m) a 
caisson is placed arout (8 m) below the mudiine. 
Pat a depth o — | 


KN) was applied to the monopode at 75 ft (23 m) above seabed. _ 


MAY 1981 
below mudline, i.e. at the bottom of the monopode (Fig.2), 
_ The monopode and caisson were modeled (finite elements) using a1 Roi 
elements, to be analyzed by ‘““DAPSYS”’ (1) as explained later, with the cf 
ay 
_ properties of steel and concrete | to represent the monopode and the caisson, _ 
Be ams Three different sets of section properties (areas and moments of : 
inertia) were used for the embedded, submerged and free-standing parts of 
the monopode. A different set of section properties was used for the elements 


K=4 


in which G = shear modulus of the soil, and taken as: 700° 


S,/2(1 + v) (Poisson’ s ratio vy 0.45); and AL = tributary length to each 


“ The spring constants were proportional to S, and were functions of the depth. - 


7 The horizontal and rocking stiffness of the caisson was also modeled with spring 


elements, where appropriate. These stiffnesses were calculated as the surface 
. stiffness of a rigid disk (4) with a 26-ft (8-m) diam and modified to | account g 
_ The system was analyzed statically for the cases with and without caisson 
“and for varying soil strength (and therefore stiffness), S,, using the computer 
« code ‘ ‘DAPSYS” (1), a general purpose finite element program developed by 
D’ Appolonia for static and dynamic analysis « of structural systems. Shear forces, 
moments, and deflections were computed at various locations along the monopode 
length forthe casesexamined. 
It should be noted that the assumption of linearity and the simplified relation- 
"ships such as given in Eq. 1 and between G and S, would not in general yield 
‘reliable results if design values are based on these. However, this simplified — 


‘procedure was adopted, for purposes of comparison only in order to assess 


the effect of the caisson. 


_ for the effect of embedment (2). A static horizontal load of 582 kips (2,590 7 


increasing values of soil shear within the range of Ss. values 


- The results obtained may tie seen i in Fig. 3, in which the ratio of the monopode 
lateral at mudline with caisson versus that without caisson is plotted 
= S,, in the range of values considered in the analysis. The following: 
are made from Fig. 3: 
d 


considered, the deflection ratio decreases monotonically, i.e., presence of 
caisson reduces the deflection at mudline. 
= The deflection ratio of 0.8 corresponding to soil shear ani equal to ics : 
ro , which is equivalent to no soil, results from the difference in deflections 


= to the change of stiffness in the cross sections with and without the caisson. 


| 
: 
— - stiffness constant of these spring elements was calcualted using Richart, et a 
: 
| 
I 
| 3 


TECHNICAL NOTES 


- 


> 


DEFLECTION RATIO AT MUDLINE — 


(WITH CAISSON/ WITHOUT CAISSON) 


. 


40) (60) (160) (140) (160) (90) 


FIG. —Lateral Pile Pile Deflection at Mudline with and w without Caisson 
nil is furthermore noted that the linear treatment of the problem, somewhat — 
underestimates the effect of the caissca, in that the G for a smaller deflection _ 
4 (with caisson) would be larger than that for a larger deflection (without caisson). a 
7 ® It may therefore be concluded that the caisson is effective in reducing the 
mudline deflections if the soil shear strength at the top layers in the vicinity ” 
_of the caisson is not very low. ys reduction by a factor of two is obtained 


a typical value of S, of 20 kN/m’ 10 


1. “DAPSYS: A Computer Code fer Soil-Water-Structure Interaction Effects,”’ D’Appo- | : 
tonia Consulting Engineers, Inc., Proprietary, Brussels, 1977. 
2. Johnson, G. R., Christiano, P. 4: and Epstein, H. I., “Stiffness Coefficients for _ 
_ Embedded Footings,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 
No. GT8, Proc. Paper 11522, Aug., 1975, pp. 789-800. 
3. Magnanini, U., ““Monopode Platform for Drilling and Producing Directional Wells in 
the Adriatic Sea in Water Depths of up to 30 Meters,” Proceedings of the Congreso 
Panamericano de Ingenieria del Petroleo, Mexico City, Mexico, 1979. 
4% Richart, F. E., Jr., Hall, J. R., Jr., and Woods, R. D., Vibrations of Soils and } Foundations, 


Prentice- “Hall, In Inc., .» Englewood Cliffs, N.J., 197 1970. 


‘LATIONSHIP 


Ralf Pek’ 


INTRODUCTION 
A description of rock-joint roughness i: is ite rucial importance in shear strength - 
“estimation. Several methods of measuring joint roughness (1,2) have been — : 
developed which involve the use of a profilograph and evaluation of the statistical _ 


4 properties of the profile. However these are time consuming, and the results _ 


are often highly dependent on on the particular profile chosen. 
_ The method presented in the following samples the roughness of an area, 
rather than a line, on the joint surface by measuring the area of shadow cast - 
_ by the asperities when the joint is illuminated by oblique unidirectional light. 
_ The data thus acquired oy be used directly it in a theoretical model to estimate 


; she ear strength. wid 


joint is by enidirectional light, such the 
= the joint plane on n the light rays is i, as shown i in | Fig. Ma). 4 


The area of shadows can then be measured as a fraction a of the total joint - 
area. This is best done on a high contrast photograph taken with the camera - 
axis perpendicular to the joint plane. In some cases it may be necessary to Sak 
Paint the joint surface white, so that a photograph with only white and black — 
tones can be obtained. 1 The fraction of black area - from such a photograph — 
is easily estimated, some possible techniques being direct measurement,  plani-- 

_ metry, electronic scanning or measurement of the amount of light which ‘passes 
through a high contrast negative slide. 

a Measuring the area of shade for various values of the angle of incidence 

i gives the roughness function a(i). This function is in a sense a — 


distribution of steepnesses and contains more information than any single 
parameter. To illustrate ‘this it is noted that, at a given normal stress, small 
‘steep asperities and larger or more numerous smooth ones may be equivalent 
in terms of their contribution to shear strength. However, as the normal stress — 
is increased, the small steep asperities will tend to be sheared off, and their 


-., "Research Asst., Dept. of Civ. Engrg., Univ. of Canterbury, Christchurch 1, New 

_ Zealand; now, Grad. Student, Mass. Inst. of Tech., Room 1-245, Cambridge, Mass. 02139. 
Note. Discussion open until October 1, 1981. To extend the closing date one month, 

a written request must be filed with the Manager of Technical and Professional Publications, 

 <.: Manuscript was submitted for review for possible publication on April 1, 1980. 

“This paper is part of the Journal of the Geotechnical Engineering Division, oe 
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TECHNICAL 2 
Light 


Direction of 


2 bed doidw af 


Lower Joint Pace under Unidirectional Light Incident at Angle i; 
of Failure; (c) Surface Tractions on Upper Joint Face at Failure 
_ contribution: to shear strength will be pares" by that of the larger or more 
numerous asperities. Therefore, unless certain similarities in the distribution 
= of steepnesses app apply, a one parameter description is specific to a given normal 


3 While this method provides a more representative and more complete descrip- — 
F tion of roughness, the main reason for proposing it, is that the optical roughness" 
data may be used directly to obtain sl shear strength estimations using the failure 


= In developing a model for joint failure the following « ae are remade: 


— 
= 
— 
| 
4 


‘The strength of the rock at the joint may be described 
by a Mohr-Coulomb failure envelope with parameters c = cohesion and @ = 
2. The effect of gaps between the joint surfaces is 
3. Failure occurs by shearing through asperities along failure planes, on which 
limiting Mohr-Coulomb stress state is developed. 
| 4. These failure planes develop in one joint face only, and are at an angle 
i to the overall joint plane orientation. (The angle i need not be assumed, = 
_ will be determined by the condition that the shear strength isa minimum.) - 
_ §. Forces are transmitted across the joint surface entirely through failure 


7 The conditions at failure are shown in Figs. 1(b) and l(c). If the area ws al 
failure planes is a fraction b of the total joint area, then equations of horizontal 


and vertical equilibrium of the joint surface shown in Fig. l(c) may be written 

S= 

‘and + cbs sin (i + cos (i+ 

in in which N and S = the normal and shear stresses transmitted across the joint 


% surface, respectively; and g = the magnitude of the noncohesional component oe 
the stress vector acting on the failure planes a as shown in in Fig. 


q from Egs. 1 and 


= cb (cosi + sini tan (i +) + Ntan(i+) 


Bq. 31 4 


$= cal + tanitan (i+ + tan (i+ 6). 


in which a = cos 


_ It is seen from Fig. 1(c) that ais the ‘horizontal projection of the area of * 


failure planes expressed as a fraction of total joint area. In addition comparison _ 
_ of Figs. 1(a) and l(c) reveals that a also equals the projection of shaded area es 
_ when the joint surface is under unidirectional illumination at the angle i. Thus _ 

a is determined “experimentally by the previously defined “roughness function 

 a(i), and Eq. 4 can be plotted on the S-N plane for various values of i as 

_ Straight lines with slope angle i + and ordinate intercept c a(i) [1 + tan 

i tan (i + )]. Since failure occurs at the angle i at which S is smallest, the 


envelope of such lines gives the shear as a of normal st stress 
__N. This process is shown in in Fig. 

g In the previous section the shear ‘Strength: was s estimated from the optical 


- may be argued that | to prong c - and is more e difficult than testing the joints 
themselves. However, in practice, the angle of sliding friction, as obtained for — 
example from a friction test on flat, saw cut rock surfaces, is a good approximation 
for the parameter it is noted that, most jomt surface 


| 
2 
qj | 


cohesion For example, using roughness data shown in Fig. 2 and ll 
. the percentage increase in shear strength S due to doubling of the cohesion - 
shows that, for normal stresses between and 1.0c, the shear strength 
doesn’t increase by more than 25%. Therefore it is sufficient to estimate c 
| applies 


A 
q 


4 


FIG. 2. —Construction of Failure Envelope 
4 defined previously it samy shown that: st 
which the unconfined ‘compressive strength of the joint wall material. 


Ls in direct shear at different normal loads. For this “purpose | it 1 was chosen to 
vs use 5-1/ 4-in. x 4-in. (130-mm x 100-mm) plaster replicas of the same perenne 


| 
web pe 4 own 
coghtant 
> | necessarv veral nt I 
| 
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—Predicted Shear 
with c_=0.1c 


SHEAR STRE 


| Shear as Function of Norma! Load and 


ranging from = 0.0006-N /c = 0.04. Confined and unconfined compression 
tests on cylinders of plaster showed that it did behave as a brittle material 
had well defined Mohr-Coulomb strength parameters c and However 
_ the use of these parameters in the foregoing model lead to an overestimation 
of shear ae The main reasons sons for t this are thought to be the ve following: t. |. 
1. The assumption 3 that the limiting Mohr-Coulomb stress state is developed 
on the failure planes is not valid for a brittle material such as plaster or rock, 
4 since it is incapable of undergoing the plastic deformations required for redistri- — 
_ buting the stresses. Therefore failure occurs progressively and the cohesion 
. c is only developed on a fraction r, , of the area of failure planes. 
There isa phase’ of sliding with loss of interlock and little e dilation before 
_ peak stress is reached. The effect of this is that the e actual area of failure 
planes is only a fraction r, of what was predicted froin the roughness description. 7 
_ 3. Asperities which are smaller than the shear displacement at peak load 
are either sheared off at an early stage or never developed due to imperfections »| 
‘in the interlock and therefore cannot contribute to shear strength. ‘This 


Tom small asperities would be predicted. 


can shown that adjusting as 1s suggested for the fir t 
and effects is equivalent to 1 to the cohesion by a reduced effective 
in which r=r neon 


9 In Fig. 3 the failure envelope has been plotted taking r = 0.1, leading to 
3 good agreement between the estimated shear strength envelope shown by the 
- continuous line, and the experimental points. This ‘suggests that, in the range 
of normal stresses tested, the correction factor r remains constant with normal 
‘stress. ‘Therefore, for any particular joint, r could be determined by testing 
the joint in direct shear at any convenient normal stress, and the result could 
_ then be used to predict the entire shear strength envelope. Further experiments 
ere show that r is constant for a wider variety of conditions. a 
Tf further refinement is desired, the third effect could be allowed for by 
‘ignoring areas of shade whose dimension in the direction of shearing is small 


compared to the displacement ‘at peak odt Yo 


A method of describing the roughness of joint surfaces has been edie 8 ; 
which samples an area as a whole, leading to a roughness function which is a 
specific to a given direction of shear. The acquired roughness data n may be i 


J 
. used i in a rational model to estimate shear strength. Although a correction factor Pa 


r must be introduced to allow for progressive failure of asperities and loss 
of interlock, this correction factor has a physical meaning and estimates of 
shear strength are not sensitive to the value of this correction factor. _ ——— 7 
The results of direct shear tests on Plaster replicas < of the same joint surface 7 
show good agreement between measured and pre predicted shear: strengths at different 7 - 
normal loads if a suitable value for r is chosen 
_ If further experiments show that r remains constant ant for a larger range’ of | 
i stresses and a variety of joint surface geometries its value will have a 
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7 _ The writer is indebted to the Central Laboratories of the New Zealand Ministry f§ 
s _of Works and Development for having provided the facilities to carry out the e 


vr aes 


Durinc EARTHQUAKES 
__ This note presents a summary of a simplified procedure (using hand ns) ng : 


for: preliminary estimates es of pile. maximum curvatures (bending stresses) during = 


a 


_* The main assumptions of the proposed procedure are: (1) Shear stresses ~ 


strains in the soil are the result of vertically propagating shear waves jo 
r underlying bedrock; and (2) piles move in phase with the soil. The first sumption 


away from epicer:ers based on evidence cited in Ref. 4. Available pi 
cited in Ref. 4 suggests that the second assumption is reasonably justified for 
— long or flexible piles. The second assumption implies minimal contribution of 
- piles to horizontal stiffness. This may result in an overestimate of pile curvature ~ 
near the pile cap since s soil-structure interaction can be expected in this zone. 
_ Thus, the proposed procedure may be used for evaluating curvatures at depth 
“and not near pile caps. For evaluating ane sections near pile caps, procedures 
as | in Ref Ref.3 may be used. 
‘The basic steps in the proposed procedure are: (1) An estimate of seismic 
shear strains using strain compatible shear moduli; and (2) an estimate of curvature | 
_as the rate of change of shear strain versus s depth. A logic diagram illustrating 
various calculation steps is shown in Fig. 1. $= 
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| 
i 
: » the importance of the situation and the results of preliminary analyses. Other = a 
mmmedspects of pile design, such as stresses induced by inertia of the structure, 
— ence, or soil failure are equally important and should be considered. Details - 
le development of the proposed procedure are presented in Ref. 4. It is 
| 
&§ 
q 


Estimate of Shear Stresses. spectrum of 
* rock- soil interface for a modal —— tatio of 0.1 and the site period, — 


Sa,, Sa, = acceleration response spectra of input motion 
orresponding to the mode i (i = 1, 2) having a period, T,, for a modal damping 


io of 0.1. The value T, ed site period, T,; Sa, is 


Me, Compute Sa, from T, 


Estimate Sap as eye | Estimate 1 and the 
distribution Tz, and 


Define basic Parameters} 
of Soil Profile includ- 
ing and small 
strain V and G max 
=) 


wit of G by entering a 


[Select reduction curves fraction (0.65 or 0.7)} 
into reduction 


based on an eye average for the period range T, = 0. T, as 
in Fig. 2; H = thickness of soil profile; o,, = total vertical stress at the base 
of soil profile; and g = acceleration of gravity, = 
The preceding two mode approximations for shear stresses is based on ‘many — ; 
numerical experiments with computer program SHAKE (7) described in Ref. os 
_ 6, The one-dimensional solution utilized in SHAKE assumes a free boundary. —_— 
’ Heavy pile caps or embedded basements may significantly violate the yn alll 
conditions in SHAKE-type models. Eqs. 1 and 2 assume the response spectrum = 
is specified at interface. If rock motion is specified at rock outcrop, the interface a ' 
spectrum can be calculated by : scaling down the outcrop spectrum (4,6). 2 
Aan Shear Strains. —Fig. 3 shows average relationships of secant shear modulus 
degradation curves for sands and clays (8) derived from laboratory tests replotted | 
in terms of shear stress, t, shear modulus at small strains, G,,,, , and modulus, 
_ compatible with shear stress, . Fig. 3 can also be in terms 


shear velocities V and Vin vax BY noting that V/V max Since the 


A max max 


sheer stress, Tz ‘computed in wed 3 is a maximum, a fraction of this maximum — 


| 
at depth, Z, beneath the surface are computed based on: 
Bes 
q 
4 


Sa, = 0.372 x 0.33 g= 0.128 


Sa, = 1.4.x 0.33 g = 0.468 


FIG. 2.—Selection of Sa, and Sa, for Typical Case » (Peak Acceleration « a, = 0.33 
g; Predominant Period Rock Motion, T, r 0.32 sec) 


Average for clays 
12 Averege for sonds 


(0.65 +z oF 0. - tz used in SHAKE) may be selected to represent the ave 
cyclic shear stress for entering Fig. 3. 


_ With the estimated maximum shear stress, «| a. at midpoint of layer 0, 
modulus, known f each soil layer, the shear strain, ¢,,is simply 
Estimate o of C urvature. —Using t the m method of difference and 


as ayer thickness as H,, the curvature, C,, can be estimated by — 


| 
Bil 
forks 
4 


— 
5 employs the maximum curvature, Ci i Eq. 5, therefore, 
is not necessarily maximum curvature > developed during ground 


A curvature time history ; analysis w would be r required to assess an analytically 
correct maximum curvature at midpoint of each layer, i. A reasonable er 


bound on C, may be established assuming zero strain in layer i— 1, and 
which n = the number of soil layers. ig 


The ability of the soil to transmit shear waves is by 


‘4 


(if the shear strength is zero, wave velocity is zero, and shear waves would - 
not pass | through | the soil layer). ‘Thus, if the calculated shear stresses are near 
the ‘undrained shear s strength of the soil, the simple analysis presented herein 

: = for ‘strains (and curvature) is no longer valid , and a detailed analysis would © 

7 be required. The modulus reduction curves similar to Fig. 3 are usually based | 

_ Example.—Assume a layered soil profile consisting of 11 layers with thickness " 
and modulus (H,,G,,,,) given in Tabie The applicable modulus 
. is: curve 2 in Fig. 3. The design rock 1 motion response spectrum at soil- rock 
- interface (bottom of layer 11) is shown in Fig. 2. The site period, 7; = usually — : 
unknown at the outset and should be estimated (see e.g., Refs. 5, 6, 10, 1 1). ] Py 
_ For this example problem, the site period 7, was estimated to be 0.45 joc. 
using the procedure described in Ref. 6. The site period is dependent on seismic “a 
- shear strain, and therefore, the characteristics of rock motion. == it 
_ Based on the preceding data, the following calculations are made: (1) Sa,/a, 

‘ = 1.39 (from spectrum in Fig. 2), in which a, = peak base rock acceleration; 
and T, = 0.3 T,-0.5 T,; (2) Sa,/a, = 17 (eye average from Fig. 2; 3) o,, 
= 3", = 0.545 x 10° psf (261 kN/m’); (4) Sa,/Sa, = 1.7/1.39 = 
(1.22 < 4; (5) 10,,(Sa, /g); and (6) t, = 0.7 x 0.546 x 10° x 1. 39 

- The rest of the calculations are shown in Table 1. The parameter + 1/G,, » 


7 (= 0.65 +,/G,,,,) for each soil layer is entered into the reduction curve od 
4 estimate G,/G,,,,, and thus G,. Strains and curvatures are then calculated using 


The curvatures computed 6 by the proposed procedure are ws — to compute 
pile bending stresses: Js 


c= 

in which d = the radius to the extreme fiber of pile « cross: section; El = yy “a 
flexural rigidity; o os = extreme fiber bending “stress; J = = = moment of inertia 
of the pile; and M = bending moment. 

.. is clear from Eq. 7 that as the pile se section increases, M (and o) increases 


“soil Thus, the most suitable pile would be a flexible pile with smallest section 
gg would satisfy other design Tequirements (e.g., , vertical | capacity), » and b have 


_ preceding ee were applied to several cases in which a solution 


| 
| 
| 
4 
4 
= 


his 


number, thickness, mid-layer, | 


4 
3 


= 


based on the SHAKE program was known. Details are presented in Re 
: - In all cases, the largest curvatures which control design are reasonably predicted. a 
py _ The stresses associated with bending moment calculated be. Eq. 7 should 
) be considered along with other sources of stresses such as ; those produced — 
by inertia of the structure and static stresses. As n noted earlier, forces induced 
inertia of the structure likely to the near pile caps. 


procedure is proposed at depth during ground 
om motion can be rationally (not empirically) estimated byh hand calculations. Results 


. using the proposed procedure are in reasonable agreement with detailed analyses 
using program SHAKE. The overall approach has not been tested with respect — 


pile based on comparison with actual measured data. With the introduction 
of the proposed procedure, the pile problem | becomes a geotechnical as well a 


to its ability to accurately predict earthquake-induced bending stresses in a 
a 
¥ 


structural problem. The character of rock ‘motion, soil profile, and soil 


‘The writer acknowledges with | the constructive comments and 


criticisms made by ASCE reviewers: helped improve the of this 

Appenoix. 

Margason, E., “Earthquake Effects on Embedded Foundations,” Associated 
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1.—Calculation of Curva- 
4 
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0.082 x 
(0.16 x 107° 
«1.14 107° 
1.71 x 107? 
6.25 x 107? 
0.95 x 107? 
10.85 x 107 
x 107 
io? 


‘f 
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a 0 


8944 | 10° 0.81 | 7,244 
8,944 | 8.69x 10° | 0.74 | 6,618 
3,354 2.59x 10% | 0.35 
(3,354 | 3.17 x 0.29 (978 
(2,674 4.61 x 10 0138 | 348 
3,002 487x10-* 0.12 | 360 
(6,708 233x10-* =| 0.38 2,549 
2,910 | 5.61x10* | 0.08 | 233 
6,708 | | 0.36 | 2,415 
= 
— 
: 
| 
| 
| 
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FINITE STRAIN CONSOLIDATION: Test 


which e ratio; k = the coefficient of permeabili ty; 
effective stress; and ya and = the unit weights of the solid and fluid, 
7 respectively. The coordinate, z, is a reduced coordinate encompassing a volume oo 
4 of solids in a volume of unit cross- “sectional area gh; between the ie 


, is related to the reduced coordinate, z, : had i 


in which e, = “the initial void ratio. Tt e upper /lower sign is ‘ities if z is 

3 Eq. 1, along with the appropriate boundary and initial conditions, forms a 

‘complete mathematical statement me of one- dimensional f finite 


aia. and permeability tests within the framework of finite strain — 
5 in od! 


_ A constant sae of strain (CRS) oie is performed by imposing a constant — 
win'o +v,, at one boundary (6,7). The other boundary is held fixed . The 
plus/minus ais is taken if v, is imposed against / with gravity. Two cases apply 
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‘TECHNICAL NO NOTES 


= one can be impervious, while t gations; or both 
One Impervious Boundary.—The test conditions for this case are schematically | 
shown in ined 1. At the impervious surface, the governing condition is ee 


GB) 
in v, = the absolute apparent velocities of the fluid and solid 
respectively. Applying this condition to Darcy’s law 


and subtracting | the result from | the equation « of ~~ ey for the bulk (solid 


in which u,, = the pore water pressure. It is noted that the consolidation process 


The imposed velocity of the is During the test, the 


remain constant but the ‘between the soil and water 7 


aa FIG. 1.—CRS Test with One Imperv: 


en, € boundary condition for an impervious surface =, 
water 


Using Eqs. 4 and 5 results in 


and the absolute ponassho velocity of the solids, v,, will be yf to v, at 

_ the soil water interface. As a consequence of | of this, ‘the apparent relati relative velocity 


_ which, along with Eq. 8, form t the boundary conditions for this | test. Pe nx 

7 srt 

It now remains to show that the conditions of Eqs. 8 and 10 conform to’ 
the CRS boundary condition when the soil is assumed to re without self-weight 

(a thin sample), and the strains are ‘infinitesimal | 6, 7). When the - ‘material is 


a without _ weight, the term y, — Ve can be taken to be zero wales 8 


- strain (5) by noting that the Eulerian cosatiaabn x, is equal to the Legrengion a 
“hth a, and that | + e is set equal to | + e,. Furthermore, it is assumed 
that me soil behaves as a linear material i in which the oeagonaael « consolidation, im 


‘ax x 


in terms of This equation reflects a m 


strain, €, for the CRS test is is expressed by (1) 


in which F(x t)= 


” Taking the derivative of this expression at the boundary x = 0 and converting | 


r to the equivalent velocity, v,, leads to Eq. 11. Thus, the moving finite ‘Strain | 
| to the equivalent condition for infinitesimal ‘Strain, 


bas zero. Both boundaries would then be governed by Ea. ll with v, set to | 


zero. The initial conditions for relaxation is the void ratio distribution at the © 
= 


test will is ina water pressure redi: redistribution within fixed 


2 
_ _ compressibility. Using Eq. 2 and the preceding, Eq. 10 red q 


the consolidation process is monotonic. This premise is violated in the relaxation 
Process. oft 19 bas ai wort od! te of tine 
Two ‘Pervious Boundaries. —If the CRS test occurs under doubly-drained 
conditions, then a "modified set of boundary conditions becomes necessary. _ 

Here we seek the relative | velocity between the fluid and the solid at each © 

~ boundary. Thus, the form of the boundary conditions follow Eq. 10, for both © 

boundaries with v, replaced by v, and v,, the apparent relative velocities “4 

: s Assume that end (2) is fixed and end (1) is moved at the imposed velocity, 

. The apparent relative velocity at the fixed boundary, v,, will be equal 

P- the velocity of the fluid surface connected to this end of the sample. The 

water surface connected to the other end of the sample will move at the same 

‘rate, v,, since the total amount of material (soil and water) in the system oa 
constant. Therefore, the apparent relative velocity at the moving boundary, _ 

= , will the between the imposed and 


Calculating the ratio, and knowing the imposed \ velocity, v 
velocities v, and vs can be calculated by Eq. 13. These velocities are = then 


ie measurement of consolidation properties entails development of both 


coefficient of k, and the void ratio. Tt is desirable to 
_ the permeability measurement during consolidation. Furthermore, for clays, a 
_ constant flow rate test (4) provides advantages in speed and accuracy which © 
make it more desirable than the usual falling head test. Each permeability test: - 
will trigger a consolidation process. This process will be most pronounced in 
_ soft materials such as marine sediments, mine wastes, and 1 dredged fil fills. ‘Thus, 
it is essential that the test enalysis couple permeability and compression testing. 7 


D De am Die Me amnpnie w while De 
| 
Furthe q 
¢ 
od) 

ii 
1 
: _ test. The relative velocities at both boundaries will be equal to the imposed 


4 flow velocity. Then, v, on the right- -hand side of Eq. 9 is replaced by y a/A, 
in which A = the cross-sectional area of the sample. . Consequently, both 
boundaries are governed by Eq. 10 with v, , replaced by q/A. The sign of q © 
will be positive at the surface where flow in input, and minus at the outflow 
1 Next, we consider a flow rate test conducted during a CRS test. In this 
case, the right-hand side of Eq. 9 becomes qQ/A at the fixed boundary, and 
. + q/A at the moving boundary. Eq. 10 represents the boundary conditions 4 : 
with the appropriate replacements. In en this test so as to avoid swelling, 
“Lastly, we examine the boundary if a defined flow rate is superim- 
posed on a standard (constant stress) oedometer test. At the boundary a 
flow is input, the right-hand side of Eq. 9 is replaced by qQ/A. Consequently, ie 
Eq. 10, with g/A replacing v,, is the condition on this surface. At the other 
_ boundary, the void ratio is specified. 


, ‘The preceding analysis has provided the boundary conditions required t to 
perform an analysis of strain controlled consolidation / permeabi lity tests where 
the finite strain consolidation theory is used. It should be noted that the 
 sheteaiies between e and o’ and between k and e are in the boundary conditions, 
as s well as the governing equations. This multiple dependency complicates the 
usec of the CRS test to determine constitutive properties. We leave an examination 
of this to a later paper. 
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STABILITY . ANAL YSIS OF EMBANKMENTS AND SLOPES 
Discussion by Robin N. Chowdhury? 
al literature shows th: that there is is coutiderable 
tance to identify and emphasize unresolved issues and unsolved problems | 
concerning the analysis of slopes. These include: (1) Initiation and growth of 
he slip surfaces; (2) stress redistribution; Be energy ; release; (4) progressive failure: 
= (5) initiation, development and rate of landslide movement. Significant aA 
progress towards understanding slope performance cannot be made without 


tackling such issues. Frequent discussion of these aspects of slopes will also 
“facilitate the development of new methods of analysis and radical modifications 
ones. Only when the need for such development is emphasized 
- it be possible to gain recognition and acceptance of new approaches. __ a. 
= ‘Major issues in slope analysis have been discussed in detail elsewhere (21,24). * 
New approaches for specific problems concerning natural and excavated slopes | 
have also been developed over the last few years (20,22,23). Herein attention 2 
- - will be given to specific issues relevant to the paper. Firstly, although stress ; 
we analysis and limit equilibrium procedures have been used widely, they do not 
generally achieve the same goals. Even where reliable data about soil properties, — 
_ initial ground stresses, and pore water pressures have been painstakingly collected 
after several decades of research, full agreement has rarely been obtained between 
limit equilibrium and stress analysis solutions to the same case histories. The > 
most famous example is that of delayed, first-time failures of excavated slopes 
The at author has demonstrated the scope an and flexibility ¢ of the w wedge ‘method 
of limit equilibrium analysis and postulated that the line of thrust need not 
_ be a criterion of acceptability in this case. Indeed the author has shown that 
_ one may obtain an acceptable solution with a ‘‘simplified’’ method which ignores | 
the line of thrust while unacceptable solutions may result al 


- does not necessarily provide the actual stresses in the slope and it is rigorous 
only within the context of the method of solution.”” In fact, as discussed later, 


; slope formation and history are not simulated. Moreover no consideration is 
given to the initial state of stress = the ground. Whether the slopes in author’s 
examples 1-4 were embankments, excavations or natural — would ke 


: 
a 
(simplified and rigorous) are subject to significant uncertainities. Note the author's: A - 
a c 2 ion in hi Onite On nace an 2 pntahle rivora jon 
statement, although welcome, does not go far enough in cautioning the reader. 
: 7 __ An important disadvantage of conventional limit equilibrium analysis is that 
Reader in Civ. Engrg., of Wollongong, P.O. Box 1144, Wollongong, N.S.W., 


difference to the calculated stresses on assumed slip 

As an example, consider first an. Se ae of specified height and 7 
inclination. The stresses within such a slope e depend to a significant extent — 
on the initial value of K ,, the ratio of lateral to vertical normal effective stresses. 
Since K , never enters into conventional limit equilibrium there is no hope that 

é calculations based on this concept can predict the actual stresses. 
_ Asasecond example, realistic stresses and deformations within an embankment afl 
a obtained by simulating the incremental process of construction. _ Any ha 


analysis which considers only the completed embankment gives mi isleading results 


~ Conventional limit equilibrium consistently ignores the incremental construction © 
Process and can never give the true stress distribution. 
4 _ Nevertheless, limit equilibrium methods continue to be popular because o 
q “historical reasons. The conditions under which these methods can provide reliable 
_ guidance to practitioners have been discussed comprehensively by the ~ 
OD. Natural and excavated slopes w with high initial stresses present special 
: problems which require radical modifications of conv conventional methods of analysis. — 
-™ writer’ s* ‘initial stress approach” (20,22) is a step in this direction and 
:: a recognition of the fact that the concept of conventional limit equilibrium = 
“is not ‘‘fundamental”’ to an understanding of slope behavior. 
s&s conclusion of the paper is that the shear strength on the internal shear 
a surfaces has an n effect on F,, the computed static factor of safety. It is already 
_known from the conventional wedge method that mobilization of strength a 
- interwedge boundaries influences the value of F. The assumption of the correct _ 
- mobilization factor requires judgement based on experience and appropriate 
references have been cited in Ref. 21. As long as there is uncertainty about — 
the mobilization of internal shear strengths, back-calculation of mobilized strength 


4 
4 of int 
J on a slip su surface must also be subject to uncertainty whatever method is used | 


- for calculation. In first-time slides which do not follow old shear surfaces, 


> the average strength along the slip surface may be much higher than the residual 
value and often close to the peak value. Therefore, prediction of field strength : 
parameters (residual or peak) from case histories of of first-time slips shou! ss = 
made with « caution for this reason also. 
— The author states on page 1516: “Even if the mass contained 4 wishin ~~ 
s surface is in a state of limiting equilibrium, the mass will not be able 
to move unless shear surfaces are formed within the body.” In reality there 
is no question of movement or development of shear sulioes as long as P, 
> 1 or K, > 0; limit equilibrium does not mean and should not be confused _ 
with critical equilibrium. (Note that the equations are based on the. 


arbitrary slices or wedges are merely part of the calculation technique. ca 

internal shear surfaces which develop when critical equilibrium is reached cannot 

a be predicted by conventional limit equilibrium. However, if a mechanism is 
= assumed, it may be simulated on a step by step basis as in the author’s ‘*Example 


_ The most important issue in conventional analysis is the one relating to factor — 
of safety, its definition and interpretation, and the dominant role it appears __ 
te play in analytical procedures. It must be emphasized that the assumption 
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weakness of any limit squiliten method. The stresses vary — patent: to 4 

point on a slip surface . The shear strength, which is a function of the normal > 

_ effective stress, also \ varies along the slip surface. The ratio of shear strength _ 

to ‘0 shear ‘stress, which governs the development of failure, is ‘never « constant 

_ along any slip surface. It is obvious that failure must start at specific overstressed | 

points within a specific overstressec tegion) before _complete can 
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20. Chowdhury, R. N., “A New Approach to nae Stability Studies,” > Department of 


\ Civil Engineering Report, University of Wollongong, Wollongong, Australia, 1977, ee 


21. Chowdhury, Slope Analysis, Elsevier Scientific Publishing ©o., msterdam, 
; 22. Chowdhury, R. N., “Analysis of the Vajont Slide—New Approach,” Rock Mechanics, 
_-:23. Chowdhury, R. N., ‘‘Propagation of Failure Surfaces in Natural Slopes,” Journa 
of Geophysical Research, Vol. 83, No. B12, 1978, pp. 5983-5988. 
Chowdhury, R. N., ‘Landslides as Natural Hazards—Emphasis on Alternative 
_ Mechanisms,’ a ’ Keynote Special Lecture, presented at the Jan., 1980, International — 
Conference on Engineering for Protection from Natural A.LT., held at 


by K. Chugh? 


— 
The wi writer would like to comment on. the nonuniqueness of wb for 
slope stability (1) Dee to inadequacy of assumption, leaving 


analysis of a natural slope the “duality of solution to the p 
‘This duality of solution is not due to the insufficiency of assumptions ow the 7 
as ohana to the limit equilibrium method of analysis. In this table, n is the: 
number of slices used to discretize a slide mass; and A stands for assumption. 
iy Also included in this table are the methods of analysis, their assumptions, 
: > unknowns, and number of equations. Since the number of equations in methods 
1, 2, and 3 exceed the number of unknowns, one has to discard one equation 
7 to achieve an equality between the number of equations and unknowns. Thus £ 
: the solution to the slope stability problem as obtained by these methods depends | 
& the equation eliminated. d. Alternatively, one could solve the overdetermined 
7 system in some statistical manner, but again the solution obtained could. not 
be unique. The writer’s understanding of methods 5 and 6 in Table 2, has © 
_ been that they make necessary and sufficient number of consistent assumptions 
(2m — 2) so as to achieve an equality between the number of equations and 


; = *Civ. Engr., U.S. Water and Power Resources Service, Engrg. and Research Center, 
P.O. Box 25007, Building 67, Denver Federal Center, Colo. 80225. 
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Sample Problem.—Fig. 9 a cross section of the riverbank slope located 
downstream of the Grand Coulee Dam in Washington. The stability status of 
a potential slide surface marked on Fig - 9 is of interest. For slope stability 
analysis, the geologic makeup of the site, , estimated strength values, surface 
‘Pressure caused by river ‘water, the Subsurface pore ¥ water 


‘interslice boundaries being to the axis. The ‘left face of the first 
at toe and the right face of the last slice at heel have zero force coincident 


Attempts to solve this problem using the computer code SSTABI (27), available 4 
at the U.S. Water and Power Resources Service, failed to converge to a solution. 


TABLE 2.—Details of Unknowns, Assumptions, an and | Equations for Slope Stability | 


Problem as Posed to Limit Equilibrium Method 


Norme! Forces 


o. of slices = 
for nel; 
general, the 
necessory and sufficient 
otassumptions 
quired for unique 
-2)-4n= 
2n-2 


linsertion in the 


oper 


slice thrustiexcluded 
om the criterio of 
ceptability ofo 
rica! solution. Thus 
0 assumption Needed 
o column @)ond no 
equation used for 
coumn@. 


| | 

SSTAB1 uses the limit equilibrium procedure of slices satisfying all conditions 

of static equilibrium for each slice; resultant interslice forces acting on each 
slice assumed parallel; the solution for factor of safety obtained by iteration 

requiring successive approximations for the factor of safety and interslice force = 

_ inclination in order to satisfy the equilibrium bounday conditions. The problem 

- was solved using the numerical-graphical solution procedure (25,26). Fig. 10 

7 is a plot of solution pairs (F,B); tan B = X/E. Clearly, the pair values of | 


: F, B) satisfying the force equilibrium requirement fall on two distinctly separate _ a 


curves: I and Il of f Fig. 10. Similarly, the pair values of (F,{) satisfying the 


“moment equilibrium requirement fall on two distinctly separate curves III and 
z IV of Fig. 10. These two pairs of curves intersect at least at two points. Each 
: is a anaeiin solution mation to the eee The numerical values for the solution 
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— parameters are (0.7957, 9.67005°) and 0. 652386, 15.2692238°) for the possible 


solutions | and 2, respectively. It is only for the objective of this discussion 
that the high order of Precision in calculated values of solution parameters — 
( F, B) is is reported. 5 
25. Chugh, A. K., dainty of Numerical Solutions for Slope Stability Problems,” 
_ International Journal for Numerical and Analytical Methods in Geomechanics, i in press. : 
26. Spencer, E., ‘‘Thrust Line Criterion in Embankment Stability Analyses,’’ Geotechnique, 
on, Wright, S. G., “SSTABI—A General Computer Program for Slope Stability Analyses,”’ 
_ Department of Civil Engineering, University of Texas at Austin, Austin, mF 1974. 
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wury for their interest st 
; - Chough points out that he has found two solutions of F for a slope stability bi 
_ > problem. In fact, in the limit equilibrium | method of analysis using the technique 


_ of slices and within the assumptions involved, there are n solutions of where a4 


: n is the number of slices (28). Of all these values of F, only the largest se 
gives shear forces at the base of the slices in a direction which is implicit 
in the limit equilibrium method. All other solutions are not ot 4 


The writer agrees with Choudhury that the limit equilibrium method of analysis 


is not suitable for determining the true stresses within an embankment or slope. 
One has to Tealize from the very beginning that the limit equilibrium method 


are actually existing on “surface are less. 

the ones adopted for the analysis. If that is the case, the stresses i Nene 
from the solution will keep the body in critical equilibrium (discusser’s term). — 
‘However, for noncircular and nonplanar surfaces, this is not sufficient for the 
body to be | able to move if the driving force is increased ‘Slightly. The writer 

therefore introduced the kinematics of the sliding blocks. 
_ Even though the method is applicable for general slope stability problems, — 
the writer intended the method to be used for analysis of actual slope failures 
rather than to calculate factors of safety for stable slopes. The method is also 
suitable for determining the € critical earthquake 
or a slope 


“Iterative Schemes for Slope Stability Analysis,” Numerical Methods 
: ~ “Lect., Dept. of Civ. Engrg., — Coll. of Science and Tech., Imperial College 
Road, London, SW7 2BU, England. pot 
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— 

OF analysis is applicable only for those surlaces which, given the might Circum- 

y stances, are likely to fail. Otherwise, Mohr-Coulomb failure law is not applicable. : 

{ 
| 


Geomechanics, C. S. Desai, Vol. ASCE, | New York, N.Y., 


Errata.—The | following corrections should be ma made to , the original paper: 


ge 1514, paragraph 4, line 3: Should read ‘‘are 2n 2, instead of “are 


Page 1514, paragraph 4, item 4, last sentence: Should be deleted. 


W, sin (b/ — a,) + R, cos, + S,,,sin@’, —a,—8,,,) —S,sin(¢’, ,-8,) 
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Luip WAVE PROPAGATION IN SATURATED AND NEARLY ad 


te 3:42 4 at SATURATED SANDS a 


"Discussion by ‘Paul F. Hadala‘ ‘and Robert F. Ballard, . Jr., Members, ASCE 

;.. The authors have suggested the | possibility c of using the in situ compressionwave 
"velocity from seismic tests as a way of determining if a cohesionless soil below 
the groundwater table is fully saturated or is just nearly saturated. Laboratory 
cyclic triaxial tests by Chaney (14) and others whose work is referenced in 
his paper indicates that, in nearly saturated | specimens, changes | in B factor 
from - ~0.90-~0.98 have a great | effect on the rate at which pore pressures 
and strains increase with number of load cycles in cyclic triaxial tests on clean 
sands. From the work of Black and Lee (4) and of Martin, Finn, and Seed — 
(9), this change in B factor corresponds to a change in degree of saturation 
for ~99.5% to ~99.9%. If the behavior of a sand specimen in the cyclic triaxial 
test is at least a qualitative indicator of the seismically induced deformations | 
in the field, then nearly saturated sands should be. expected to behave better 
than fully saturated sands, other things being equal. Thus, it is desirable to | 
determine if sand in the field is so nearly saturated that its B factor would — 7 

_ “March, 1980, by N. Bataan Allen, F. E. Richart, Jr., and R. D. Woods (Proc. — 


“Acting Asst. Chf., Geotechnical Lab., U.S. Army Engr. Experiment 


om of Engrs., P.O. Box 631, Vicksburg, Miss. 39180. 

*Research Geophysicist, Earthquake Engrg. and Div., "Geotechnical Lab. 

U. S. Army Engineer Waterways Experiment Station, Corps of Engrs., P.O. Box al, 
Vicksburg, Miss. 39180. 
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be 0. 97 or more. The large difference | p- -wave at 99. 5% 

2 99.9% saturation demonstrated by the authors’ data and by theory in Fig. at 
7 will not be masked by the errors normally associated with urhele, cross © 


moderate prord below a a known phreatic surface | are so ‘© nearly saturated that 
_ their B factors are 0.97 or more. This index is likely to be quite relevant in 
the evaluation of seismically- -induced strains in medium density sands 
The writers are aware of several cases where low P-wave velocities — 
measured down to some depth below the existing water table and where 
supplementary data suggest that the material should not have been — 
3 _ The first of these cases is the 40-ft (12-m) high Barkley Dam in Kentucky. 
 ThevU. s. Army Engineer Waterways Experiment Station ran consistent refraction, 
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of the adi in 1977 (pending publication). The seismic profile ‘determined from 
these tests and the groundwater level determined from piezometers and wells _ 
in the same vicinity are both shown in Fig. 10. The dom has had at least 
a few feet of water against it continuously since 1967. The depth at which 

i the seismic velocity increases from 2,300 ) fps to 6,100 fps (701 m/s-1859 m/s) 

7 ' a under the center line of the dam coincided closely with water levels measured — { 


in borings at this site prior to the construction of the dam (20). In this case, 
: materials which have been continually submerged for a decade and through 
which there has been continuous seepage still have a low seismic eg and s _ 
The second case is a high- explosive | field test site in Arizona (17) soa - 
‘in the valley of the Bill Williams River, an intermittent stream. The valley — 
fill is predominately sands and gravels. Only after local heavy rains (which 
: ‘occur rarely) or after releases from Alamo Dam, several miles upstream, does _ 
water flow on the surface. Meissner wave-front and cross hole seismic tests & 


were run at the site. in February and gave consistent results. In June, 
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refraction seismic tests were conducted in same ne area. Fig. the 
 eheneved ground-water levels and the seismic velocities determined in each case. 
While the ground-water level rose significantly in the 4-1/2 months between 
- surveys, the depth at which the P-wave velocity jumped to approx 5,000 fps 

— (1,524 m/s) rose by a a lesser amount which is consistent with the calculated 

_ effect of the pore pressure increase in a material originally about 98.5% saturated. oat 
= examples of low P-wave velocities below the ground-water level (15, 16,18) 

have been noted in the past few years. The two examples described are based 
- well- documented in situ measurements, support the theory and laboratory 


high P-wave velocity layer is cavelineaballr can sometimes be wrong. ee 
Shot ang Adis! other Greciions of 
can De actotamodated under conditions of 
Chaney, R. C., “Saturation Effects on the ‘Cyclic Strength of Sands,”” Earthquake 
ke Engineering a and Soil Dynamics, Vol. 1, ASCE, New York, N.Y., pp. 342-358. 
7 15. Curro, J. R., Jr., Horn, J. E., and Ballard, R. F., Jr., “In Situ Seismic Investigation, 
* Fort Peck Dam, Montana,’’ Miscellaneous Paper S-75-4, U.S. Army Engineer Weer 
ways Experiment Station, CE, Vicksburg, Miss., Feb., 1975. 
Ae «16. Curro, J. R., Jr., Skoglund, G. R., and Ballard, R. F., Jr., ““Geophysical Investigation, = 
= _ Prado Dam and Mentone Damsite, California,”’ Miscellaneous Paper S-75-6, U.S. 


Army Engineer District, Los Angeles, Calif., June, | 
Jackson, A. E., Jr., Ballard, R. F., Jr., and Phillips, B. R., “Geotechnical Investigation 
for MISERS BLUFF II: Results from the Subsurface Exploration Programs,”’ U.S. | 
-* Army Engineer Waterways Experiment Station, Vicksburg, Miss., Apr., 1977 (unpub- 
18. Jackson, A. E., Jr., Ballard, R. F., Jr., and Curro, J. R., SJr., , “Material Property 
on Investigation for Pre-Dice Throw I and II; Results from the Subsurface Exploration 
o eel U.S. Army Engineer Waterways Experiment Station, Vicksburg, Miss. 
4 Martin, G. R., Fiann, W. D. L., and Seed, H. B., “Effects of System Compliance ; “) 
on Liquefaction Tests,” Journal of the Geotechnical Engineering Division, ASCE, 4 
104, No. GT4, Proc. Paper 13667, Apr., 1978, pp. 463-479. 
“*Project Data for Periodic Inspection, Barkley Project, Cumberland River, Kentucky,” 
of the Army, Nashville District, Nashville, Tenn., Mar., 1976. 
Closure by N. Foster Alllen,° F. E. Richart, Jr..’F. ASCE, 


writers thank Hadala and Ballard for their ally the paper 
_ and for the valuable information contained in the discussion. Carefully obtained — > 
field data are vital to ‘effective site evaiuation, and the two cases described 
in detail and the others referenced give strong indication that incomplete saturation 
_ may occur frequently in natural or constructed deposits, even after long periods 


of low pressure submersion. As they note, this incomplete saturation may also - 


Asst. Project Engr., E. D’Appolonia Consulting Engrs., Pittsburgh, Pa. 
-7W. J. Emmons Prof. of Civ. Engrg., Univ. of Michigan, Ann Arbor, Mich. 48109. aol 
Prof. of Civ. Engrg., Univ. of Ann Arbor, Mich. 48109. 
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. quite significant i in the seismic seismic beh ior of cohesionless soils, and it indicates — — 


_ The final point made by the discussers is also most significant. The occurrence is 
- of the piezometric surface, in incompletely saturated soils, at elevations different 
_ from that identified by the P-wave velocity for saturated materials may lead 
to unpleasant surprises during construction. Interpretations of refraction seismic tae 
_ tests should always be cross-checked with water table level evaluations. o> 
_ Field verification of trends indicated by laboratory tests are most welcome, = 
and the discussion by H Hadala and Ballard represents a significant ant contribution — og 


‘Discussion by Antonio Gens,? Matthew Syn mes,° and David Hight” 

The authors have focussed attention on the rotations of principal stress 
lirections which take place in ‘many field loading or unloading situations and 
rou examined some of the effects of principal stress rotation on the senpenee 
of sands in a skillfully devised new apparatus, the directional shear cell (DSC). 
In reviewing the alternative laboratory apparatuses which have been used to 
examine the effects of rotation of principal stress directions, the authors comment 
on some of the shortcomings of the apparatus for testing hollow cylinders - : 
soil under combined axial load, torque, and internal and external radial pressures. 
In particular, they raise objections to the restrictions which they feel must 
be imposed in rotating principal stress directions under conditions of plane strain 
in the hollow cylinder apparatus (HCA); they consider that such we in 


a 


a restriction which fixes one of the no- extension directions as horizontal. 


| can only be achieved by ensuring no change in radius of the cylindrical sample, 


The writers recommend that less severe restrictions than ‘‘no change in sample 
radius” should be considered for maintaining plane strain conditions in torsional 
; tests on hollow cylindrical samples. Examine the definitions for strain in terms _ 
of cylindrical coordinates for situations independent of z and 6—radial strain: ins 

| “April, 1980, by J. Robin F. Arthur, Ken S. Chua, Treve Dunstan, and Juan I. Rodriguez 

Research Asst., Dept. of Civ. Engrg., Imperial Coll. of Science and Tech., Imperial 
College Road, London SW7 2BU, England. 
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and circumferential strain: 


= the radial, vertical and angular respectively. 

_ in a hollow cylindrical sample, the application of torque causes rotation of 
the principal stresses in concentric surfaces. The condition of plane strain then © 
implies no strain in the radial direction, i.e., €, = 0. This is satisfied by u— 
= independent of r, so that uniform radial translation of th the cylinder wall 
is the requirement “for plane strain. The condition of zero radial translation 
- satisfies the general requirement, but as the authors | point out, it is compatible 

with only one particular direction of principal stresses, 

_ By allowing uniform radial displacement, other directions of no-extension — 
lines and principal stresses can be accommodated under conditions of plane 

Sf strain, but only at the expense of a variation in €, which may cause nonuniformities — 
o in 1 the internal stresses. These variations | are a function of the geometry of 

5 the hollow cylindrical sample and the imposed stress state. By adopting a sample 


geometry with large radius and thin wall, the nonuniformities produced by radial nm 
translation can be minimized; they become more acceptable than the nonuni- — 


_ formities associated with a test in the HCA in which no radial displacement | 


With a carefully selected | sample geometry, the HCA offers considerable - 7 
advantages over the simple shear apparatus for rotation of principal stress 
directions under plane strain conditions, since, not only is the principal stress 7 
direction controlled directly, but the requirement for provision of complementary 
shear stresses is also automatically satisfied. The DSC is not subject to the 
- stress inhomogeneities which result from the curvature of samples i in the HCA, S&S 
oz it may involve greater operational difficulties. It is fair to point out that " 
in the HCA, a plane strain test can only be carried out by ensuring no change | 
in wall thickness; in other words, only average plane strain conditions can be 
tates In common with most testing apparatus, end restraints may affect 
- necessary free radial displacement, although these may be reduced in the = 


The authors additionally consider the effect of changes in the radial intermediate ve 
principal stress from that required for plane strain in the HCA. They argue a 
that the resulting strain in the radial direction simultaneously induces a circum- 
ferential strain which, together with the current value for the dilation angle, 

; Yas determines the orientation | of the major principal strain-rate direction and, bs 
in turn, ‘the major p principal stress direction. They: conclude that the major principal - a — 

stress direction is always effectively determined by the radial stress and if 
_ this direction does not coincide with that applied by stress comet at the ends, | 
an inhomogeneous strain response might be expected. 
id The writers disagree. A change in radial stress is ikea not only by | 

‘radial and circumferential strain, but also by changes in vertical and shear strain 

_ determined by the material’s constitutive law. The circumferential strain-r Tate 

p and ype angle are not sufficient, then, to define completely the strain increment 


tensor, so that the restriction on the major principal strain-rate direction inf erred 


|. 
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In fact, generalized ‘stress can be applied to hollow 
4 soil samples with no restraints on their strain response other than those resulting _ : 
from sample shape and end fixity, and both these can be reduced to acceptable _ 
‘a limits, as discussed previously. The principal stress directions are Gen: determined a 
_ by the combination of shear, vertical, and circumferential stresse sid edi q 
* In the HCA a wide range of stress states can be applied to produce rotations 
of principal stress directions so that it appears to have an important role in 
simulating, at the appropriate levels, the stress paths which are followed, e.g., _ 
by elements beneath ocean structures. In view of the dearth of knowledge 

7 concerning the effects of principal stress rotations on the behavior of soil and = 
of the different advantages and disadvantages of the HCA and DSC, it is possible y 
to regard the two types of apparatus as complementary. Research into ‘these 
effects is underway at Imperial College, London, England, using a recently 
- constructed large HCA. The work will be coordinated with that of the authors. _ 
if observations of the same effects are made in the two essentially different 4 
testing machines, this will give additional confidence in the results from both. th. 


_ The authors are of course absolutely right: rotation of principal stress directions 
| is a factor that needs to be taken account of in accumulating data and developing 
numerical models for the ‘complete behavior of soil. Indeed, , the writer suspects - 
that because of the paucity of good data in this. area, most numerical predictions 
€. g., finite element analyses) are made on the bases that rotation of principal 
stress directions does not per se have any effect on the behavior of soil; the 
y authors’ tests in which continuous rotation of axes can occur at constant principal . 
: show that th this may not be a conservative assumption. 
The authors discuss briefly the capabilities of the c< conventional axisymmetric 
- triaxial, and the true triaxial apparatus. It does not seem to the writer to be 
helpful to say that these are only able to apply sudden 90° jumps of principal 
stress direction—since in fact, paths in these tests move smoothly through stress 
_ space, albeit confined to a single plane in the conventional triaxial test. The — 
triaxial apparatus provides a freedom to to explore general principal | stress 
and strain states which is a facet of response a: as important as (though + el 
less neglected than) the need to study the effects of rotation of principal axes. 
_In the conventional triaxial apparatus, the 90° switch of principal occurs as i. 
the stress state passes through the isotropic state—and it is possible that some > 
’ effets of material response may be associated with passage through this condition a 
tea deviator stress. This restriction does not apply in the true triaxial apparatus — 
and it is not clear that something ‘special happens e every time a stress path 
crosses one of the dividing radii in the 7-plane shown in Fig. 12. 
7 yt the authors dismiss these other testing apparatus rather rapidly 
o in their understandable enthusiasm to present their fascinating directional shear 
cell. This device seems to produce surprisingly uniform deformation c nditions 
for a flexible boundary apparatus, but it is difficult to see how the same freedom = 


to apply controlled rotations of principal stress directions could be achieved wal 
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with ‘rigid own s of the to be both wenemit: 
7 shear stresses and also to change in length. As presently constructed, shear 
_ Stresses can only be applied i in one direction on each face—so that the simple 
_ shear apparatus perhaps has an advantage where the effects of m major cycles 
f _ The authors’ suggestion that the simple | shear apparatus can only apply sudden 7 
_ jumps in direction of principal stress in cyclic tests is based on an assumption — 
that the restrictions of their qualitative model of sand response (3) apply to - 
the whole of the deformation that occurs. In fact, it appears that there is a 2 
transition region (Fig. as the sand structure adjusts from a 
f developed plastic flow (which might be expected to fit in with the model — 4 
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FG. rr Plane into with Different Relative Magnitudes of 


: o{,o3,ando4; GH Shows Path of Conventional Triaxial Test Passing from Compression 
of / Arthur, et t al.) in one direction to a state of developed plastic | flow in the a 
- Opposite direction. This transitional stage is probably a manifestation of the 
: 2 sort of phenomena that are observed in the authors’ tests, both the tests - 
with a jump in principal stress directions and those with continuous rotation _ 
- - of the directions. The sand structure, which is to say the assemblage of groups - 
« particle contacts, develops as a result of some applied stress state and this 
structure will inevitably appear weak to alternative probing stress states. The 
authors do not show the volumetric strains in the continuous rotation tests—is 
= continuing straining a shearing phenomenon, with little further change in— 
volume occurring after the first few cycles? = 
The application of the model of Arthur, et al. (3) to the simple shear apparatus 


implies that a considerable « — of restraint is provided by the apparatus — 


Wood and Budhu,  (28)) that: fine sand is much less restricted by the hen 

boundaries in its ability to develop ruptures than is coarse sand. No doubt 
in practice the macroscopic behavior of a body of sand will be governed by 

the Tandom collection of and unruptured regions; a simple shear 


rere 


. 13 —Orientations of Principal Axes during Cyclic Simple Shear Test on a 


Sand; a Is Measure of Shear Strain; |, & and x are, Respectively, Orientations 
to Vertical of Major Principal Stress, Major Principal Strain Increment, and Major 


“large enough in comparison with the particle size may be able to reflect this 
macroscopic | behavior more than one which is more restrictive. tae PO ee 
However, it is certainly not the writer's intention ‘to set up as an apologist 
for the simple shear apparatus. The writer has spent “much time recently trying 
to persuade of the care in of simple shear an 


; to the deformation of the sample. It seems from radiographic observations of 
| 
— 
| | 
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father than the provi informatio and Budhu, (28)}. simple 
ps a apparatus, suitably analyzed, can provide some information about the 
behavior of soils under conditions of rotation of directions of principal axes. — 
& There may be certain practical cases where the simple shear apparatus is directly _ 
~ relevant—e. g., a gravity offshore platform supported on a thin deforming soil 
layer and subjected to lateral wave loading—but the directional shear cell 
a described by the authors offers greater freedom, greater control, , and greater 
sophistication in the hunt for the effect of the missing parameter. The writer 
4 _ forward to seeing the results of tests on other soils, and perhaps at higher | 
stress levels, from the directional shear cell. od: booW 
bas bose bores vino svad weal seodT notation 
= 28. Wood, D. M., and Budhu, M., “‘The Behaviour of Leighton Buzzard Sand in Cyclic 
_ Simple Shear Tests,” Proceedings of the International Symposium on Soils Under 
7 Cyclic and Transient Loading, G. N. Pande and O. “3 Zienkiewicz, eds., Vol. 1, 


. A. Balkema, Rotterdam, Netherlands, 1980, pp. 9-21. ee ee 
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is agreement between us on the major issue, but there are points which require 
Ss He suggests that the emphasis we placed on the restriction to sudden jainp 
“ changes: of principal stress direction of 90° in conventional and true triaxial _ 
apparatus is unhelpful. We are unrepentant because this restrictive feature has" 
been greatly neglected, and changes which appear smooth in the context of < 
= Fig. 12 may not be in terms of stress-strain behavior. We do of course 
: agree that there need not be any sudden change in soil behavior when such > 
aninterchange of principal stresses occurs because this behavioral depends 
the deposition and stress-strain history of the soil. 
We have perhaps placed too much emphasis on the uniformity strain 
_ distribution achieved in the directional shear cell. The data we have suggest 
" that from 15 mm in from the sample faces the strain distribution is as uniform _ 
| that in the middle third of a triaxial compression sample ‘sheared within © 
enlarged lubricated end plattens. Those who have studied strains in such — a 
will realize that appreciable variations occur even under these rather ideal 
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aximum rotation of major or principal : stress direction — 
can be achieved in the simple shear apparatus (SSA) (DSC 90°; SSA estimated 
_ for dense sand at 110°). However we maintain our simplifying assertion that 
the strain restraint imposed in cyclic loading by the SSA may be regarded ~ 
as leading to sudden jump rotations. This is based on the very simple assumption 
of coincidence of axes of stress and strain rate, rather the somewhat speculative. 
— model of the “first writer, et al. (3). Wood’s suggestion of a transition zone 
is quite acceptable to us as it appears to imply similar phenomena to those i 
_we have observed in our own jump rotation tests carried out in the directional 
qa ‘Wood raises the important question of volumetric strain in our continuous — 
rotation tests. These tests have only beeen carried out on dense sand and, 
- in tests in which shear strain continued, continuing dilation was observed. The ” 
_ tests were stopped at the limits imposed by the cell so the end of dilation — 
~ could not be observed. In practical terms in the sea-bed, this dilation els 


be ‘waneinted into a a negative- pore water pressure of uncertain and 


Finally it was ‘not our to suggest that the simple shear 


is always inappropriate for cyclic testing but ‘merely to Gow ehention to a 

Our colleagues at Imperial College point out that, in suggesting 

oe was an essential Strain restraint in hollow cylinder testing, we ignored 
the possibility of varying the radius” of the soil cylinder during shear strain 
as well as the cylinder thickness. If the radius could be varied freely there 
4 would be true stress control with no strain restraint and only a residual difficulty 
in attempting plane strain tests. However, the requirement to apply shear at 
- the ends of the cylinder is incompatible with freedom for the cylinder to change - 
ee at the ends. This restraint will inevitably lead to inhomogeneous s ner 


“away from the ends. The Imperial College design i is ingenious in using a i 
- evlindr at the ends to steer most of the deformation away from the ends 
_ to the zone of greatest freedom from the restraint. Thus, the strain inhomogeneity | 
_ is “‘guided’’ in a very desirable way, , but there may well be problems at larger — 
strains when the radius in the center of the cylinder has increased relative 
to the ends. As we suggested in our paper, there is a need in such potentially — 
inhomogeneous situations for detailed strain measurements over the deforming 
"portion of the sample. At least some measurements at different radii within 
_ Naturally, we look forward to coordinating our research and gaining the 
q “advantages and chall challenges of data two different 
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DISCUSSION 


Discussion Erwin Kessler? 


es ve 
higher than those given by Frank and Varaskin e. we = 
“ton; drop height = 20 m; distance = 40 m; ont particle velocity = 


(corresponding toa scaledenergy factorof 11), 


The following measurements by ourselves at the base of a building, 500 - 


_ from: the densification | site and “tl weight = = 12 t ton; drop height = 18 m; 


factor ‘of 0.74), also yield a , point in the fe can of Fig. 5, that is significantly 


_ below the indicated range of particle velocities; 
and Varaskin, “Verdichten von Biden durch dynamische ‘Einwirking 
Baumaschine + Bautechnik, Heft 


er 170 by Robert G. Lukas,’ M . ASCE 


writer appreciates the comments of E. measured particle 
—- velocities that are much lower than would be predicted by Fig. 5. Unfortunately, 
: Kessler does not describe the type. of soil or fill deposit that was. densified. 
id The measured particle velocity is not only a function of the scaled energy _ 
- factor but the type and relative density of the soil deposit being improved. saute 
- _ Since the original article was published, additional data have been obtained 
: relating measured particle velocity to the scaled energy factor for a very stiff i 
_ sandy clay fill and loose decomposed refuse. These data are reproduced on _ a 
Fig. In addition, the relationship between particle velocity. and the scaled 
energy factor for a loose sand as reported by Leonards (4) is also shown. | 
_ The data reported by Kessler plot below the loose decomposed garbage fill 
‘es, indication that the deposit that was densified must have been extremely 
- loose in order to obtain such low particle velocity readings. One of the readings — 
taken by Kessler may not be appropriate since it was taken 500 m trom the 


“April, 1980, by Robert G. Lukas (Proc. Paper 15376). 
— *Civ. Engr., Heierli Ag, Postfach 8033, 8006 Zurich, Switzerland. 
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a densification site and the measured ed particle velocity | may have been too hie 
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441, line 3: Should read “and. decrease instead of ‘ “and 
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DISCUSSION 


The — correction should be made to the original paper: 


Page 1363, Fig. 9, equation for. M...: 
= —w/8 L’’L (Center Lift) instead of M,,,, = 9/8 (Edge Lift) 
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Original papers should be submitted in triplicate to the Manager of juideibe ond Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. Authors must indicate the a es 
Technical Division or Council, Technical Committee, Subcommittee, and Task Committee (if any) - 
to which the paper should be referred. Those who are planning to submit material will al a 
the review and publication procedures by complying with the following basic ee i 
1. Titles must have a not 50 and af Pe 


i 2. The € manuscript oo a original ribbon copy an and two duplicate copies) should be double- -spaced 
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subject cannot be adequately presented within the 10,000-word limit, the paper should be accompanied 7 
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published by the American National Standards Institute and to the Authors’ Guide to the Publications 


ee Figures should be drawn in black ink, at a size that, with a 50% reduction, would have 
a published width in the Journals of from 3 in. (76 mm) to 4-1/2 in. (110 mm). The lettering | 
an be legible at the reduced size. Photographs should be submitted as glossy prints. Explanations — 
descriptions must be placed in text rather than within the figure. 


8. Tables should be typed (an original ribbon copy and two duplicates) on one side of 8- 1/2-in. 


a (220-mm) by 1 1-in. (280-mm) paper. An explanation of each table must appear in the text. ees hs on 
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end of the paper, or preceding the Appendix.—Notation, as an Appendix.—References. — iy fi, eo 
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